Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,731 +C 17.8% intergenic (+543/‑388) RCAP_rcc01306 → / → RCAP_rcc01307 Cas2 family CRISPR‑associated protein/hypothetical protein

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,395,7311.→C17.8% 154.5 ‑1.9 ‑0.0258intergenic (+543/‑388)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/27):  ref base (113/99):  total (132/126)
Fisher's exact test strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 274 total aligned reads displayed.

GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAC  >  NC_014034/1395631‑1395831
                                                                                                     |                                                                                                           
 ggccGCGCCCCCGGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                            <  1:598615/98‑1 (MQ=37)
 cggcGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                            <  1:3325631/98‑1 (MQ=255)
 cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                                                            >  1:14741/1‑101 (MQ=255)
  gTCGCCCCCCCCGCGGGGGCGTGGATCAAAACATCTGCATGTAGTGGGCGGGCAGTCCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                                                           <  1:848499/101‑1 (MQ=37)
    cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.ggga                                                                                                         >  2:3414626/1‑99 (MQ=255)
           cccGCGGGGGGGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                  <  1:2698946/101‑1 (MQ=255)
           cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                  <  1:3439648/101‑1 (MQ=255)
             cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.g..g                                                                                             <  2:2499199/101‑1 (MQ=255)
                gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCTCAGTCGATTCCCCCGCGGGGGCGTGGATCGAAACT.GTTATGGATGT.gactc                                                                                            >  2:2072258/1‑97 (MQ=37)
                gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gcgc                                                                                          >  2:3200016/1‑101 (MQ=255)
                     cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACg                                                                                     >  2:3742661/1‑101 (MQ=255)
                       tGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..a                                                                                 >  2:1629877/1‑101 (MQ=255)
                       gggATCGAACCATCTGCACGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..a                                                                                 <  2:1913152/100‑1 (MQ=255)
                        ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aa                                                                                <  1:3194024/101‑1 (MQ=255)
                         gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                               <  2:3930694/101‑1 (MQ=255)
                         gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                               >  2:1187649/1‑101 (MQ=255)
                         gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAT.GGTATCGGTCT.G..GCGCAAACGG..aag                                                                               >  2:3092418/1‑100 (MQ=255)
                             gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCt                                                                           <  2:2609382/101‑1 (MQ=255)
                             gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCt                                                                           <  2:4028524/101‑1 (MQ=255)
                                  aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                    >  1:3550556/1‑101 (MQ=255)
                                  aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgt                                                                    <  2:591493/101‑1 (MQ=255)
                                   tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.gcgtc                                                                   >  1:919872/1‑101 (MQ=255)
                                      gCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                >  2:1465308/1‑101 (MQ=255)
                                      gCATGTAGCCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                <  2:1876006/101‑1 (MQ=255)
                                       cATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTc                                                               <  1:443112/101‑1 (MQ=255)
                                        gtGTGGTCGGGGGGGAGTGGGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                              <  1:2382559/100‑1 (MQ=37)
                                        aggTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                              <  1:979659/99‑1 (MQ=255)
                                        aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                              >  2:1290095/1‑101 (MQ=255)
                                        aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                              >  2:1561747/1‑101 (MQ=255)
                                          gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                            >  1:179460/1‑101 (MQ=255)
                                               cggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgcg                                                       >  1:275299/1‑101 (MQ=255)
                                                ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                      >  1:2788767/1‑101 (MQ=255)
                                                ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                      >  2:303445/1‑101 (MQ=255)
                                                 gcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggg                                                     >  1:3469124/1‑101 (MQ=255)
                                                   gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                   >  2:1402118/1‑101 (MQ=255)
                                                   gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                   <  1:4122783/101‑1 (MQ=255)
                                                       aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                               <  1:1015688/101‑1 (MQ=255)
                                                        gtgcGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                              <  1:2952267/98‑1 (MQ=255)
                                                        ggtCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                              <  2:3451297/99‑1 (MQ=255)
                                                         ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                             >  1:1496272/1‑101 (MQ=255)
                                                         ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                             >  2:2350300/1‑101 (MQ=255)
                                                          tcaagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTg                                               >  1:3764702/9‑101 (MQ=2)
                                                           tgcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGa                                             >  2:1027987/10‑101 (MQ=0)
                                                           cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           >  1:1820771/1‑101 (MQ=255)
                                                           cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           >  2:896153/1‑101 (MQ=255)
                                                           cGACCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  2:2914618/101‑1 (MQ=255)
                                                             cgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  1:3028020/94‑1 (MQ=0)
                                                             agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCGGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  2:3673078/96‑1 (MQ=2)
                                                              gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                          <  2:486512/95‑1 (MQ=0)
                                                              gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                           >  2:4109434/5‑101 (MQ=2)
                                                              gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  1:2638359/97‑1 (MQ=2)
                                                               attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                          >  1:2833475/4‑101 (MQ=2)
                                                                gctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  1:4030793/97‑1 (MQ=2)
                                                                 tcAGTCCCGCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCGGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  2:1266001/100‑1 (MQ=2)
                                                                 gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                     <  1:3056906/101‑1 (MQ=255)
                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       <  1:1148354/101‑1 (MQ=2)
                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       <  1:2056282/101‑1 (MQ=2)
                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                       <  1:1991216/101‑1 (MQ=2)
                                                                  cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                    >  2:647281/1‑101 (MQ=255)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCta                                   >  2:2195220/1‑100 (MQ=255)
                                                                   aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                   >  1:1573329/1‑101 (MQ=255)
                                                                   aGTCGCTCCCACCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                   >  1:3659981/1‑101 (MQ=255)
                                                                    gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                  <  2:3280583/101‑1 (MQ=255)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                   >  2:1769359/1‑98 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACCtat                                  >  1:3529207/1‑99 (MQ=2)
                                                                      cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                  <  1:1116532/101‑3 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCG..CCGCTCTGGGT.AAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGGGGATCGAGAcaccc                                  >  1:693912/1‑97 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                  <  1:2735186/101‑5 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                  <  1:566754/101‑5 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACCtatg                                 >  2:3782631/1‑98 (MQ=2)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                               <  1:1629986/101‑1 (MQ=255)
                                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCgcggg                               >  2:2319886/1‑100 (MQ=255)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                 <  2:2949159/101‑6 (MQ=2)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                 <  2:3125093/101‑6 (MQ=2)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCCTGGATCGAAAcacccg                                 <  1:3695981/101‑6 (MQ=2)
                                                                        cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatga                                <  2:21273/101‑5 (MQ=2)
                                                                         tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                >  1:3393258/1‑95 (MQ=2)
                                                                         tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                               >  1:1389987/1‑96 (MQ=2)
                                                                          ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                               <  1:3149155/101‑8 (MQ=0)
                                                                          ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                            >  1:2411606/1‑101 (MQ=255)
                                                                           cccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                              <  1:4177721/101‑9 (MQ=2)
                                                                             cccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgcc                            <  1:4109695/101‑11 (MQ=2)
                                                                                  gggggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                    <  1:1967890/101‑1 (MQ=255)
                                                                                  gggggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                    >  1:3849061/1‑101 (MQ=255)
                                                                                   ggggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                   <  1:1430911/101‑1 (MQ=255)
                                                                                      gCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGa                <  2:982545/101‑1 (MQ=255)
                                                                                       cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGTTCGAAACCTGCGGCGTGAGTTGCCAGGGAt               >  1:1653899/1‑101 (MQ=255)
                                                                                       cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGGGGCGTGAGTTGCCAGTGAt               >  2:721981/1‑101 (MQ=255)
                                                                                         tGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg             <  1:2475855/101‑1 (MQ=255)
                                                                                          ggATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa            >  1:631523/1‑101 (MQ=255)
                                                                                           gATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGAc           >  2:806233/1‑101 (MQ=255)
                                                                                           gATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGAc           >  1:2993949/1‑101 (MQ=255)
                                                                                             tCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc         <  1:129554/101‑1 (MQ=255)
                                                                                             tCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc         <  1:1892411/101‑1 (MQ=255)
                                                                                                aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt      >  2:2835321/1‑101 (MQ=255)
                                                                                                aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt      >  2:1083631/1‑101 (MQ=255)
                                                                                                   cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa   >  2:3592223/1‑101 (MQ=255)
                                                                                                   cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa   <  1:1757267/101‑1 (MQ=255)
                                                                                                   cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa   <  1:22412/101‑1 (MQ=255)
                                                                                                    t.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc  <  2:2741296/101‑1 (MQ=255)
                                                                                                     |                                                                                                           
GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAC  >  NC_014034/1395631‑1395831

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