Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,743 G→C 15.4% intergenic (+555/‑376) 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,7430G→C15.4% 111.9 ‑0.9 ‑0.2292intergenic (+555/‑376)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (19/26):  ref base (123/124):  total (142/150)
Fisher's exact test strand distribution p-value = 4.18e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.78e-01

Only 100 of 308 total aligned reads displayed.

CGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCC  >  NC_014034/1395644‑1395843
                                                                                                    |                                                                                                          
cgcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg..g                                                                                                         >  2:2836856/1‑101 (MQ=255)
 gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                        >  1:1235156/1‑101 (MQ=255)
 gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                        <  1:3567444/101‑1 (MQ=255)
 gcaggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                        >  1:2892064/4‑101 (MQ=255)
  cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gcg                                                                                                       >  2:1161170/1‑101 (MQ=255)
   gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gcgc                                                                                                      <  1:5029480/101‑1 (MQ=255)
    ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCa                                                                                                     >  2:3915629/1‑101 (MQ=255)
      ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCaaa                                                                                                   >  2:521171/1‑101 (MQ=255)
          tGGCTCGAAACATCTGCTTGTAGGCGGCGGGGAGTGGGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAGCTG..GCGCAAACGG..a                                                                                             <  1:334782/101‑1 (MQ=25)
           ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aa                                                                                            >  2:170361/1‑101 (MQ=255)
           ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aa                                                                                            <  1:4171172/101‑1 (MQ=255)
            ggtCGAAGCATCTGGTGGTGGGCGGCGGGCGGTTGGAGCCGGAGTGGCGCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                           <  1:3611064/99‑1 (MQ=21)
            gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                           >  1:138248/1‑101 (MQ=255)
            gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..aaa                                                                                           >  2:3603901/1‑101 (MQ=255)
             aTCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAg                                                                                          <  2:1310250/101‑1 (MQ=255)
               cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTc                                                                                        <  2:2886316/101‑1 (MQ=255)
                gAAACATCTGCATGTAGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGCCt                                                                                       >  2:2138198/1‑101 (MQ=255)
                gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                                       >  2:1940270/1‑101 (MQ=255)
                 aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                     >  1:2734176/1‑101 (MQ=255)
                 aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                     <  1:89681/101‑1 (MQ=255)
                  aaCATCTGCATGTGGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.t                                                                                   <  1:2106180/101‑1 (MQ=255)
                   aCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tc                                                                                  >  1:1928921/1‑101 (MQ=255)
                    cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcg                                                                                 >  1:2817565/1‑101 (MQ=255)
                     ctCTGCATGTGGCGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                <  2:1548907/100‑1 (MQ=255)
                      tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                               >  2:2678765/1‑101 (MQ=255)
                      tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGCCT.g.ccgtc                                                                               >  2:3679848/1‑101 (MQ=255)
                        tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                             >  2:2260291/1‑101 (MQ=255)
                           aTGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                                          >  1:3816673/1‑101 (MQ=255)
                           aTGTAGTCGGCGGGAAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                                          <  1:3124245/101‑1 (MQ=255)
                            ggTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AACGTCT.G.TCGTCGCTccc                                                                         <  1:4413511/100‑1 (MQ=255)
                             gTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                        <  1:1161287/101‑1 (MQ=255)
                             gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                        <  2:4599327/101‑1 (MQ=255)
                               cggcggcggGCAGTTCGACCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACT.GGTTTGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                      <  1:798152/98‑1 (MQ=38)
                                 tgggggggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgc                                                                    <  2:4502216/96‑1 (MQ=255)
                                   ggcggGCAGTGCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                                  <  2:3864247/101‑1 (MQ=255)
                                      gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                               >  1:3318497/1‑101 (MQ=255)
                                      gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                               >  1:4274480/1‑101 (MQ=255)
                                       ggCGGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGc                                                              <  2:3378504/101‑1 (MQ=255)
                                          aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                           <  1:323472/101‑1 (MQ=255)
                                            ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                                         <  1:4930463/101‑1 (MQ=255)
                                             tCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                        >  1:2637204/1‑101 (MQ=255)
                                             ccGAGCCGCCGTCGCCCCCCCCGCGGGGGCGTTGATCGAAACT.GGTATGGATCTG..GCCCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                        <  1:2138346/100‑1 (MQ=38)
                                              caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTgg                                                          <  1:4882990/94‑1 (MQ=2)
                                              cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                       <  2:721739/101‑1 (MQ=255)
                                               gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                        <  1:98905/93‑1 (MQ=0)
                                               aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGa                                                         <  1:3662728/95‑1 (MQ=2)
                                                agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                        >  2:1988707/6‑101 (MQ=2)
                                                 gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                      <  1:3977772/95‑1 (MQ=0)
                                                 gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.GGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                       >  2:506593/5‑101 (MQ=2)
                                                 gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                    <  2:4364384/101‑1 (MQ=255)
                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                      <  1:60664/98‑1 (MQ=2)
                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                      <  1:4809219/98‑1 (MQ=2)
                                                   ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                     <  1:2681320/99‑1 (MQ=2)
                                                    gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                 >  2:273323/1‑101 (MQ=255)
                                                    gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                 >  2:4111666/1‑101 (MQ=255)
                                                    ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                   <  2:4923984/98‑1 (MQ=2)
                                                      ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                 >  2:2289756/2‑101 (MQ=2)
                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                  <  2:4756481/101‑1 (MQ=2)
                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                               <  2:1488763/101‑1 (MQ=255)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                                 >  1:1683291/1‑100 (MQ=2)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                              >  2:4007226/1‑101 (MQ=255)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                              <  2:3744672/101‑1 (MQ=255)
                                                       gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                              >  1:1314846/1‑101 (MQ=255)
                                                        tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcac                                                <  1:2535894/101‑3 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                               <  1:1849072/101‑4 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                              <  2:4289950/101‑3 (MQ=2)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                            >  1:452418/1‑101 (MQ=255)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                            <  1:4564855/101‑1 (MQ=255)
                                                         cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGGCT.G.TCCTCGCTCCCCCCGCGGGGGCGGGGCTCGAAAGCTgtgg                                            >  2:4184789/1‑101 (MQ=25)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                              >  2:4984721/1‑97 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                             <  1:155601/101‑4 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                             <  1:1057341/101‑4 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                             >  1:378423/1‑98 (MQ=2)
                                                          gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggt                                           <  2:4565843/101‑2 (MQ=255)
                                                           cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                             <  2:4893434/101‑6 (MQ=2)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgtatgag                                           <  1:2066794/101‑8 (MQ=2)
                                                            tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                           >  1:2439499/1‑96 (MQ=2)
                                                            cccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                            <  2:2086059/100‑7 (MQ=0)
                                                             ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                           <  2:4338879/101‑8 (MQ=0)
                                                             ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                           >  2:1887599/1‑94 (MQ=0)
                                                             ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCTGCGGCGTg                                        >  1:2194961/1‑101 (MQ=255)
                                                              cccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCGCGGGGGCGTGGATCGAAACCCGCGGCGGGAg                                      >  2:3491625/1‑101 (MQ=255)
                                                                cccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGt                                     <  1:4623008/101‑1 (MQ=255)
                                                                 ccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                                    <  2:3792730/101‑1 (MQ=255)
                                                                  cgcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAgatc                                   >  2:5017295/1‑98 (MQ=255)
                                                                   gcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGc                                  <  2:3127920/101‑1 (MQ=255)
                                                                     gggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                                >  1:2319675/1‑101 (MQ=255)
                                                                          cGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                           <  1:2384191/101‑1 (MQ=255)
                                                                            tGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                         >  1:728176/1‑101 (MQ=255)
                                                                              gATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGACAGAc                       <  2:4273559/101‑1 (MQ=255)
                                                                               aTCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                      >  1:2985520/1‑101 (MQ=255)
                                                                               aTCGAAAAT.GGAATGGATCAG..GCGCAAACGG..TAAGTCT.G.TCGTCGGTCCCCCCGCGGGGGCGTGGATCGAAGCACGCGACGGGAGGTCCCAGCGATAGagg                      >  2:2443276/1‑99 (MQ=16)
                                                                                   aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCGGTGATAGACGCCGt                  >  1:4973808/1‑101 (MQ=255)
                                                                                   aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt                  >  2:3326507/1‑101 (MQ=255)
                                                                                     aCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg                <  2:2958186/101‑1 (MQ=255)
                                                                                      cT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa               >  1:4036796/1‑101 (MQ=255)
                                                                                      cT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa               >  2:2219080/1‑101 (MQ=255)
                                                                                           tATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc           >  1:3443522/1‑101 (MQ=255)
                                                                                                    g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc  >  2:4673375/1‑101 (MQ=255)
                                                                                                    g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc  <  1:3064132/101‑1 (MQ=255)
                                                                                                    |                                                                                                          
CGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCC  >  NC_014034/1395644‑1395843

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG