Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 904,926 Δ1 bp 96.6% intergenic (+30/‑166) RCAP_rcc00821 → / → murB hypothetical protein/UDP‑N‑acetylmuramate dehydrogenase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_014034904,9260T→.96.6% 4.0 0.0 ‑0.4116intergenic (+30/‑166)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (64/48):  ref base (1/1):  total (66/50)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.61e-01

Only 100 of 144 total aligned reads displayed.

TGGCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..GC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGC  >  NC_014034/904826‑905020
                                                                                                    |                                                                                                    
tGGCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCctct                                                                                                      <  2:261440/101‑1 (MQ=255)
 gccAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCtt                                                                                                     <  1:1086432/99‑1 (MQ=255)
  gCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCCCCttt                                                                                                    >  1:2023479/1‑101 (MQ=255)
  gCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttt                                                                                                    <  2:1724527/101‑1 (MQ=255)
  gCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttt                                                                                                    <  2:3499531/101‑1 (MQ=255)
   cgggggggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCtttt                                                                                                   <  1:3333330/96‑1 (MQ=255)
   cAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCtttt                                                                                                   <  1:2286923/101‑1 (MQ=255)
   cAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCtttt                                                                                                   <  2:1125433/101‑1 (MQ=255)
    gggggggCTGAGCGTGCTGCGCTGGCCGGGGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttttg                                                                                                  <  1:3269240/97‑2 (MQ=255)
    aGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttttg                                                                                                  <  1:1637061/101‑2 (MQ=255)
    aGGCGGGCGGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttttg                                                                                                  <  2:1376454/101‑2 (MQ=255)
     ggCGGGCTGGGCGTGCTGCGCGGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCttttgg                                                                                                 <  1:996369/101‑3 (MQ=255)
       ggggTTGATCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGCCGCTGCCCCCCTCttttggct                                                                                               <  2:3366491/100‑5 (MQ=38)
        gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.CTTTG.GCT..g                                                                                          >  1:1725835/1‑101 (MQ=255)
          gCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTTTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCGGCCGCCCTCTTTTGg.gtgcg                                                                                           >  1:2278930/1‑97 (MQ=255)
           cTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCGGGCCCCCTCTTTTGg.gtgcgc                                                                                          >  2:1738570/1‑96 (MQ=255)
           cTGAGCGTGCTGCGCTGGACGGTGCGGCGGCTCTTTGCCGGGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCCTCTTTTTgctgcgc                                                                                           >  2:3534609/1‑95 (MQ=38)
            tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGGTGCCCCCCTC.TTTTG.GCT..Gcgcc                                                                                      >  1:3692708/1‑98 (MQ=255)
            tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgcc                                                                                      >  2:3799087/1‑98 (MQ=255)
            tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCCGCCCCCCTC.TTTTG.GCT..Gcgcc                                                                                      <  2:4405218/101‑4 (MQ=255)
              aGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTgctgcgccga                                                                                        >  2:2158273/1‑92 (MQ=255)
               gCGTGCTGCGCTGGCCGGTGCGGGCGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaa                                                                                   <  1:4391059/101‑7 (MQ=25)
               gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaa                                                                                   <  2:4298826/101‑7 (MQ=255)
                cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaaa                                                                                  <  1:543397/101‑8 (MQ=37)
                 gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaaat                                                                                 <  1:2864022/101‑9 (MQ=37)
                    ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCa                                                                              >  1:2996326/1‑101 (MQ=25)
                     ggcgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAc                                                                             <  1:1592884/100‑1 (MQ=25)
                         cTGGCCGGTGCGGTGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAccccc                                                                         >  1:3649430/1‑101 (MQ=21)
                         cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAccccc                                                                         >  2:116562/1‑101 (MQ=21)
                          tGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAcccccc                                                                        >  2:2118434/1‑101 (MQ=21)
                            gCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.gg                                                                     >  1:1830108/1‑101 (MQ=21)
                             ccGGTGCGGCGGCTCTTTTCCGCGATGCTCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTGG.CTTC.GC..CAAAATTCACCCCCC.GGc                                                                    <  2:422682/101‑1 (MQ=18)
                               ggTGCGGCGGCGCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGt                                                                  <  2:4606809/101‑1 (MQ=21)
                                    ggcggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTg                                                             <  1:1331149/101‑1 (MQ=21)
                                      ggggtctTTGCCGCGATGCGCCGGCGCGGCGGAGGGCGGGGCAGGGGGGCGCGGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCa                                                           <  1:4541392/97‑1 (MQ=17)
                                        gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAcc                                                         >  1:55546/1‑101 (MQ=21)
                                        gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAcc                                                         >  2:1479203/1‑101 (MQ=21)
                                           ctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAgg                                                      <  1:2501909/101‑1 (MQ=21)
                                             ttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaa                                                    <  2:2882401/101‑1 (MQ=21)
                                              tGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaa                                                   >  2:4529002/1‑101 (MQ=21)
                                                   cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGcc                                              >  2:2801335/1‑101 (MQ=21)
                                                   cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGcc                                              >  2:3971655/1‑101 (MQ=21)
                                                    gATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGccc                                             <  1:3283003/101‑1 (MQ=21)
                                                     aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCg                                            >  1:2693186/1‑101 (MQ=21)
                                                      tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                           <  2:140051/101‑1 (MQ=21)
                                                      tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                           >  1:3447674/1‑101 (MQ=21)
                                                      tGCACCGGCGCGACTGAGGGCGGGGCAGAGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                           >  1:2438607/1‑101 (MQ=21)
                                                         cccGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGt                                        <  2:374482/100‑1 (MQ=21)
                                                              cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CGAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                   <  1:2113328/101‑1 (MQ=21)
                                                               gcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTcc                                  >  2:1756770/1‑101 (MQ=21)
                                                                cgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAATATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCt                                 <  1:3948389/101‑1 (MQ=21)
                                                                cgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.CGCGCATTTGCACCACGAAAAAGCCCGGCGTTTTTCCt                                 >  1:7330/1‑101 (MQ=18)
                                                                cgACTGAGGGCGGCGCAGGAGGGCGCTGCCTCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCCTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCt                                 <  1:440122/101‑1 (MQ=18)
                                                                 gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTg                                >  2:3213770/1‑101 (MQ=21)
                                                                 gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCACGGCGTTCTTCATg                                >  2:433502/1‑101 (MQ=18)
                                                                 gACTGAGGGCGGGGCA.GGGGGCGCGGCCCCCCCCCTTTGG.CCTC.GC..CAAAAATTCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTgg                               <  2:2221254/101‑1 (MQ=25)
                                                                  aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGGAAAACGCCGGCGTTTTTCCTgg                               >  2:740969/1‑101 (MQ=18)
                                                                    tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGGGC..CAAAATTCACCCCCC.GGCGTATTTGCACCACGAAAAAGCCCGGCGTTTTTCCTGGTc                             >  2:1522456/1‑101 (MQ=21)
                                                                    tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTc                             >  1:732878/1‑101 (MQ=21)
                                                                     gAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCg                            >  1:1305927/1‑101 (MQ=21)
                                                                       gggCGGGGCGGGGGGACGCTCCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaa                          <  2:1446298/101‑1 (MQ=21)
                                                                       gggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaa                          >  2:1875790/1‑101 (MQ=21)
                                                                       gggCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..GCGGCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGaa                          <  2:1757080/101‑1 (MQ=21)
                                                                              gCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGcc                   <  2:2071951/101‑1 (MQ=21)
                                                                               cAGGGGGGCGCTTCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCg                  <  2:359727/101‑1 (MQ=21)
                                                                                aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                 >  2:2368886/1‑101 (MQ=25)
                                                                                aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                 >  1:2557930/1‑101 (MQ=25)
                                                                                 ggggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggg                >  2:4305059/1‑101 (MQ=25)
                                                                                  gggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggg               <  1:3349761/101‑1 (MQ=25)
                                                                                  gggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggg               >  1:3282318/1‑101 (MQ=25)
                                                                                   ggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACCCCggtcg              >  1:2552987/1‑98 (MQ=21)
                                                                                    gggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggggg             <  1:2266876/101‑1 (MQ=25)
                                                                                     ggCGCTGCCCCCCTCTTTTGGCTCC..GC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa            <  1:3534867/101‑1 (MQ=25)
                                                                                     ggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa            <  1:3868298/101‑1 (MQ=25)
                                                                                     ggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa            <  1:3837191/101‑1 (MQ=25)
                                                                                     ggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGa            >  1:3763516/1‑101 (MQ=25)
                                                                                      gcgcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAg           <  1:1522565/101‑1 (MQ=25)
                                                                                      gcccTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTGCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAg           <  1:4032029/98‑1 (MQ=21)
                                                                                       cgcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt          >  2:110867/1‑101 (MQ=37)
                                                                                       cgcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt          >  1:1726610/1‑101 (MQ=37)
                                                                                       cgcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt          >  1:3426220/1‑101 (MQ=37)
                                                                                        gcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTc         >  1:140895/1‑101 (MQ=37)
                                                                                         cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc        >  1:1080083/1‑101 (MQ=37)
                                                                                         cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc        >  1:2123970/1‑101 (MQ=37)
                                                                                         cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc        >  2:2596578/1‑101 (MQ=37)
                                                                                         cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc        <  2:1726491/101‑1 (MQ=37)
                                                                                         cTGCCCCCCT.ATTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCAACGGGAAAAAGAACGGCGTTNNNNNTGGTCGTAATACACCGGGGGGAggtc        >  1:2393677/1‑98 (MQ=12)
                                                                                          tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg       >  2:3719575/1‑101 (MQ=37)
                                                                                          tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg       <  2:2692990/101‑1 (MQ=37)
                                                                                          tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg       >  1:39017/1‑101 (MQ=37)
                                                                                           gCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc      >  2:1924402/1‑101 (MQ=37)
                                                                                           gCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGCCAAAAAACGCCGGGGGGGGTCCGc      >  1:4033166/1‑101 (MQ=25)
                                                                                            ccccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt     <  2:807826/101‑1 (MQ=37)
                                                                                            ccccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt     >  2:1549671/1‑101 (MQ=37)
                                                                                             cccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt    <  1:694534/101‑1 (MQ=37)
                                                                                             cccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt    <  1:4526498/101‑1 (MQ=37)
                                                                                             cccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt    <  1:3799372/101‑1 (MQ=37)
                                                                                              ccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTg   >  1:3595983/1‑101 (MQ=37)
                                                                                              ccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGACCGCTTg   <  2:444012/101‑1 (MQ=37)
                                                                                                 ccttttttG.GCC..GCCGCAAAATTCCCCCCCCGGGCTTTTTTCCACAAGAAAAAACCCGGGCGTTTTCCCTGGCCAAAATCCGCCGGGGGGAGTCCGCTTGc  <  2:1604907/98‑1 (MQ=255)
                                                                                                    |                                                                                                    
TGGCAGGCGGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..GC..CAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGC  >  NC_014034/904826‑905020

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