Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 2,185,930 +C 87.6% intergenic (+203/‑302) mtlK → / → RCAP_rcc02028 mannitol 2‑dehydrogenase/FeoA family protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140342,185,9301.→C87.6% 30.5 ‑1.6 ‑0.889intergenic (+203/‑302)mtlK/RCAP_rcc02028mannitol 2‑dehydrogenase/FeoA family protein
Reads supporting (aligned to +/- strand):  new base (32/46):  ref base (2/9):  total (34/55)
Fisher's exact test strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.61e-01

Only 100 of 106 total aligned reads displayed.

GCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCC..GCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCAC  >  NC_014034/2185830‑2186027
                                                                                                       |                                                                                                      
gCGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCCAAGGCCTTG.TTTTGATGGGACAAGAAAAGCCGAA.CCG.c                                                                                                       >  1:31410/1‑101 (MQ=25)
 cGGGGGCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGGCTTG.TTTTGATCGCACAAAAAAAGCCGAA.ACggc                                                                                                       >  2:3291219/1‑99 (MQ=21)
  gggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.ccc                                                                                                     >  1:4855030/1‑101 (MQ=37)
   tctgCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cccc                                                                                                    <  2:2830474/98‑1 (MQ=37)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGGACAAGAAACGCCGAA.CCG.cccc                                                                                                    >  2:5186685/1‑101 (MQ=37)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cccc                                                                                                    <  1:1566177/101‑1 (MQ=37)
   ggggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cccc                                                                                                    <  2:1876699/101‑1 (MQ=37)
    gggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCCCAAGAAACGCCGAA.CCG.ccccc                                                                                                   <  2:5930044/101‑2 (MQ=37)
    gggCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.ccccc                                                                                                   <  1:3143140/101‑2 (MQ=37)
    gatcTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.ccccc                                                                                                   <  2:1005038/98‑2 (MQ=37)
      gcgccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctt                                                                                                <  1:2312366/94‑1 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.ccccaga                                                                                                 >  1:2830576/1‑98 (MQ=37)
      gCTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctt                                                                                                <  2:5545759/101‑1 (MQ=37)
       cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.cttt                                                                                               <  1:835836/101‑2 (MQ=37)
       cTCCGCCCCCGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCCAA.CCG.CCCC.cttt                                                                                               >  2:5491341/1‑100 (MQ=37)
        tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctttc                                                                                              <  1:103540/96‑3 (MQ=37)
        tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctttc                                                                                              <  2:6510797/96‑3 (MQ=37)
        tccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctttc                                                                                              <  1:1395103/96‑3 (MQ=37)
        tCCT.CCCCCCCCCCCCCGGCCTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.ctttc                                                                                              <  2:2283069/101‑3 (MQ=255)
         ccgcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCG.cccccagatc                                                                                              >  1:1005073/5‑95 (MQ=38)
           gcccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCAcc                                                                                            >  1:2537826/3‑101 (MQ=255)
            ccccccacacccccGGCTTATTTACCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCAcc                                                                                            <  2:6925011/92‑1 (MQ=37)
            cccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCt                                                                                           <  2:6312669/100‑1 (MQ=38)
             ccccGCA.CCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCtt                                                                                          <  2:2719504/101‑1 (MQ=38)
                 cgcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGc                                                                                       >  1:3950526/5‑101 (MQ=255)
                  gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCt                                                                                      <  2:5404053/98‑1 (MQ=38)
                  gcacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCt                                                                                      <  2:4368177/98‑1 (MQ=38)
                   cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGGCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTc                                                                                     >  2:1308890/3‑101 (MQ=38)
                   cacccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTc                                                                                     <  2:1450717/99‑1 (MQ=39)
                    acccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCa                                                                                    >  2:4295559/2‑101 (MQ=39)
                    acccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCa                                                                                    <  1:536204/100‑1 (MQ=39)
                     cccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaa                                                                                   <  2:6731258/101‑1 (MQ=255)
                      ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaaa                                                                                  <  1:3788312/101‑1 (MQ=255)
                      ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaaa                                                                                  <  1:4025265/101‑1 (MQ=255)
                      ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaaa                                                                                  <  2:2811932/101‑1 (MQ=255)
                      ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaaa                                                                                  <  2:6039369/101‑1 (MQ=255)
                      ccccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATAGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCaaa                                                                                  >  1:4187494/1‑101 (MQ=255)
                       cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAat                                                                                 <  1:6162744/101‑1 (MQ=255)
                       cccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAat                                                                                 >  1:930103/1‑101 (MQ=255)
                        ccGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCCAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAata                                                                                >  2:4554925/1‑101 (MQ=255)
                         cGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGG..CAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATcc                                                                             <  2:5230863/101‑1 (MQ=255)
                          ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATc                                                                              <  1:4889155/101‑1 (MQ=255)
                          ggCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATc                                                                              <  1:6031272/101‑1 (MQ=255)
                           gCGTATTTATCCCAAGAAAACGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATCTcc                                                                             <  1:5603659/101‑1 (MQ=255)
                           gCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATcc                                                                             <  1:1448778/101‑1 (MQ=255)
                           gCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATcc                                                                             >  1:2329222/1‑101 (MQ=255)
                             gTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCg                                                                           <  1:4295781/101‑1 (MQ=255)
                              tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCgg                                                                          <  1:877913/101‑1 (MQ=255)
                              tATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCgg                                                                          <  2:6304778/101‑1 (MQ=255)
                               aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCgcg                                                                         >  1:5196694/1‑99 (MQ=255)
                               aTTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCcgg                                                                         >  2:5313318/1‑101 (MQ=255)
                                tttATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCgggg                                                                        >  1:4244709/1‑101 (MQ=255)
                                   aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCggggggc                                                                     >  1:2203777/1‑99 (MQ=38)
                                   aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCgggggcc                                                                     >  1:7028227/1‑99 (MQ=39)
                                   aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTc                                                                     >  1:813274/1‑101 (MQ=255)
                                   aTCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTc                                                                     <  2:4780087/101‑1 (MQ=255)
                                    tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTcc                                                                    <  2:290266/101‑1 (MQ=255)
                                    tCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGCGGGGcc                                                                    >  2:5413002/1‑101 (MQ=255)
                                     cccAAGGAAAAGCGGCGGCGTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                                   <  2:3274710/101‑1 (MQ=38)
                                     cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                                   >  1:950348/1‑101 (MQ=255)
                                     cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                                   >  2:4800631/1‑101 (MQ=255)
                                     cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                                   >  1:2551045/1‑101 (MQ=255)
                                     cccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCg                                                                   <  2:1529421/101‑1 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                                  <  2:6481762/101‑1 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                                  <  2:5638197/101‑1 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                                  >  1:2681067/1‑101 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCgg                                                                  >  1:4810337/1‑101 (MQ=255)
                                      ccAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCAAAAGAAACGCCGAACCCC.CCCC.TTTCACCTTGGCTCAAATATCCCGGGGGGCCgg                                                                  >  1:3308552/1‑101 (MQ=38)
                                         aGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCggggg                                                               >  1:5735974/1‑101 (MQ=255)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                              <  1:5915894/101‑1 (MQ=255)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                              <  2:3308436/101‑1 (MQ=255)
                                          gAAAAAGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGc                                                              <  1:5010120/101‑1 (MQ=255)
                                            aaaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCgg                                                            >  2:6450543/1‑101 (MQ=255)
                                             aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                           <  2:4091212/101‑1 (MQ=255)
                                             aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                           >  1:580685/1‑101 (MQ=255)
                                             aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                           >  2:5356732/1‑101 (MQ=255)
                                             aaaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGa                                                           <  2:2575137/101‑1 (MQ=255)
                                              aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAg                                                          <  1:4249842/101‑1 (MQ=255)
                                              aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAg                                                          <  2:3261211/101‑1 (MQ=255)
                                              aaGCAGCGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAg                                                          >  2:387688/1‑101 (MQ=255)
                                               agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGTTGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..c                                                       <  1:4555164/101‑1 (MQ=37)
                                               agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..c                                                       <  1:4800874/101‑1 (MQ=255)
                                               agcagcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..c                                                       >  1:6706418/1‑101 (MQ=255)
                                                   gcGGCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTTACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..ccccc                                                   <  2:3388154/101‑3 (MQ=37)
                                                    cGGCTTGCCCCGACGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccg                                                  <  2:6559502/101‑4 (MQ=25)
                                                     ggCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccgg                                                 >  2:5097781/1‑97 (MQ=25)
                                                     ggCTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccgg                                                 >  2:1095687/1‑97 (MQ=25)
                                                       cTTGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGGGGGG..Ccc..gcggcg                                             >  2:399100/1‑97 (MQ=25)
                                                        ttGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgg                                              >  2:1013453/1‑94 (MQ=25)
                                                         tGCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcggg                                             >  2:5563753/1‑93 (MQ=21)
                                                          tccccGACGGGCTTG.TTTTGATGGCAGAAGAAACGTCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATAGCCCCGGCGGCCAGGGGCGGCG..cccctggcggag                                            >  2:4624593/2‑92 (MQ=255)
                                                          gCCCCGAAGGCCTTG.TTTTGATGGCACAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgggg                                            >  1:190550/1‑92 (MQ=21)
                                                           ccccGAAGGCCTTG.TTTTGATGGCCCAAGAAACGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..cccccggcgggtc                                           <  2:7027853/101‑11 (MQ=18)
                                                                                         aaaCGCCGAA.CCGCCCCC.TTTCCCCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcc            <  2:4176276/101‑1 (MQ=2)
                                                                                           ccGCCGAC.CCGCCCCC.TTTCCCCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGCC..CCCCGGGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc          <  1:6478459/100‑1 (MQ=2)
                                                                                           aCGCCGAC.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCGGGGGGCGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACAcccc          <  1:6207504/101‑1 (MQ=2)
                                                                                            cGCCGAA.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCa         <  2:5735672/101‑1 (MQ=0)
                                                                                              ccGAC.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGTTCCGGGGGCGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAgc       <  2:2203682/101‑1 (MQ=0)
                                                                                              ccGAC.CCGCCCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAgc       <  1:1013488/101‑1 (MQ=2)
                                                                                                    ccG.CCCCTTTTCGCCTTGGCTCAAATATCCCGGGGGTCGGGGGGGGGAGCCCCC..GGCGGGTCGCCTCGGGCA..AAGCCCGAGGCGAAACACCCCAGCGCCac  <  2:5196429/101‑1 (MQ=2)
                                                                                                       |                                                                                                      
GCGGGGGCTCCG.CCCCGCACCCCCCGGCGTATTTATCCCAAGAAAAAGCAGCGGCTTGCCCCGAAGGCCTTGTTTTTGATGGCACAAGAAACGCCGAA.CCG.CCCC.TTTCACCTTGGCTCAAATATCCCGGGGGTCCGGGGGCGGAG..CCC..GCCGGGTCTCCTCCGGGAGCAAGGCCGAGGCGAGACACCCCAGCGCCAC  >  NC_014034/2185830‑2186027

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG