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% 41.3 110intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (55/40):  ref base (1/13):  total (57/53)
Fisher's exact test strand distribution p-value = 3.45e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.33e-01
*NC_014034904,9342.→C100.0% 45.0 110intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (55/40):  ref base (2/13):  total (57/53)
Fisher's exact test strand distribution p-value = 1.63e-03
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 5.70e-01

Only 100 of 118 total aligned reads displayed.

GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..G.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025
                                                                                                      ||                                                                                             
gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..g                                                                                              <  2:3240911/101‑1 (MQ=255)
   cTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgc                                                                                          >  1:1446889/1‑99 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGGCCCCCTC.TTTTG.GCT..G.cgcc                                                                                         >  2:2621868/1‑98 (MQ=38)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcc                                                                                         <  2:1582500/101‑4 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcc                                                                                         <  2:4347192/101‑4 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcc                                                                                         >  2:319084/1‑98 (MQ=255)
    ggggCGTGCTGCGCGGGCCGGTGCGGCGCCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcc                                                                                         <  1:13756/98‑4 (MQ=37)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcca                                                                                        <  2:693550/101‑5 (MQ=255)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcca                                                                                        <  2:2440821/101‑5 (MQ=255)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgcca                                                                                        <  1:1926765/101‑5 (MQ=255)
      aGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaa                                                                                       <  2:2112757/101‑6 (MQ=255)
      aGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaa                                                                                       <  2:4088139/101‑6 (MQ=255)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaaaa                                                                                     <  2:906949/101‑8 (MQ=37)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaaaat                                                                                    <  2:3399475/101‑9 (MQ=37)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaaaat                                                                                    <  2:171278/101‑9 (MQ=37)
         gTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..G.cgccaaaat                                                                                    <  1:395957/101‑9 (MQ=37)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTc                                                                                    <  2:1672097/101‑1 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTc                                                                                    <  1:1895865/101‑1 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTc                                                                                    >  1:557274/1‑101 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCa                                                                                   <  1:4153060/101‑1 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCa                                                                                   >  2:3388094/1‑101 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAc                                                                                  <  2:2944676/101‑1 (MQ=25)
             tgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAc                                                                                  <  1:3246206/101‑1 (MQ=25)
               cgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCCAAATTCAccc                                                                                >  1:1195311/1‑101 (MQ=25)
               cgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAccc                                                                                <  2:1607839/101‑1 (MQ=25)
               cgcCGGACGGTGCGGTGGCTCTTTGCCGCCATGCCCCGGAGCGACTGAGGGCGGGGCACGGGTGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAccc                                                                                <  1:3687503/101‑1 (MQ=12)
                gcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAcccc                                                                               <  1:378125/101‑1 (MQ=25)
                  tGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCAcccccc                                                                             <  1:24874/101‑1 (MQ=21)
                  tGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCC.CCTTTTG.GCTGCG.CCAAAAATCAcccccc                                                                             >  1:167696/1‑101 (MQ=21)
                   ggCCGGTGCGGCGGCTCTTTGCCGCCATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.g                                                                           <  2:2168952/101‑1 (MQ=21)
                     ccGGAGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.CGc                                                                         >  2:157225/1‑101 (MQ=18)
                        gTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTa                                                                      >  1:3159432/1‑101 (MQ=21)
                         ggcggcgCCTCTTTGCCGCGATGCGCCGGCGCGAGGGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTAt                                                                     <  1:895076/100‑1 (MQ=17)
                          gcggcgGCTCTTTGCCGCTATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTGCTGC..G.CCCAAATTCACCCCCC.CGCGGAtt                                                                    >  1:826138/1‑101 (MQ=18)
                           cggcggCTCTTTGCCGCGATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTAttt                                                                   <  2:2649842/101‑1 (MQ=18)
                           cggcggCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTAttt                                                                   >  1:3639948/1‑101 (MQ=21)
                           cggcggCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTAttt                                                                   >  1:4179475/1‑101 (MQ=21)
                            ggcggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTg                                                                  <  2:5477/101‑1 (MQ=21)
                              cggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTGG.CTTC.G.CCAAAATTCAACCCCC.GGCGTATTTGCa                                                                >  1:2804400/1‑101 (MQ=18)
                              cggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCTCCCCCC.GGCGTATTTGCa                                                                <  1:3062422/101‑1 (MQ=21)
                              cggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCa                                                                <  2:304854/101‑1 (MQ=21)
                              cggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCa                                                                >  2:1488810/1‑101 (MQ=21)
                               ggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCAc                                                               >  1:2222468/1‑101 (MQ=21)
                                 ctctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGGAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTCTTTGCACCa                                                             >  1:3458335/1‑101 (MQ=21)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAg                                                            >  1:502551/1‑101 (MQ=21)
                                     ttGCCGCGATGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGaa                                                         <  1:2696799/101‑1 (MQ=21)
                                     ttGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCCCCCCCC.GGCGTATTTGCACCAGGaa                                                         <  2:283737/101‑1 (MQ=21)
                                      tGCCGCGGTGCCCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaa                                                        <  1:1584287/101‑1 (MQ=18)
                                      tGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaa                                                        >  1:1123130/1‑101 (MQ=21)
                                       gCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGaaaa                                                       <  2:3944538/101‑1 (MQ=21)
                                        ccGCTATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCCCCCCCC.GGCGTATTTGCCCCAGGaaaaa                                                      <  2:2649621/101‑1 (MQ=18)
                                            gATGCACCGGCGCGACTGAGGGCGGGGCAAGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGccc                                                  <  2:1886978/101‑1 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCg                                                 >  2:1209026/1‑101 (MQ=21)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCg                                                 >  1:1442764/1‑101 (MQ=21)
                                              tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                                <  1:1708218/101‑1 (MQ=21)
                                              tGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCgg                                                >  2:3038996/1‑101 (MQ=21)
                                               gCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGc                                               >  1:4065951/1‑101 (MQ=21)
                                                     gcgcGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGttttt                                         >  1:1685855/1‑101 (MQ=21)
                                                      cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                        >  1:2598423/1‑101 (MQ=21)
                                                          aCTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTgg                                    >  1:824331/1‑101 (MQ=21)
                                                           cggAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCGAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGt                                   <  2:1368954/99‑1 (MQ=18)
                                                            tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTc                                  >  2:3543946/1‑101 (MQ=21)
                                                                 gCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAAt                             <  2:2890473/101‑1 (MQ=21)
                                                                 gCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAAt                             <  2:2663521/101‑1 (MQ=21)
                                                                  cGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATa                            <  2:3099902/101‑1 (MQ=21)
                                                                  cGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATa                            >  1:2280249/1‑101 (MQ=21)
                                                                    gggCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACg                          <  2:4157842/101‑1 (MQ=21)
                                                                      gCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGcc                        >  2:2208460/1‑101 (MQ=21)
                                                                      gCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGcc                        >  2:3973903/1‑101 (MQ=21)
                                                                       cAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCg                       <  1:3280828/101‑1 (MQ=21)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      <  2:3872172/101‑1 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      >  1:2936171/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      >  1:3342301/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgg                      <  2:3413829/101‑1 (MQ=25)
                                                                         ggggggCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggg                     >  1:3637746/1‑101 (MQ=25)
                                                                          gggggCGCTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCgggg                    >  2:2250113/1‑101 (MQ=25)
                                                                              gcgcTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAg                >  2:2710536/1‑101 (MQ=25)
                                                                               cgcTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt               >  2:2549927/1‑101 (MQ=37)
                                                                               cgcTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGt               >  1:1835143/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             >  1:3358973/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             <  2:557222/101‑1 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTcc             >  2:1331138/1‑101 (MQ=37)
                                                                                  tcccccccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCggcgggtgaccc            >  2:3631711/3‑95 (MQ=25)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg            >  1:2288450/1‑101 (MQ=37)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCg            >  2:1827058/1‑101 (MQ=37)
                                                                                  ccccccccTTTTTTTG.GCT..GCCCAAAAATCCCCCCCCGGGTGTTTTTGCACCAGAAAAAACCCCGGTTTTTTCCCTGGTCAAAATCCGCCGGGGGGAGTcc             <  2:794650/99‑1 (MQ=255)
                                                                                   gCCCCCCTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGc           >  2:1374464/1‑101 (MQ=37)
                                                                                   gCCCCCCTC.TTGTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACGAGGAAAAAGCCCGGCGTTTTTCGTGGTCGAAATACGCCGGGGGGAGTCCGc           <  2:1040422/101‑1 (MQ=25)
                                                                                    ccccccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          <  1:1786442/101‑1 (MQ=37)
                                                                                    ccccccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          >  2:2716944/1‑101 (MQ=37)
                                                                                    ccccccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCt          >  1:1535675/1‑101 (MQ=37)
                                                                                     cccccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCtt         <  2:199429/101‑1 (MQ=37)
                                                                                       cccTC.TTTTG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGc       >  2:2132817/1‑101 (MQ=37)
                                                                                          ctcttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgg     >  1:2360253/4‑101 (MQ=37)
                                                                                          ctcttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCgc     >  1:361501/4‑100 (MQ=37)
                                                                                           tcttttG.GCTGCG.CCAAAATTCCCCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGa    <  1:2442299/99‑1 (MQ=37)
                                                                                            cttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   <  2:2761896/100‑1 (MQ=255)
                                                                                            cttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc   >  1:1538238/2‑101 (MQ=255)
                                                                                             ttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg  >  1:1821753/1‑101 (MQ=255)
                                                                                             ttttG.GCTGCG.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg  >  1:1190071/1‑101 (MQ=255)
                                                                                                      ||                                                                                             
GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..G.CCAAAATTCACCCCCC.GGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGGTCGAAATACGCCGGGGGGAGTCCGCTTGCGGACG  >  NC_014034/904834‑905025

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 24 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG