Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,734 T→C 16.0% intergenic (+546/‑385) 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,7340T→C16.0% 112.3 ‑0.7 ‑1.6302intergenic (+546/‑385)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (25/20):  ref base (117/118):  total (154/148)
Fisher's exact test strand distribution p-value = 5.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 2.60e-02

Only 100 of 313 total aligned reads displayed.

TCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTC  >  NC_014034/1395634‑1395834
                                                                                                     |                                                                                                          
cccacgcccccGCGGGGGCGGGGATGGAAACATCTGCATGTGGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                                                            <  1:4345951/95‑1 (MQ=37)
  gccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTAt                                                                                                          <  2:1111510/98‑1 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTAt                                                                                                          >  2:4451473/1‑101 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                        >  1:2558229/1‑101 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                        >  2:333000/1‑101 (MQ=255)
    cccccccGGGGGGGCTTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                        <  1:4133412/100‑1 (MQ=255)
          cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg..g                                                                                                <  2:3622294/101‑1 (MQ=255)
           gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                               <  2:2848268/101‑1 (MQ=255)
           gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                               <  1:3921331/101‑1 (MQ=255)
           gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                               <  1:2287525/101‑1 (MQ=255)
            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGGTCCCCCCGCGGGTGCGTGGATCGAAACT.GGTATGGATCTG..gcg                                                                                              >  1:303260/1‑101 (MQ=255)
            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gcg                                                                                              <  2:4578539/101‑1 (MQ=255)
            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gcg                                                                                              >  1:74763/1‑101 (MQ=255)
            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gcg                                                                                              >  1:1083018/1‑101 (MQ=255)
              ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCa                                                                                            <  1:1247097/101‑1 (MQ=255)
                ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCaaa                                                                                          <  2:3854951/101‑1 (MQ=255)
                  cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACg                                                                                        >  2:3679177/1‑101 (MQ=255)
                  cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACg                                                                                        <  2:873761/101‑1 (MQ=255)
                   gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACgg                                                                                       >  1:2239688/1‑101 (MQ=255)
                   gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACgg                                                                                       >  1:1524427/1‑101 (MQ=255)
                      gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                  >  2:3093844/1‑101 (MQ=255)
                      gATCGAAACATCCGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                  <  1:2834193/101‑1 (MQ=255)
                          gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                              >  2:1822788/1‑101 (MQ=255)
                          gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                              >  2:945078/1‑101 (MQ=255)
                          gAAACATCTGCACGTAGGCGGCGGGCAGTGCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                              <  2:1484410/101‑1 (MQ=255)
                           aaaCATCTGCATGTAGGCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                            <  1:874158/101‑1 (MQ=255)
                            aaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.t                                                                          <  1:4085980/101‑1 (MQ=255)
                             aCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tc                                                                         <  1:2466967/101‑1 (MQ=255)
                               aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcag                                                                       >  1:713618/1‑99 (MQ=255)
                               aTCTGCAGGGAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                       <  1:1693849/101‑1 (MQ=255)
                                 cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                     >  1:4527098/1‑101 (MQ=255)
                                 cTGCATGTAGGCGGCGGGCAGTCCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                     <  2:1982267/101‑1 (MQ=255)
                                  tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                    >  1:4704638/1‑101 (MQ=255)
                                  tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                    <  2:4083736/101‑1 (MQ=255)
                                    cATGTAGTCGGCTGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTc                                                                  >  1:4575174/1‑101 (MQ=255)
                                    cATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTc                                                                  <  1:333034/101‑1 (MQ=255)
                                      tGTGGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccc                                                                <  1:1008536/101‑1 (MQ=255)
                                      tGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccc                                                                >  2:680513/1‑101 (MQ=255)
                                       gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                               >  2:3359987/1‑101 (MQ=255)
                                       gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                               <  1:3277560/101‑1 (MQ=255)
                                          gTCGGCGGGCATTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcg                                                            <  1:1788660/101‑1 (MQ=255)
                                          gTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcg                                                            >  1:4006729/1‑101 (MQ=255)
                                           tCGGCGGCCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgc                                                           <  1:163999/101‑1 (MQ=255)
                                             ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                         <  2:2652036/101‑1 (MQ=255)
                                               cggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgggg                                                       >  1:3352971/1‑101 (MQ=255)
                                                gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                      >  2:1097349/1‑101 (MQ=255)
                                                gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGAATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                      >  2:3932533/1‑101 (MQ=255)
                                                 ggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGc                                                     <  1:2345679/101‑1 (MQ=255)
                                                  gCATTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCg                                                    <  1:4289933/101‑1 (MQ=255)
                                                      ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  1:3178278/101‑1 (MQ=255)
                                                       tggAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:2940689/99‑1 (MQ=255)
                                                        cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                              <  2:4006487/101‑1 (MQ=255)
                                                        cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                              <  1:3745380/101‑1 (MQ=255)
                                                         aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  1:1176314/95‑1 (MQ=2)
                                                          aGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                            <  1:3386006/101‑1 (MQ=255)
                                                           gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  1:3651983/95‑1 (MQ=0)
                                                           gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                              <  1:632628/97‑1 (MQ=2)
                                                            attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  2:2278787/98‑1 (MQ=2)
                                                            attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  1:2960804/98‑1 (MQ=2)
                                                             gctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           >  2:1940721/5‑101 (MQ=2)
                                                              tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           >  2:2237326/2‑101 (MQ=2)
                                                              gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                        <  2:4133997/101‑1 (MQ=255)
                                                               cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACCt                                       >  2:2350673/1‑101 (MQ=255)
                                                               cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                       <  2:3525284/101‑1 (MQ=255)
                                                                aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         >  2:4616800/1‑101 (MQ=2)
                                                                aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         >  2:2697263/1‑101 (MQ=2)
                                                                   cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                      >  1:2873885/1‑98 (MQ=2)
                                                                   cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                     >  1:3610317/1‑99 (MQ=2)
                                                                   cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                     >  1:2039795/1‑99 (MQ=2)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                     <  1:1569710/101‑5 (MQ=2)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                    <  2:1695853/101‑4 (MQ=2)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                    >  2:1155676/1‑98 (MQ=2)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                  >  1:4606182/1‑101 (MQ=255)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                  >  2:2501717/1‑101 (MQ=255)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAATGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                    <  1:955707/101‑6 (MQ=2)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatga                                   >  2:2164446/1‑97 (MQ=2)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggcg                                 >  1:2727116/1‑101 (MQ=255)
                                                                      tCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGt                                >  1:4366109/1‑101 (MQ=255)
                                                                       ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                  >  2:3118724/1‑94 (MQ=0)
                                                                       ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                  >  2:2576631/1‑94 (MQ=0)
                                                                       ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                               >  2:1333160/1‑101 (MQ=255)
                                                                       ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                               >  2:2779105/1‑101 (MQ=255)
                                                                         ccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgc                                >  2:1068479/1‑92 (MQ=2)
                                                                            cgcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTg                          >  1:513622/1‑101 (MQ=255)
                                                                               gggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                       >  1:2453674/1‑101 (MQ=255)
                                                                               gggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGGTGccg                       >  2:3272562/1‑100 (MQ=255)
                                                                                   gCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGa                   >  2:1571275/1‑101 (MQ=255)
                                                                                     gTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..CAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                 >  2:2409430/1‑101 (MQ=255)
                                                                                     gTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                 >  2:49737/1‑101 (MQ=255)
                                                                                       ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa               >  1:4399762/1‑101 (MQ=255)
                                                                                           cGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGcc           <  1:3511853/101‑1 (MQ=255)
                                                                                            gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg          <  1:999017/101‑1 (MQ=255)
                                                                                            gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGGGAGTTGCCAGTGATAGACGCCg          >  2:776075/1‑101 (MQ=255)
                                                                                             aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt         >  2:2997030/1‑101 (MQ=255)
                                                                                              aaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTg        >  2:1400372/1‑101 (MQ=255)
                                                                                               aCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg       <  2:4302151/101‑1 (MQ=255)
                                                                                                 t.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc     >  2:4102548/1‑101 (MQ=255)
                                                                                                 t.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc     <  2:2441124/101‑1 (MQ=255)
                                                                                                 t.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc     >  1:1942408/1‑101 (MQ=255)
                                                                                                     tATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc  <  2:3740830/101‑1 (MQ=255)
                                                                                                     |                                                                                                          
TCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTC  >  NC_014034/1395634‑1395834

Base quality scores:  ATCG  < 3 ≤  ATCG  < 20 ≤  ATCG  < 28 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG