Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 904,934 +GC 100% intergenic (+38/‑158) RCAP_rcc00821 → / → murB hypothetical protein/UDP‑N‑acetylmuramate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034904,9341.→G100.0% 60.0 157intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (72/65):  ref base (4/16):  total (76/81)
Fisher's exact test strand distribution p-value = 7.83e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.84e-01
*NC_014034904,9342.→C100.0% 60.0 157intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (72/65):  ref base (4/16):  total (76/81)
Fisher's exact test strand distribution p-value = 7.83e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.21e-01

Only 100 of 169 total aligned reads displayed.

GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTGGCT..G..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025
                                                                                                     ||                                                                                                
gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..g                                                                                                 <  1:522312/101‑1 (MQ=255)
gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..g                                                                                                 <  1:2094628/101‑1 (MQ=255)
 ggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..c                                                                                              >  1:1134854/1‑101 (MQ=255)
 ggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..c                                                                                              <  2:5272508/101‑1 (MQ=255)
   cTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgc                                                                                            >  2:5375651/1‑99 (MQ=255)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgcca                                                                                          <  1:955043/101‑5 (MQ=255)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgcca                                                                                          <  1:2601811/101‑5 (MQ=255)
       gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgccaaa                                                                                        <  1:4298836/101‑7 (MQ=255)
       gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgccaaa                                                                                        <  2:4681084/101‑7 (MQ=255)
       gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgccaaa                                                                                        <  1:1143380/101‑7 (MQ=255)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgccaaaa                                                                                       >  2:3774536/1‑94 (MQ=37)
          tgctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..cgccaaaatt                                                                                     <  1:2449985/101‑10 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTc                                                                                   >  2:4909067/1‑101 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTc                                                                                   <  2:2426117/101‑1 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCa                                                                                  <  1:1767440/101‑1 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAc                                                                                 <  2:3292211/101‑1 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAc                                                                                 >  1:5432075/1‑101 (MQ=25)
              gcgcTGGCCGGTGCGGCGGCTTTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAcc                                                                                >  1:3287348/1‑101 (MQ=21)
              gcgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAcc                                                                                >  1:5517984/1‑101 (MQ=25)
              gcgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAcc                                                                                <  2:453255/101‑1 (MQ=25)
               cggtGGCCGGTGCGGCGCCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAccc                                                                               <  2:1869315/98‑1 (MQ=21)
                gcgggCCGGGGGGGGGCCTTTTGGCGGGGATCCGCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAcccc                                                                              <  1:1167552/98‑1 (MQ=2)
                 cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAccccc                                                                             <  2:583419/101‑1 (MQ=25)
                       ggTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCCGCGt                                                                       >  1:2335112/1‑101 (MQ=21)
                         tGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTAt                                                                     <  2:5687281/101‑1 (MQ=21)
                         tGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCAACCCCCGGCGGAt                                                                     >  1:5534224/1‑101 (MQ=21)
                         ggcggcgGCTTTTGGCCGCGATGCACGGGCGCGACTGAGGGCGGGGCAGGGGGACGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTAt                                                                     <  2:5010099/100‑1 (MQ=18)
                            ggcgtttctTTGCCGCGAGGGACAGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCTCCCCCCGGCGTATTTg                                                                  <  2:2311579/95‑1 (MQ=17)
                            ggcggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTg                                                                  >  1:4046096/1‑101 (MQ=21)
                             gcggcTCTTTGCCGCGATGCGCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGc                                                                 <  1:480374/101‑1 (MQ=18)
                              cggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCa                                                                >  2:3549860/1‑101 (MQ=21)
                                gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCT..G..CGGC.AAAATTCACCCCCCGGCGTATTTGCAcc                                                              <  1:3832861/101‑1 (MQ=21)
                                 ctctTTGCCGCGATGCATCGGCGCGACTGAGGCCGGGGCAGGGGGGCGCGGCCCCCCTC.TTTTGGCTGCG..C..C.AAACGTAACCCCCCGGCGCATTTGCACCa                                                             >  2:3676418/1‑101 (MQ=18)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAg                                                            <  2:666840/101‑1 (MQ=21)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAg                                                            >  2:2377827/1‑101 (MQ=21)
                                   ctTTGCCGCGGTGCACTGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAgg                                                           <  1:147078/101‑1 (MQ=18)
                                    tttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGa                                                          >  2:2190187/1‑101 (MQ=21)
                                    tttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGa                                                          >  2:4091781/1‑101 (MQ=21)
                                      tGCCGCGATGCACCGCCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGaaa                                                        <  2:3010441/101‑1 (MQ=18)
                                       gcggcgATGCCCCGGGGCGACTGAGGTCGGGGCGGGGGGGCGCGGCCCCCCTC.TTTTGGCT..GC.C..C.AAAAATTCCCCCCCGGCGTATTTGCACCAGGaaaaa                                                      <  2:829823/98‑1 (MQ=16)
                                       gCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGaaaa                                                       >  2:5665693/1‑101 (MQ=21)
                                       gCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCCGCGTATTTGCACCAGGaaaa                                                       >  1:1488293/1‑101 (MQ=21)
                                        ccGCGGTGCACCTGCGCGACTGAGGGCGGGGCAGGTGGGCGCTGCCCCCCTC.TTTTGGCT..GCTC..C.AAAATTCTCCCCCCGGCGTATTTGCACCAGGaaaaa                                                      <  1:908970/101‑1 (MQ=18)
                                        ccGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGaaaaa                                                      <  2:3246143/101‑1 (MQ=21)
                                           cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTCTTTGCACCAGGAAAAAGcc                                                   >  2:2725095/1‑101 (MQ=21)
                                            gATGGACCGGCGCGACTGAGTGCGGGGCAGGGGCGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCCCCCCCCGGCGTATTTGCACCAGGAAAAAGccc                                                  <  1:1548111/101‑1 (MQ=18)
                                             aTGTACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                 <  1:799938/101‑1 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                 >  2:3303502/1‑101 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                 >  2:1529971/1‑101 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGCGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                 <  2:3961107/101‑1 (MQ=21)
                                              tggaCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCgg                                                <  1:2196170/98‑1 (MQ=18)
                                               gCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGc                                               >  1:492972/1‑101 (MQ=21)
                                                   cGGCGCGACTGGGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                           <  1:2231628/101‑1 (MQ=18)
                                                   cGGCGCGACTGAGGGCGGGGGAGGGGGGCGCTTCCCCCCTC.TTTTGGCT..GCCC..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                           <  1:5614486/101‑1 (MQ=21)
                                                   cGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                           >  2:4752275/1‑101 (MQ=21)
                                                     gcgcGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttttt                                         >  2:963612/1‑101 (MQ=21)
                                                       gcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTcc                                       <  1:3034500/101‑1 (MQ=21)
                                                          aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTgg                                    >  2:2343992/1‑101 (MQ=21)
                                                          aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTgg                                    >  1:4104220/1‑101 (MQ=21)
                                                          aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGGTTTTCCTgg                                    >  1:1834249/1‑101 (MQ=18)
                                                           cTGAGGGCGGGGCAGGGTGGCCCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                   <  1:4908116/101‑1 (MQ=21)
                                                           cTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                   >  2:1091195/1‑101 (MQ=21)
                                                           cTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                   <  2:4055050/101‑1 (MQ=21)
                                                                ggCGGGGCAGGGGGGCGCTGCTCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaaa                              >  1:3284376/1‑101 (MQ=21)
                                                                ggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaaa                              >  2:2222240/1‑101 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCTCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATAc                           <  2:2032584/101‑1 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATAc                           <  2:1966245/101‑1 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGGTTTTCCTGGTCGAAATAc                           >  1:4077299/1‑101 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCCGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATAc                           >  1:1572025/1‑101 (MQ=21)
                                                                    gggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACg                          <  2:3089487/101‑1 (MQ=21)
                                                                     ggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGc                         >  2:1429558/1‑101 (MQ=21)
                                                                     ggCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGc                         <  1:3148845/101‑1 (MQ=21)
                                                                      gCAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGcc                        >  2:736497/1‑101 (MQ=21)
                                                                       cAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCg                       <  1:5271844/101‑1 (MQ=21)
                                                                       cAGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGCATTTGCACCAGG.AAAAGCCCGGCGTTTTTCTTGGTCGAAATACGCCgg                      >  2:1452935/1‑101 (MQ=18)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      >  2:4331018/1‑101 (MQ=25)
                                                                         ggggggCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggg                     <  2:4358810/101‑1 (MQ=25)
                                                                          gggggCGCTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggg                    <  2:1261250/101‑1 (MQ=25)
                                                                              gcgcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAg                <  1:3655754/101‑1 (MQ=25)
                                                                               cgcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAgg               >  2:2635387/1‑100 (MQ=25)
                                                                               cgcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAgc               >  1:3674772/1‑100 (MQ=25)
                                                                               cgcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt               <  1:5714489/101‑1 (MQ=25)
                                                                                gcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCCCCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc              <  1:1530064/101‑1 (MQ=37)
                                                                                gcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc              >  2:4327751/1‑101 (MQ=37)
                                                                                gcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc              <  1:2981585/101‑1 (MQ=37)
                                                                                gcTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc              >  2:1651366/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCCCCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             <  1:444896/101‑1 (MQ=25)
                                                                                  tGCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg            <  1:2377961/101‑1 (MQ=37)
                                                                                   gCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc           >  1:5661407/1‑101 (MQ=37)
                                                                                   gCCCCCCTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc           >  2:1397248/1‑101 (MQ=37)
                                                                                    ccccccTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          <  1:3774779/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         >  2:134521/1‑101 (MQ=38)
                                                                                     cccccTC.TTTTGGCTGCG..C..C.AAAATTCACCCAGCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  1:4503802/101‑1 (MQ=25)
                                                                                     cccccT.TTTTTGGCT..G..CG.CAAAAATCCCCCCCCCGGCGTATTTGCACAGGAAAAAAGCCCGGCGTTTTTCGTGGCCAAAATCCGCCGGGGGGAGTCCGCtt         <  2:2904153/101‑1 (MQ=255)
                                                                                       cccTC.TTTTGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGc       <  2:1794067/101‑1 (MQ=37)
                                                                                          ctcttttGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg     >  1:5234252/4‑101 (MQ=37)
                                                                                          ctcttttGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg     <  1:2725244/98‑1 (MQ=38)
                                                                                           tcttttGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    >  1:1279760/3‑101 (MQ=37)
                                                                                            cttttGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGAGGAGTCCGCTTGCGGAc   >  2:5195239/2‑101 (MQ=38)
                                                                                             ttttGGCTGCG..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg  >  1:5724090/1‑101 (MQ=255)
                                                                                                     ||                                                                                                
GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTGGCT..G..C..C.AAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 31 ≤  ATCG  < 39 ≤  ATCG