Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,214,790 +A 5.9% intergenic (‑46/+12) dctQ3 ← / ← dctP3 TRAP C4‑dicarboxylate transport system permease subunit DctQ/TRAP C4‑dicarboxylate transport system permease subunit DctP

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,214,79010.→A 5.9% 2.3 0.0 ‑0.034intergenic (‑46/+12)dctQ3/dctP3TRAP C4‑dicarboxylate transport system permease subunit DctQ/TRAP C4‑dicarboxylate transport system permease subunit DctP
Reads supporting (aligned to +/- strand):  new base (2/0):  ref base (31/1):  total (33/1)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

TCGCGATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAACCTC  >  NC_014034/3214690‑3214888
                                                                                                               |                                                                                                           
tcgcGATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGagcc                                                                                                              <  1:3891883/101‑1 (MQ=255)
 cgcgATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGcgcc..........g                                                                                                   >  2:981178/1‑101 (MQ=255)
 cgcgATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGcgcc..........g                                                                                                   >  2:3446218/1‑101 (MQ=255)
 cgcgATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGcgcc..........g                                                                                                   >  2:2527735/1‑101 (MQ=255)
 cgcgATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGcgcc..........g                                                                                                   >  1:968369/1‑101 (MQ=255)
   tgATGAGCGCGGCAATCAGCCCCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGcgcc..........gtc                                                                                                 <  1:2910237/100‑3 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGggcc..........gtca                                                                                                >  2:2910024/1‑98 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtcc                                                                                                <  1:3755856/101‑1 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtcc                                                                                                <  2:4173278/101‑1 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtcc                                                                                                <  2:2874870/101‑1 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtcc                                                                                                <  2:1862244/101‑1 (MQ=255)
    gATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtcc                                                                                                <  2:1920485/101‑1 (MQ=255)
     aagAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccg                                                                                               >  2:1496258/3‑101 (MQ=255)
     aTGAGGGCGCCAATCAGCCCCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccg                                                                                               <  2:2063117/101‑1 (MQ=255)
     aTGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccg                                                                                               >  1:2351580/1‑101 (MQ=255)
     aTGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccg                                                                                               <  1:4269174/101‑1 (MQ=255)
      tGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccgg                                                                                              <  2:1230578/101‑2 (MQ=255)
      tGAGCGCGGCAATCAGCACCTCGTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccgg                                                                                              <  1:4576686/101‑2 (MQ=255)
       gAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggc                                                                                             <  1:3994664/101‑3 (MQ=255)
       gAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggc                                                                                             >  1:1902890/1‑99 (MQ=255)
       gAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggc                                                                                             <  2:734447/101‑3 (MQ=255)
       gAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggc                                                                                             >  2:4245261/1‑99 (MQ=255)
        aGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........GGGCGcc..........gtccggcc                                                                                            >  1:2706266/1‑98 (MQ=255)
          cgcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgcc                                                                                          <  2:3110074/101‑6 (MQ=255)
          cgcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgcc                                                                                          <  1:4693784/101‑6 (MQ=255)
          cgcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgcc                                                                                          <  1:2075188/101‑6 (MQ=255)
          cgcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgcc                                                                                          <  2:2628321/101‑6 (MQ=255)
          cgcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........ggccggcgcc                                                                                          >  1:4817398/1‑96 (MQ=255)
           gcgGCAATCAGCCCCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccg                                                                                         <  2:4213547/101‑7 (MQ=255)
           gcgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGGGcc..........gtccggcgccc                                                                                         >  2:3173918/1‑95 (MQ=255)
            cgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgc                                                                                        >  2:1185227/1‑94 (MQ=255)
            cgGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGACGTGGAAGGAAAGGGAAAGGG..........GGGAGac..........gtccggcgccgc                                                                                        >  2:2692723/1‑94 (MQ=255)
              gCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgccc                                                                                      <  1:3503875/101‑10 (MQ=255)
              gCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgccc                                                                                      >  1:2224238/1‑92 (MQ=255)
              gCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgccc                                                                                      <  1:573178/101‑10 (MQ=255)
              gCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgccc                                                                                      <  1:2904766/101‑10 (MQ=255)
               cAATCAGCACCTCTTCGGCCCGGTCGAGGATGCTCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgcccg                                                                                     <  1:4324324/101‑11 (MQ=255)
               cAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGcc..........gtccggcgccgcccg                                                                                     >  2:1237620/1‑91 (MQ=255)
                                            aTGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGa                                                                  >  1:538671/1‑101 (MQ=255)
                                                 cagcaTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGGcggcg                                                             >  2:1066504/1‑101 (MQ=255)
                                                  agcaTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGAcggcgg                                                            >  1:4034047/1‑101 (MQ=255)
                                                  agcaTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGCCGTCCGGCGCAGCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGAccgcgg                                                            >  1:1164736/1‑101 (MQ=255)
                                                   gcaTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACggcggc                                                           >  1:3579014/1‑101 (MQ=255)
                                                    caTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGAcggnnngn                                                          >  1:3583035/1‑96 (MQ=255)
                                                    caTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGgcggcg                                                          >  2:4249343/1‑101 (MQ=255)
                                                    caTGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGgcggcg                                                          >  1:2935588/1‑101 (MQ=255)
                                                      tGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGAt                                                        >  2:1330855/1‑101 (MQ=255)
                                                          ggcCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGGTGGtt                                                    >  2:982609/1‑101 (MQ=255)
                                                          ggcCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGtt                                                    >  2:3326722/1‑101 (MQ=255)
                                                          ggcCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGtt                                                    >  2:1107327/1‑101 (MQ=255)
                                                          ggcCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGCCGTCCGGCGCAGCCCGCCGGTCTTATTCCGCCGTCGCGGACTTGACGGCGGCGGTgggt                                                    >  1:2754679/1‑99 (MQ=255)
                                                           gcCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGCCGCGGCCTTGACGGCGGCGATGGgtt                                                   >  1:1060573/1‑101 (MQ=255)
                                                            cccccGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTc                                                  >  2:4635117/1‑101 (MQ=255)
                                                             ccccGCTTAGGTCGTGGAAGGAAAGGGAAAGGGCGGCGCCGTCCGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGGTTCg                                                 >  2:4272794/1‑101 (MQ=255)
                                                                                              cggcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGTAcc                >  1:4001930/11‑101 (MQ=255)
                                                                                              cggcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCAcc                >  2:158113/11‑101 (MQ=255)
                                                                                               ggcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCg               >  2:2496459/10‑101 (MQ=255)
                                                                                               ggcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCg               >  2:1417779/10‑101 (MQ=255)
                                                                                               ggcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCg               >  1:810700/10‑101 (MQ=255)
                                                                                               gcccccgcccggccCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCg               <  2:4033744/92‑1 (MQ=255)
                                                                                                gcgccgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCgg              <  1:1417906/93‑1 (MQ=255)
                                                                                                 ccccgccgggcgCC..........CCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCggg             <  1:864727/93‑1 (MQ=255)
                                                                                                    cgtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTc          >  2:4818285/5‑101 (MQ=255)
                                                                                                    cgtccgccgcc..........ccCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTc          <  2:457735/94‑1 (MQ=255)
                                                                                                     gtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGGGGTGGTTTCGGCGCCGATCCCCCCCGGCATTTTGTCCTTCACCCGGGGCa         >  2:4300181/4‑101 (MQ=255)
                                                                                                     gtccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCa         >  2:3403782/4‑101 (MQ=255)
                                                                                                      tgctgcgcC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAg        >  2:2656707/5‑101 (MQ=255)
                                                                                                      tccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAg        >  2:2071995/3‑101 (MQ=255)
                                                                                                      tccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAg        >  1:4861566/3‑101 (MQ=255)
                                                                                                      tccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAg        >  2:2237100/3‑101 (MQ=255)
                                                                                                       cggccgcc..........ccCCGGAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTccat       <  2:2598028/97‑4 (MQ=255)
                                                                                                       ccggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGa       >  2:3138709/2‑101 (MQ=255)
                                                                                                        cggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGaa      >  1:309946/1‑101 (MQ=255)
                                                                                                        cggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGaa      >  2:1229403/1‑101 (MQ=255)
                                                                                                        cggcgCC..........CCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGaa      <  2:2357909/101‑1 (MQ=255)
                                                                                                         ggcgCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAAc     >  1:4996265/1‑101 (MQ=255)
                                                                                                         gccgcc..........ccCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAAc     <  2:3588565/99‑1 (MQ=255)
                                                                                                          gcgcC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGCCAGAAcc    >  2:2280543/1‑101 (MQ=255)
                                                                                                           cgcc..........gcCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAACct   <  1:1229520/101‑1 (MQ=255)
                                                                                                            gcc..........gccCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCTTCATCCCCCGGGGTCAGAACctc  >  2:1464178/1‑101 (MQ=255)
                                                                                                            gcc..........gccCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAACctc  <  1:4818612/101‑1 (MQ=255)
                                                                                                               |                                                                                                           
TCGCGATGAGCGCGGCAATCAGCACCTCTTCGGCCCGGTCGAGGATGCGCAGCATGGCGGCCCCCGCTTAGGTCGTGGAAGGAAAGGGAAAGGG..........CGGCGCC..........GCCCGCAGGTCTTATTCCGCCGTCGCGGCCTTGACGGCGGCGATGGTTTCGGCGCCGATCCGCTCGGCCATTTCGTCATGCACCGGGGTCAGAACCTC  >  NC_014034/3214690‑3214888

Base quality scores:  ATCG  < 3 ≤  ATCG  < 18 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG