Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 5.3% intergenic (+510/‑421) 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,6980A→. 5.3% 27.0 ‑1.4 0.0227intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/4):  ref base (101/113):  total (109/118)
Fisher's exact test strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 248 total aligned reads displayed.

CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798
                                                                                                      |                                                                                                       
cccccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCa                                                                                                         <  2:490829/98‑1 (MQ=255)
ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCa                                                                                                         <  2:194520/101‑1 (MQ=255)
  accTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGt                                                                                                       <  1:1040768/101‑1 (MQ=255)
  accTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGt                                                                                                       >  1:3439271/1‑101 (MQ=255)
   ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTc                                                                                                      <  1:2441117/101‑1 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  1:578040/1‑101 (MQ=255)
     ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGc                                                                                                    <  2:3519710/101‑1 (MQ=255)
       gCGGATTGTCGCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  <  1:754805/101‑1 (MQ=255)
       ccGGATTGTCCCTTGCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  <  2:4190437/100‑1 (MQ=255)
        cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                                 <  2:84113/101‑1 (MQ=255)
        cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcc                                                                                                 >  1:4014569/1‑101 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               <  1:427003/101‑1 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               <  1:2058529/101‑1 (MQ=255)
          gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                                                               <  2:3064234/101‑1 (MQ=255)
             tGTCACTTCCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                            <  2:3835910/101‑1 (MQ=255)
               tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                                                          >  1:3501222/1‑101 (MQ=255)
               tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                                                          <  1:500662/101‑1 (MQ=255)
                 aCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                                        >  1:1931550/1‑101 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                                                      >  2:712163/1‑101 (MQ=255)
                   tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCCCggggg                                                                                      >  1:1964492/1‑101 (MQ=255)
                    ttCGCGGGTCGAAGCGTCGCCCCCCCCGCGGGGACGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGc                                                                                     <  1:3852222/101‑1 (MQ=255)
                      cgcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCgg                                                                                   >  2:4280485/1‑100 (MQ=255)
                       tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                                                   >  2:124383/11‑101 (MQ=2)
                       gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                                                  <  2:2236862/101‑1 (MQ=255)
                        cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  2:2745645/101‑1 (MQ=255)
                        cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  1:3678340/101‑1 (MQ=255)
                         ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGGGTggg                                                                                >  2:1004670/1‑100 (MQ=255)
                          tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                                >  1:2761306/8‑101 (MQ=0)
                          tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                                                >  2:1717133/8‑101 (MQ=0)
                          gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  2:3974134/101‑1 (MQ=255)
                          gggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  2:3264323/101‑1 (MQ=255)
                            gCGAGGCGTCCCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  2:1581173/101‑1 (MQ=255)
                            gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:1096445/101‑1 (MQ=255)
                            gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  1:3910691/101‑1 (MQ=255)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCCGCTCCA.GCCGCAGTCGCTCCCCCCCCGGGGGCGGGGATCGaa                                                                           >  1:4089719/1‑101 (MQ=37)
                              gAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           <  1:2983183/101‑1 (MQ=255)
                                 gCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                         >  1:217397/1‑101 (MQ=2)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      >  1:311163/1‑101 (MQ=255)
                                    gcGCCCCCCCCGGCGGGGGGTGGATCAACACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                    <  1:32959/100‑1 (MQ=25)
                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTAt                                                                  >  2:2973063/1‑101 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 >  2:3039144/1‑101 (MQ=255)
                                        tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TGttcgg                                                                >  2:3422239/1‑97 (MQ=255)
                                            cccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  1:2566959/101‑1 (MQ=255)
                                                cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGgcg                                                        >  1:2722511/1‑101 (MQ=255)
                                                 gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACT.GGTATGGATCTGgcgc                                                       >  2:562734/1‑101 (MQ=255)
                                                   gggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaa                                                     <  2:282795/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    <  2:1531056/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    <  1:3362370/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    <  1:621574/101‑1 (MQ=255)
                                                      cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACg                                                  >  1:2661199/1‑101 (MQ=255)
                                                       gtagATGGAAACATATGCATGTGGCCGGGGG.GGAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 <  1:1193797/98‑1 (MQ=25)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  2:1450670/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  1:3768994/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  2:1547835/1‑101 (MQ=255)
                                                       gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACgg                                                 >  1:2037351/1‑101 (MQ=255)
                                                        tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              >  2:1281276/1‑101 (MQ=255)
                                                        tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..a                                              <  2:1743662/101‑1 (MQ=255)
                                                           aTCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAg                                           >  1:2134865/1‑101 (MQ=255)
                                                            tCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGt                                          <  2:321850/101‑1 (MQ=255)
                                                             cGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTAGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTc                                         <  1:594622/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  1:1306888/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  1:1438693/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        <  2:1869812/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  2:3981513/1‑101 (MQ=255)
                                                               aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       >  2:1691194/1‑101 (MQ=255)
                                                                 aCACCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtc                                     <  2:2702594/101‑1 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    >  2:4056125/1‑101 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    <  1:3632364/101‑1 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCAGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    >  1:3958474/1‑101 (MQ=255)
                                                                   ctcggCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  2:4141045/97‑1 (MQ=255)
                                                                    tCTGGAGGTGGTGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  <  2:4341927/101‑1 (MQ=37)
                                                                    tCTGCATGTGGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  <  1:2672198/101‑1 (MQ=255)
                                                                    tCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  <  1:3172250/101‑1 (MQ=255)
                                                                    tCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgtc                                  <  2:2453155/101‑1 (MQ=255)
                                                                      tGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGc                                >  1:1766536/1‑101 (MQ=255)
                                                                       tgaTGTAGGGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGAGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCt                               <  1:1000128/99‑1 (MQ=255)
                                                                       gcctGTAGGCGGCGG.GCAGTGCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCt                               <  1:2262801/98‑1 (MQ=255)
                                                                          tGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccc                            <  1:3039365/101‑1 (MQ=255)
                                                                           gTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccc                           <  2:2937660/101‑1 (MQ=255)
                                                                             aGTCGGGGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  1:1415543/101‑1 (MQ=255)
                                                                             aGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         >  1:2850836/1‑101 (MQ=255)
                                                                              gTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcg                        >  2:819203/1‑101 (MQ=255)
                                                                                cggcgg.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcgcg                      >  2:3784859/1‑101 (MQ=255)
                                                                                  gcgg.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggg                    <  2:311128/101‑1 (MQ=255)
                                                                                    gg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggggg                  >  1:816004/1‑101 (MQ=255)
                                                                                       gCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCg                <  2:4077986/101‑1 (MQ=255)
                                                                                          gtgcGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             <  2:1927522/98‑1 (MQ=255)
                                                                                            tggA.GCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  1:2309088/99‑1 (MQ=255)
                                                                                             cGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCCCGGGGGCGTGGATc          >  1:1668730/1‑101 (MQ=255)
                                                                                              caagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGACGGTGGTCGCTCCCCCCGCGGGGGCGTgg             >  2:1136513/8‑101 (MQ=2)
                                                                                                agattgaGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGAt           <  2:25814/95‑1 (MQ=2)
                                                                                                agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGAt           <  2:488347/96‑1 (MQ=2)
                                                                                                 gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  1:4329476/101‑1 (MQ=255)
                                                                                                 gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:445303/101‑1 (MQ=255)
                                                                                                 gCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:883851/101‑1 (MQ=255)
                                                                                                   cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     >  2:1798745/1‑101 (MQ=255)
                                                                                                   cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     <  1:406260/101‑1 (MQ=255)
                                                                                                    tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  2:2010199/100‑1 (MQ=2)
                                                                                                     cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt   >  1:1203010/1‑101 (MQ=255)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg  <  1:1383700/101‑1 (MQ=255)
                                                                                                      |                                                                                                       
CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798

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