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% 39.6 164intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (84/66):  ref base (1/13):  total (85/79)
Fisher's exact test strand distribution p-value = 4.16e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.00e-01
*NC_014034904,9342.→C100.0% 39.4 164intergenic (+38/‑158)RCAP_rcc00821/murBhypothetical protein/UDP‑N‑acetylmuramate dehydrogenase
Reads supporting (aligned to +/- strand):  new base (84/66):  ref base (1/13):  total (85/79)
Fisher's exact test strand distribution p-value = 4.16e-04
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.71e-01

Only 100 of 177 total aligned reads displayed.

GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..GC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGACGGGGGGC  >  NC_014034/904834‑905031
                                                                                                      ||                                                                                                    
gggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCC.CCTTTTG.GCT..g                                                                                                     >  1:1681769/1‑101 (MQ=255)
 ggCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gc                                                                                                    >  1:2429349/1‑101 (MQ=255)
    tGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgcc                                                                                                 <  1:1186027/101‑4 (MQ=255)
     gAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgcca                                                                                                <  2:3171964/101‑5 (MQ=255)
      aGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaa                                                                                               <  1:1588462/101‑6 (MQ=255)
       gCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaa                                                                                              <  2:363125/101‑7 (MQ=255)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaaa                                                                                             <  2:282475/101‑8 (MQ=37)
        cGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCT..Gcgccaaaa                                                                                             <  2:3242659/101‑8 (MQ=37)
          tgctgcGCTGGACGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTTGGGGCGGGGGAGGGGGGCGCTGCCCCCCTC.GTTTG.GCT..Gcgctaaaatt                                                                                           >  1:3904403/1‑92 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTc                                                                                          <  1:3853261/101‑1 (MQ=25)
           gctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGAGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTc                                                                                          <  1:4606020/101‑1 (MQ=25)
            ctgcGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCa                                                                                         <  1:2890388/101‑1 (MQ=25)
               cgcTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGGCTGCCCCCCTC.TTTTG.GCTGCGC..CACAATTCAccc                                                                                      >  1:2200259/1‑101 (MQ=25)
                 cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGGACTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTTAccccc                                                                                    >  1:4007642/1‑101 (MQ=21)
                 cTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAccccc                                                                                    <  1:1151868/101‑1 (MQ=25)
                  tGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAcccccc                                                                                   <  2:2433260/101‑1 (MQ=21)
                  tGGCCGGTGCGGCGGCTCTTTGCCGCGATGAACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CCAAAATCAcccccc                                                                                   >  1:4226205/1‑101 (MQ=21)
                  tGGCCGGGGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCAcccccc                                                                                   <  1:2250421/101‑1 (MQ=21)
                      cGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCg                                                                               <  2:216038/101‑1 (MQ=21)
                       ggTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGt                                                                              >  2:2758757/1‑101 (MQ=21)
                       ggTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGt                                                                              >  2:356191/1‑101 (MQ=21)
                       ggTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGt                                                                              >  1:2843892/1‑101 (MQ=21)
                       ggTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCC.CCTTTTG.GCTGCGC..CAAAATTCCCCCCCCGGCGt                                                                              >  2:4046564/1‑101 (MQ=21)
                         tGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CCAAATTTACCCCCCGGCGCAt                                                                            >  2:2244705/1‑101 (MQ=21)
                         tGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGGAt                                                                            >  2:1288684/1‑101 (MQ=21)
                         tGCGGCGGCGCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTAt                                                                            <  2:286612/101‑1 (MQ=21)
                          gcggcgGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGATGGCCCCCTCTTTTTGCTGC..GC..CAAAATTCCCCCCCCGGCGTGtt                                                                           >  1:4478073/1‑101 (MQ=255)
                           cggcggCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTAttt                                                                          <  1:214760/101‑1 (MQ=21)
                           cggccgCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCCGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTAttt                                                                          <  1:1269193/101‑1 (MQ=18)
                            ggcggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTg                                                                         <  2:226556/101‑1 (MQ=21)
                            ggcggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCAC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTg                                                                         >  2:546623/1‑101 (MQ=21)
                             gccgcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGc                                                                        <  2:2984096/98‑1 (MQ=18)
                               ggcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCAc                                                                      >  1:3087841/1‑101 (MQ=21)
                               ggcTCTTGGCCGCGAGGCGCCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCTCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCAc                                                                      <  1:3572400/101‑1 (MQ=18)
                                gcTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCAcc                                                                     <  2:4154842/101‑1 (MQ=21)
                                 ctctTTGCCGCGATGCACCGGCGCGACGGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCa                                                                    <  2:723162/101‑1 (MQ=18)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCCCCAg                                                                   >  1:1764819/1‑101 (MQ=18)
                                  tctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAg                                                                   >  1:947129/1‑101 (MQ=21)
                                   ctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAgg                                                                  >  1:996684/1‑101 (MQ=21)
                                   ctTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAgg                                                                  >  2:4434194/1‑101 (MQ=21)
                                      tGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGaaa                                                               <  1:4295321/101‑1 (MQ=21)
                                       gCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCCGCCGATTTGCACCCGGgaaa                                                              >  2:1643011/1‑101 (MQ=17)
                                         cgcgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAg                                                            <  2:4463484/101‑1 (MQ=21)
                                           cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGcc                                                          >  2:3704511/1‑101 (MQ=21)
                                           cgATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGcc                                                          >  2:699604/1‑101 (MQ=21)
                                            gATGCACCTGCGTGACTGAGGGTTGGGCGGGGGGGCGCTGTCGCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCACGAAAAAGccc                                                         <  2:1893506/101‑1 (MQ=12)
                                             aTGCAGCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                        <  1:3719299/101‑1 (MQ=18)
                                             aTGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCg                                                        >  1:274218/1‑101 (MQ=21)
                                               gCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGc                                                      >  2:70581/1‑101 (MQ=21)
                                               gCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGc                                                      <  2:4369140/101‑1 (MQ=21)
                                                   cGGCGCGGCTGAGCGCGGGGCAGGGGGGCGCTGCCCCCCTC.TGTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                                  <  2:1562359/101‑1 (MQ=18)
                                                   cGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttt                                                  <  1:535047/101‑1 (MQ=21)
                                                    ggCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGtttt                                                 >  2:917811/1‑101 (MQ=21)
                                                     gcgcGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGttttt                                                <  2:3744865/101‑1 (MQ=21)
                                                      cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                               <  2:4144997/101‑1 (MQ=21)
                                                      cgcgACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTc                                               >  1:612552/1‑101 (MQ=21)
                                                         gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTg                                            >  2:3773198/1‑101 (MQ=21)
                                                         gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTg                                            >  2:3151128/1‑101 (MQ=21)
                                                         gACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTg                                            >  2:2307710/1‑101 (MQ=21)
                                                           cTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.t                                         >  1:1026588/1‑101 (MQ=21)
                                                            tGAGGGCGGGGCCGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.Tc                                        >  1:2562507/1‑101 (MQ=21)
                                                            tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.Tc                                        >  2:1734657/1‑101 (MQ=21)
                                                            tGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.Tc                                        >  2:1917820/1‑101 (MQ=21)
                                                                  cGGGGCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATa                                  >  2:382991/1‑101 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATAc                                 >  2:830847/1‑101 (MQ=21)
                                                                   ggggCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATAc                                 <  2:2611650/101‑1 (MQ=21)
                                                                    gggCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACg                                <  1:1862578/101‑1 (MQ=21)
                                                                     ggCAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGc                               >  2:1394614/1‑101 (MQ=21)
                                                                       cAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTTG.TCGAAATACGCCg                             >  1:1607611/1‑101 (MQ=21)
                                                                       cAGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCg                             <  1:4433644/101‑1 (MQ=21)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgg                            >  1:3063381/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgg                            >  1:4501335/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgg                            >  1:2445052/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgg                            >  2:1679965/1‑101 (MQ=25)
                                                                        aGGGGGGCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgg                            >  2:1480114/1‑101 (MQ=25)
                                                                          gggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgggg                          >  1:3713839/1‑101 (MQ=25)
                                                                          gggggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgggg                          >  1:3276711/1‑101 (MQ=25)
                                                                            gggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGTGTTTTTCCTGG.TCGAAATACGCCgggggg                        <  2:3904161/101‑1 (MQ=25)
                                                                            gggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCgggggg                        <  1:1327064/101‑1 (MQ=25)
                                                                             ggCGCTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGa                       >  1:1014990/1‑101 (MQ=25)
                                                                             cgatcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGa                       <  1:3450892/97‑1 (MQ=25)
                                                                               cgcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGt                     <  2:1468955/101‑1 (MQ=25)
                                                                                gcTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTc                    >  2:4524558/1‑101 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTcc                   <  2:3584934/101‑1 (MQ=37)
                                                                                 cTGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTcc                   >  1:3242218/1‑101 (MQ=37)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTcca                  >  2:2613606/1‑100 (MQ=37)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCg                  <  1:1405022/101‑1 (MQ=37)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCg                  <  1:4550607/101‑1 (MQ=37)
                                                                                  tGCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCg                  >  1:3926318/1‑101 (MQ=37)
                                                                                   gCCCCCCTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGaagatcg                 >  2:3450680/1‑95 (MQ=25)
                                                                                     cccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCtt               >  1:790627/1‑101 (MQ=37)
                                                                                     cccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCtt               <  1:2758921/101‑1 (MQ=37)
                                                                                     cccccTC.TTTTG.GCT..GCGACAAAATCCACCCCCCGGCGTATTTGCACCAGGAAAAACCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCtt               <  2:1796961/101‑1 (MQ=21)
                                                                                      ccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTg              >  2:2020917/1‑101 (MQ=37)
                                                                                      ccccTC.TTTTG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTg              <  1:1286321/101‑1 (MQ=37)
                                                                                           tcttttG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGa          <  2:1126535/99‑1 (MQ=37)
                                                                                            cttttG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc         <  1:917888/100‑1 (MQ=255)
                                                                                            cttttG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGAc         >  2:3340020/2‑101 (MQ=255)
                                                                                             ttttG.GCTGCGC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGACg        >  1:2266999/1‑101 (MQ=255)
                                                                                                  gggt..GC..CAAAGTTCCCCGACCGGAGTATTTGCACCA.GAAAAAACCCTGCCGTTTGCCTGGATCAAAACAGGCGGGGGGGAGTCCGCTTGCGGACGGGGGGc  <  2:4007400/98‑1 (MQ=12)
                                                                                                      ||                                                                                                    
GGGCTGAGCGTGCTGCGCTGGCCGGTGCGGCGGCTCTTTGCCGCGATGCACCGGCGCGACTGAGGGCGGGGCAGGGGGGCGCTGCCCCCCTCTTTTTG.GCT..GC..CAAAATTCACCCCCCGGCGTATTTGCACCAGGAAAAAGCCCGGCGTTTTTCCTGG.TCGAAATACGCCGGGGGGAGTCCGCTTGCGGACGGGGGGC  >  NC_014034/904834‑905031

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