Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,737 Δ1 bp 7.5% intergenic (+549/‑382) 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,7370G→. 7.5% 64.4 ‑0.7 ‑0.1305intergenic (+549/‑382)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (14/9):  ref base (149/133):  total (163/142)
Fisher's exact test strand distribution p-value = 5.19e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.63e-01

Only 100 of 315 total aligned reads displayed.

TCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCT  >  NC_014034/1395638‑1395837
                                                                                                    |                                                                                                          
cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                           <  2:92621/100‑1 (MQ=255)
 ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGa                                                                                                          >  2:2570285/1‑101 (MQ=255)
  cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                         >  1:3163244/1‑101 (MQ=255)
   ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                                                                        >  1:373867/1‑101 (MQ=255)
     ccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg                                                                                                      >  1:2046872/1‑101 (MQ=255)
       gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..gc                                                                                                  <  1:3068577/101‑1 (MQ=255)
       gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACT.GGTAGGGATCTG..gc                                                                                                  >  2:374939/1‑101 (MQ=255)
          ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCa                                                                                               >  2:2883898/1‑101 (MQ=255)
           gggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCaa                                                                                              >  2:3653403/1‑101 (MQ=255)
            ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCaaa                                                                                             >  2:2667554/1‑101 (MQ=255)
             gCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAAc                                                                                            <  1:3781842/101‑1 (MQ=255)
              cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACg                                                                                           >  1:4098536/1‑101 (MQ=255)
              cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACg                                                                                           >  1:2124910/1‑101 (MQ=255)
              cGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACg                                                                                           >  2:2753365/1‑101 (MQ=255)
               gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACgg                                                                                          >  2:4544386/1‑101 (MQ=255)
                tGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..a                                                                                       >  2:3619179/1‑101 (MQ=255)
                     cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTc                                                                                  >  2:2500292/1‑101 (MQ=255)
                      gAAACATCTGCATGTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                                 <  1:2078454/101‑1 (MQ=255)
                      gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCt                                                                                 >  1:1538630/1‑101 (MQ=255)
                       aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                               >  2:3580190/1‑101 (MQ=255)
                          cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTcagatc                                                                             >  1:2410716/1‑96 (MQ=255)
                           ctCTGGATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                          <  2:761454/100‑1 (MQ=255)
                           aTCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                          <  2:2213704/101‑1 (MQ=255)
                            tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                         <  2:4617803/101‑1 (MQ=255)
                             cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                        <  2:4412411/101‑1 (MQ=255)
                              tGCATGTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                       <  2:4596815/101‑1 (MQ=255)
                                   gTAGTCGTCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                  >  2:3977951/1‑101 (MQ=255)
                                    tAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccccc                                                                 >  1:4517277/1‑101 (MQ=255)
                                     aggcggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                <  2:1188966/98‑1 (MQ=255)
                                     aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                <  1:4613562/101‑1 (MQ=255)
                                     aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                <  1:1423169/101‑1 (MQ=255)
                                      gTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcg                                                               >  1:2490861/1‑101 (MQ=255)
                                       tCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgc                                                              <  2:163898/101‑1 (MQ=255)
                                          gcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggg                                                           >  1:1304446/1‑101 (MQ=255)
                                          gcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GATATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggg                                                           >  2:2951661/1‑101 (MQ=255)
                                           ggggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgggg                                                          <  1:1489895/100‑1 (MQ=255)
                                            gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTAT.G.TCGTCGCTCCCCCCGCggggg                                                         >  1:2703832/1‑101 (MQ=255)
                                             ggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATGTG..GCGCAAACGG..AAAGTCT.G.TCCTCGCTCCCCCCGCGGGGGc                                                        >  2:3451924/1‑101 (MQ=255)
                                             ggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGc                                                        >  2:4349142/1‑101 (MQ=255)
                                              gCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCg                                                       <  2:1530540/101‑1 (MQ=255)
                                              gCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCg                                                       >  1:1737181/1‑101 (MQ=255)
                                                atttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                     <  2:3239287/99‑1 (MQ=255)
                                                   tcaagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTg                                                     >  2:605887/9‑101 (MQ=2)
                                                    tgcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                   >  1:4388154/10‑101 (MQ=0)
                                                    caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTgg                                                    <  1:2860000/94‑1 (MQ=2)
                                                    caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTgg                                                    <  2:3194019/94‑1 (MQ=2)
                                                    caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTgg                                                    <  1:3209641/94‑1 (MQ=2)
                                                    cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                 <  1:479351/101‑1 (MQ=255)
                                                     gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                >  1:1950425/1‑101 (MQ=255)
                                                     gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                >  2:1354112/1‑101 (MQ=255)
                                                      ggCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                               <  2:2121320/100‑1 (MQ=255)
                                                      cgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                 <  2:3489250/94‑1 (MQ=0)
                                                      cgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                 <  1:4289853/94‑1 (MQ=0)
                                                       gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                >  2:4185670/7‑101 (MQ=0)
                                                        ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                             <  1:3945996/101‑1 (MQ=255)
                                                        ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGGCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                             <  2:1446838/101‑1 (MQ=255)
                                                        attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                <  1:3400467/98‑1 (MQ=2)
                                                         ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                               >  1:2635635/3‑101 (MQ=2)
                                                         gctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                              >  2:2798113/5‑101 (MQ=2)
                                                         cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaaa                                            >  1:3542527/1‑100 (MQ=255)
                                                          gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                           <  1:4137029/101‑1 (MQ=255)
                                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                           <  1:1767462/100‑1 (MQ=2)
                                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAAcc                                           >  1:2750144/2‑101 (MQ=2)
                                                            aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                         <  1:943617/101‑1 (MQ=255)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                           <  2:1076338/101‑2 (MQ=2)
                                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                          >  2:2458298/1‑101 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGACG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                         >  1:566111/1‑98 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGNCGTGGATCGAAACCtat                                        >  2:1960121/1‑99 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                        >  2:3713949/1‑99 (MQ=2)
                                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGGGGCTCGAAAcccat                                        >  2:3882/1‑98 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                        >  1:416133/1‑97 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgcatg                                       <  2:3949448/101‑6 (MQ=2)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                     <  1:814475/101‑1 (MQ=255)
                                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                     <  2:2184236/101‑1 (MQ=255)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                       <  1:4326623/101‑6 (MQ=2)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatga                                      >  2:733649/1‑97 (MQ=2)
                                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggcg                                    <  2:501049/101‑1 (MQ=255)
                                                                  tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                     >  1:4630914/1‑96 (MQ=2)
                                                                       ccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                              >  1:2558966/1‑101 (MQ=255)
                                                                          cgGGGGAGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGcc                           <  1:745423/101‑1 (MQ=255)
                                                                            ggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                         <  1:1906833/101‑1 (MQ=255)
                                                                            ggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                         >  2:2027370/1‑101 (MQ=255)
                                                                             gggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGt                        <  2:1737062/101‑1 (MQ=255)
                                                                             gggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGt                        <  1:3635752/101‑1 (MQ=255)
                                                                                 gTGGATCGAAACT.GGTATGGATCTG..GCGCAAACCT..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                    <  1:4005472/101‑1 (MQ=255)
                                                                                  tGGATCGAAACT.GGTATGGATCTG..GCGCAAACGA..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                   <  1:4037410/101‑1 (MQ=255)
                                                                                   ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                  >  2:330087/1‑101 (MQ=255)
                                                                                   ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                  >  2:2343658/1‑101 (MQ=255)
                                                                                     aTCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                <  1:3305823/101‑1 (MQ=255)
                                                                                       cGAAANT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGcc              <  1:1402384/101‑1 (MQ=255)
                                                                                       cGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGcc              >  2:1788846/1‑101 (MQ=255)
                                                                                       cGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGcc              >  1:1031767/1‑101 (MQ=255)
                                                                                        gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg             >  1:2809837/1‑101 (MQ=255)
                                                                                          aaCT.GGTATGGAGCTG..GGGCAAACGG..AAAGTCT.G.TCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTg           <  2:4098250/101‑1 (MQ=255)
                                                                                            cT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa         <  1:1083238/101‑1 (MQ=255)
                                                                                             t.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc        <  2:4258161/101‑1 (MQ=255)
                                                                                               gatATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACg       <  1:2951890/99‑1 (MQ=255)
                                                                                                  aTGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCg    >  1:2348573/1‑101 (MQ=255)
                                                                                                   tGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGc   <  1:3494717/101‑1 (MQ=255)
                                                                                                    ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCt  >  2:1117856/1‑101 (MQ=255)
                                                                                                    |                                                                                                          
TCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCT  >  NC_014034/1395638‑1395837

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