Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 3,609,145 G→A 1.3% intergenic (+84/‑87) clpB → / → dinB chaperone ClpB/DNA polymerase IV

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140343,609,1450G→A 1.3% 4.0 ‑0.4 ‑0.0437intergenic (+84/‑87)clpB/dinBchaperone ClpB/DNA polymerase IV
Reads supporting (aligned to +/- strand):  new base (4/2):  ref base (238/193):  total (242/195)
Fisher's exact test strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 456 total aligned reads displayed.

GCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCGCTTT  >  NC_014034/3609048‑3609244
                                                                                                       |                                                                                                    
gcgTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAg                                                                                                      <  1:5960285/101‑1 (MQ=38)
  gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                                    <  2:3907768/101‑1 (MQ=38)
  gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                                    <  1:6916960/101‑1 (MQ=38)
  gTGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGat                                                                                                    >  2:1356175/1‑101 (MQ=37)
   tGCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGata                                                                                                   <  2:2164628/101‑1 (MQ=38)
    gCTGAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGatat                                                                                                  >  1:7017510/1‑101 (MQ=37)
       gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                                               >  2:5377706/1‑101 (MQ=37)
       gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                                               >  2:5651371/1‑101 (MQ=37)
       gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                                               >  2:1417602/1‑101 (MQ=37)
       gAACTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATccc                                                                                               >  2:6163919/1‑101 (MQ=37)
          cTGATACGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCggg                                                                                            >  1:6729935/1‑101 (MQ=37)
               aCGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                                                       <  1:5536689/101‑1 (MQ=37)
               aCGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGa                                                                                       <  2:4368958/101‑1 (MQ=37)
                cGGGAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGaa                                                                                      >  1:5584373/1‑101 (MQ=25)
                 gggAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGGGAAt                                                                                     >  2:3278700/1‑101 (MQ=25)
                  ggAGTGGGGGGCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAAtt                                                                                    >  2:2825769/1‑101 (MQ=25)
                     gtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                                                                 <  1:3561066/98‑1 (MQ=37)
                     gtggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTg                                                                                 <  2:6813769/98‑1 (MQ=37)
                      tggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                                                                >  1:947206/3‑101 (MQ=37)
                      tggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTgg                                                                                <  2:6860260/99‑1 (MQ=38)
                       ggggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgc                                                                               >  1:417483/2‑101 (MQ=37)
                         ggggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGgcgc                                                                             >  1:5884698/1‑101 (MQ=37)
                           ggCTCTGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAg                                                                           >  2:3736803/1‑101 (MQ=37)
                           ggCTCTG.CCCCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaag                                                                            >  2:4173107/1‑93 (MQ=2)
                            gCTCTG.CCCCCCGGCTTCGCTCCCCCCGGGATATTT.GCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGTGAAttggccaaagg                                                                           >  2:3115225/1‑92 (MQ=2)
                             ctctGCCCCCC..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGcc                                                                         >  2:2652525/1‑101 (MQ=37)
                              tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                                        <  1:715498/96‑1 (MQ=255)
                              tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                                        >  1:609605/6‑101 (MQ=255)
                              tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                                        >  1:1794536/6‑101 (MQ=255)
                              tctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAAATATCCCGGGGGTGAATTTTGGCGCAGCCa                                                                        <  1:3995438/96‑1 (MQ=25)
                               ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                                       >  1:348309/5‑101 (MQ=255)
                               ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                                       <  2:5933695/97‑1 (MQ=255)
                               ctgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaa                                                                       >  1:326681/5‑101 (MQ=255)
                                tgcccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCaaa                                                                      <  1:3159304/98‑1 (MQ=255)
                                  cccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                                    <  2:3998698/100‑1 (MQ=255)
                                  cccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGAATTGCCCATTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAag                                                                    <  1:2800169/100‑1 (MQ=37)
                                   ccccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga                                                                   <  1:2480788/101‑1 (MQ=255)
                                    cccc..GCTTCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAaga.g                                                                 <  1:271870/101‑1 (MQ=255)
                                            ttCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggca                                                           <  1:2337481/101‑4 (MQ=255)
                                             tCGCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.gggggcag                                                          >  2:1072397/1‑97 (MQ=255)
                                               gCTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcg                                                         <  1:482084/101‑1 (MQ=25)
                                                cTCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgc                                                        <  1:3385870/101‑1 (MQ=21)
                                                 tCCCCCCGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcc                                                       >  1:5421856/1‑101 (MQ=21)
                                                  ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAgcgcac                                                      >  2:660182/1‑99 (MQ=21)
                                                  ccccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccc                                                      >  1:6245099/1‑101 (MQ=21)
                                                   cccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGcccc                                                     <  1:191119/101‑1 (MQ=21)
                                                    ccccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGccccc                                                    <  1:126944/101‑1 (MQ=21)
                                                     cccGGGATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCt                                                   <  2:6247335/101‑1 (MQ=21)
                                                        gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc                                                >  1:1520467/1‑101 (MQ=21)
                                                        gggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccc                                                >  2:4811669/1‑101 (MQ=21)
                                                         ggATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcccc                                               >  1:753979/1‑101 (MQ=21)
                                                          gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc                                              <  1:3123230/101‑1 (MQ=18)
                                                          gATATTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTccccc                                              <  1:568616/101‑1 (MQ=18)
                                                           atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTcccccc                                             >  2:6950741/1‑100 (MQ=18)
                                                           atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                                             >  2:2061872/1‑101 (MQ=21)
                                                           atatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCa                                             >  1:1054010/1‑101 (MQ=21)
                                                            tatTTGG.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCAc                                            <  2:5455841/101‑1 (MQ=18)
                                                                tgg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGcc                                        >  1:6505339/3‑101 (MQ=255)
                                                                 gg.ggCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCg                                       >  1:5748314/2‑101 (MQ=255)
                                                                      cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                                   >  2:880871/1‑101 (MQ=255)
                                                                      cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                                   <  1:1677278/101‑1 (MQ=255)
                                                                      cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                                   <  2:1782931/101‑1 (MQ=255)
                                                                      cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                                   >  1:4303931/1‑101 (MQ=255)
                                                                      cAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTc                                   >  1:4992827/1‑101 (MQ=255)
                                                                        aGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCt                                 <  1:6638669/101‑1 (MQ=255)
                                                                         ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg                                >  1:1456380/1‑101 (MQ=255)
                                                                         ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg                                >  2:2710303/1‑101 (MQ=255)
                                                                         ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTg                                >  1:4646291/1‑101 (MQ=255)
                                                                         ggTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGCCCtc                                >  1:2215499/1‑100 (MQ=255)
                                                                          gTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGc                               >  2:5737793/1‑101 (MQ=255)
                                                                           tAAAATCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCg                              <  2:2263175/101‑1 (MQ=255)
                                                                            aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGt                             >  1:2924968/1‑101 (MQ=255)
                                                                            aaaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGt                             >  1:3877919/1‑101 (MQ=255)
                                                                             aaaGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGtt                            >  1:2589692/1‑101 (MQ=255)
                                                                               atcATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAg                          <  2:2229935/99‑1 (MQ=255)
                                                                                 cATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTCCGTTAGat                        >  2:5835734/1‑100 (MQ=255)
                                                                                 cATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGGCCTGCGTTAGaa                        <  1:6352093/101‑1 (MQ=255)
                                                                                  aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt                       <  2:3352453/101‑1 (MQ=255)
                                                                                  aTTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAAt                       <  2:6496053/101‑1 (MQ=255)
                                                                                    tGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg                     <  1:1754669/101‑1 (MQ=255)
                                                                                    tGCCCTTTATCTTGCCTCAGATAGCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTg                     <  2:7033611/101‑1 (MQ=255)
                                                                                     gCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGa                    >  2:2758044/1‑101 (MQ=255)
                                                                                     gCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGa                    >  1:2139708/1‑101 (MQ=255)
                                                                                       ccTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATc                  >  2:4754939/1‑101 (MQ=255)
                                                                                       ccTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATc                  <  2:5046025/101‑1 (MQ=255)
                                                                                        cTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCa                 >  1:3043416/1‑101 (MQ=255)
                                                                                         tttATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCAc                >  1:113358/1‑101 (MQ=255)
                                                                                          ttATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCAcc               <  1:3736953/101‑1 (MQ=255)
                                                                                          ttATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCAcc               >  2:1936485/1‑101 (MQ=255)
                                                                                           tATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCa              <  2:6385867/101‑1 (MQ=255)
                                                                                               ttGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGcc          >  2:4539171/1‑101 (MQ=255)
                                                                                               ttGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGcc          <  2:3453394/101‑1 (MQ=255)
                                                                                                 gCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCcg        >  1:455225/1‑101 (MQ=255)
                                                                                                 gCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCcg        >  1:1215490/1‑101 (MQ=255)
                                                                                                  ccTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCcgc       >  2:1562243/1‑101 (MQ=255)
                                                                                                    tCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCgcgc     >  1:5199490/1‑101 (MQ=255)
                                                                                                     cAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCCCt    >  1:5348826/1‑101 (MQ=255)
                                                                                                      aGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCGCtt   <  2:2527899/101‑1 (MQ=255)
                                                                                                      aGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCGCtt   <  1:4105815/101‑1 (MQ=255)
                                                                                                       gATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGAGGGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCGCttt  <  2:2398314/101‑1 (MQ=255)
                                                                                                       |                                                                                                    
GCGTGCTGAACTGATACGGGAGT.GGGGGCTCTG.CCCCC..GCTTCGCTCCCCCCGGGATATTT.G.GGCAAGGTAAAAGCATTGCCCTTTATCTTGCCTCAGATATCCCGGGGGTGAATTTTGGCGCAGCCAAAAGA.GGGGCAGCGCCCCCTCCCCCACGGCCGGGTCCTGCGTTAGAATTGATCACCATGCCCGCGCTTT  >  NC_014034/3609048‑3609244

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