Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,109 +C 1.6% intergenic (+48/‑123) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1091.→C 1.6% 6.8 ‑1.8 ‑0.0359intergenic (+48/‑123)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (4/1):  ref base (150/203):  total (155/204)
Fisher's exact test strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.31e-01

Only 100 of 418 total aligned reads displayed.

AAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.T.GGGGGCTCTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGG  >  NC_014034/3609011‑3609208
                                                                                                       |                                                                                                    
 aGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.g                                                                                                     <  1:3431926/101‑1 (MQ=39)
   gTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.ggg                                                                                                   >  2:2898043/1‑100 (MQ=39)
    tctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTAAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggc                                                                                                  <  1:5338962/101‑3 (MQ=38)
     ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggca                                                                                                 >  1:657471/1‑98 (MQ=38)
     ctcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggca                                                                                                 >  1:6771019/1‑98 (MQ=38)
      tcGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTT.g.gggcaa                                                                                                <  2:5536852/101‑5 (MQ=38)
       cGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGGCCTGCCCCCCGCCTCGCTCCCCCCCGGGTAGTT.gtggggag                                                                                                >  2:1640834/1‑94 (MQ=255)
          ccagccagCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTa                                                                                             <  1:3757739/101‑1 (MQ=38)
           cagccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCCCTGCCCCCCCCCCCGCTCCCCCCCGGGGATTT.t.gggcgcagtga                                                                                           >  2:5677019/1‑92 (MQ=255)
            ccccagcGCGTCCCCACCCGACCTCGCGTGCCGAACTGATCCCGGAGTT.GGGGGCTCTGCCCCCCCCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaa                                                                                           <  1:6931588/99‑1 (MQ=16)
             gccagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTaaaa                                                                                          <  2:1302728/101‑1 (MQ=38)
               cagcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGc                                                                                        >  2:4097309/1‑101 (MQ=37)
                agcGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCa                                                                                       >  2:4935777/1‑101 (MQ=37)
                 gcgcGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCAt                                                                                      <  1:322771/101‑1 (MQ=38)
                      tCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGccc                                                                                 >  1:5650200/1‑101 (MQ=37)
                       ccAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCt                                                                                >  2:6428781/1‑101 (MQ=37)
                        cAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCtt                                                                               >  2:1880511/1‑101 (MQ=37)
                         aGACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCttt                                                                              <  2:122037/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                             <  2:6723260/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                             <  2:538812/101‑1 (MQ=38)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                             >  1:5933999/1‑101 (MQ=37)
                          gACCGAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTa                                                                             <  2:5133349/101‑1 (MQ=38)
                              gAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                                                         >  2:6912648/1‑101 (MQ=37)
                              gAACCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCtt                                                                         <  2:2733003/101‑1 (MQ=38)
                                aCCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                       <  2:1733129/101‑1 (MQ=38)
                                aCCGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGc                                                                       <  1:1382435/101‑1 (MQ=38)
                                 ccGCGCGTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGcc                                                                      <  2:3603023/101‑1 (MQ=38)
                                    cgcgTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCa                                                                   <  2:3773739/101‑1 (MQ=38)
                                      cgTGCTGAACTGATACGGGAG.TGGGGGGCTTTGCCCCCCGCTTCGCTCCCCCCCGGATATTTGG.GGCAAGGGAAAAGCGTTGCCCTTTATCTTGCCTCTGa                                                                 >  1:3376358/1‑101 (MQ=255)
                                      cgTGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGa                                                                 >  1:1496114/1‑101 (MQ=37)
                                        tGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                               <  2:5335168/101‑1 (MQ=38)
                                        tGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                               >  1:5206920/1‑101 (MQ=37)
                                        tGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                               <  2:5769768/101‑1 (MQ=38)
                                        tGCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                               >  1:2527993/1‑101 (MQ=37)
                                         gCTGAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                              <  2:5489854/101‑1 (MQ=38)
                                            gAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAANCATTGCCCTTTATCTTGCCTCAGATATccc                                                           >  1:1782995/1‑101 (MQ=37)
                                            gAACTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                           <  2:4203396/101‑1 (MQ=37)
                                             aaCTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCg                                                          >  1:1883687/1‑101 (MQ=37)
                                               cTGATACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg                                                        >  2:2041333/1‑101 (MQ=37)
                                                  aTACGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGt                                                     <  2:2821528/101‑1 (MQ=38)
                                                    aCGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                   <  1:4409267/101‑1 (MQ=37)
                                                    aCGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                   <  2:4368958/101‑1 (MQ=37)
                                                    aCGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                   >  2:6527791/1‑101 (MQ=25)
                                                    aCGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCCGGATATTTGG.GGCAAGGTAAACGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                   >  1:299732/1‑101 (MQ=37)
                                                     cGGGAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                  >  1:5584373/1‑101 (MQ=25)
                                                      gggAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                                 <  2:2372058/101‑1 (MQ=37)
                                                      gggAG.TGGGGGGCTCTGCCCCCCGCTTCGCTCCACCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAt                                                 >  1:5021957/1‑101 (MQ=25)
                                                           gtggggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                             <  1:3561066/98‑1 (MQ=37)
                                                           gtggggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                             <  2:6813769/98‑1 (MQ=37)
                                                            tggggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                            >  1:947206/3‑101 (MQ=37)
                                                             ggggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                                           >  2:6351769/2‑101 (MQ=37)
                                                              gggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg                                          >  1:4826564/1‑101 (MQ=37)
                                                              gggggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcg                                          <  1:3138438/101‑1 (MQ=37)
                                                                gggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCa                                        <  2:84358/101‑1 (MQ=38)
                                                                 ggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                       <  2:3953384/101‑1 (MQ=38)
                                                                 ggCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                       <  2:33441/101‑1 (MQ=38)
                                                                  gCTCTGCCCCCCGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGc                                      <  1:2885509/101‑1 (MQ=37)
                                                                     ctgccccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                   >  1:2001603/5‑101 (MQ=255)
                                                                       gccccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                 >  1:5702382/3‑101 (MQ=255)
                                                                       gccccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                 <  2:5924319/99‑1 (MQ=255)
                                                                       gccccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaaa                                 <  1:4273671/99‑1 (MQ=255)
                                                                         cccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                               >  2:2237457/1‑101 (MQ=255)
                                                                          ccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGTTGAATTTTGGCGCAGCCAAAaga.g                             >  1:735810/1‑101 (MQ=255)
                                                                          ccccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                             <  1:271870/101‑1 (MQ=255)
                                                                          ccccGCTTCGCTCCCCCCGGGATATTT.GAGGCGAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                             >  1:6272088/1‑101 (MQ=255)
                                                                            ccGCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggg                           <  1:5970334/101‑1 (MQ=255)
                                                                              gCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.ggggg                         >  2:3435436/1‑100 (MQ=255)
                                                                               cTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggc                        >  1:3467647/1‑99 (MQ=255)
                                                                                ttCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggca                       >  2:4580443/1‑98 (MQ=255)
                                                                                  cGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcagc                     <  2:2577181/101‑6 (MQ=255)
                                                                                   gCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                     >  2:2986459/1‑101 (MQ=25)
                                                                                    cTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgc                    >  2:6045013/1‑101 (MQ=21)
                                                                                      ccccgcGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                  <  1:2255014/101‑1 (MQ=21)
                                                                                      ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgcac                  >  2:660182/1‑99 (MQ=21)
                                                                                      ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                  >  1:4672587/1‑101 (MQ=21)
                                                                                      ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                  >  1:6128991/1‑101 (MQ=21)
                                                                                       cccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcccc                 >  1:4097975/1‑101 (MQ=21)
                                                                                        ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                <  1:126944/101‑1 (MQ=21)
                                                                                        ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                <  2:208535/101‑1 (MQ=21)
                                                                                         cccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCt               <  2:6247335/101‑1 (MQ=21)
                                                                                         cccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCt               <  2:21756/101‑1 (MQ=21)
                                                                                         cccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCt               <  2:1454125/101‑1 (MQ=21)
                                                                                          ccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTc              <  2:4826321/101‑1 (MQ=21)
                                                                                           cGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcc             >  1:6372303/1‑101 (MQ=21)
                                                                                            gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc            >  1:3932691/1‑101 (MQ=21)
                                                                                            gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc            >  2:6996182/1‑101 (MQ=21)
                                                                                            gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc            >  2:2681984/1‑101 (MQ=21)
                                                                                            gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc            <  1:7082524/101‑1 (MQ=21)
                                                                                              gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc          <  1:4141356/101‑1 (MQ=18)
                                                                                              gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc          <  1:3123230/101‑1 (MQ=18)
                                                                                               atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa         >  1:1054010/1‑101 (MQ=21)
                                                                                               atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa         <  1:2363994/101‑1 (MQ=18)
                                                                                               atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa         >  1:5477109/1‑101 (MQ=18)
                                                                                                tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc        >  2:4662338/1‑101 (MQ=21)
                                                                                                tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc        <  2:3132588/101‑1 (MQ=18)
                                                                                                 atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg       <  1:1669874/96‑1 (MQ=255)
                                                                                                 atttgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACg       >  1:1127044/6‑101 (MQ=255)
                                                                                                  tttgg.ggCAAGGTAACAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACgg      <  1:3278816/97‑1 (MQ=18)
                                                                                                    tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc    >  2:6193548/3‑101 (MQ=255)
                                                                                                    tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc    >  1:6505339/3‑101 (MQ=255)
                                                                                                       |                                                                                                    
AAGGTCTCGGCCAGCCAGCGCGTCCAGACCGAACCGCGCGTGCTGAACTGATACGGGAG.T.GGGGGCTCTG.CCCCCGCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGG  >  NC_014034/3609011‑3609208

Base quality scores:  ATCG  < 3 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 35 ≤  ATCG  < 41 ≤  ATCG