Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,670 C→G 5.3% intergenic (+482/‑449) 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,6700C→G 5.3% 18.4 ‑1.4 ‑0.0224intergenic (+482/‑449)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (8/4):  ref base (103/109):  total (111/113)
Fisher's exact test strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.34e-01

Only 100 of 247 total aligned reads displayed.

GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCT  >  NC_014034/1395570‑1395769
                                                                                                    |                                                                                                     
gCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGc                                                                                                       >  1:1071265/1‑101 (MQ=255)
  cccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCAt                                                                                                     <  1:816131/100‑1 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                    >  1:3264555/1‑101 (MQ=255)
   ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                    <  1:2323485/101‑1 (MQ=255)
    cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                   <  2:3419076/101‑1 (MQ=255)
    cccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                   <  1:4009075/101‑1 (MQ=255)
     ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCCTGTa                                                                                                  >  2:4003822/1‑101 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                 <  1:2772661/101‑1 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                 >  1:157611/1‑101 (MQ=255)
      cccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAg                                                                                                 <  1:618687/101‑1 (MQ=255)
       ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                <  1:3693159/101‑1 (MQ=255)
        cgcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTc                                                                                               <  1:228663/101‑1 (MQ=255)
         gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                              >  1:4144220/1‑101 (MQ=255)
         gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                              <  1:482649/101‑1 (MQ=255)
              ggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggcggc                                                                                         <  1:2572949/101‑2 (MQ=255)
               gCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.Gc                                                                                       <  2:3340658/101‑1 (MQ=255)
                cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCa                                                                                      >  1:3908300/1‑101 (MQ=255)
                 gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAg                                                                                     >  2:2535801/1‑101 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGCG.GCact                                                                                    >  1:453205/1‑99 (MQ=255)
                  tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGt                                                                                    >  1:3747002/1‑101 (MQ=255)
                   ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGtt                                                                                   >  2:2352148/1‑101 (MQ=255)
                    gATCGAAACCACCTTGCGGATTGGCACTTTCGCGGTGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTc                                                                                  <  2:527116/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCg                                                                                 >  1:674239/1‑101 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCg                                                                                 <  2:4143959/101‑1 (MQ=255)
                     aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGAGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCg                                                                                 >  2:1146289/1‑101 (MQ=255)
                      gcGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGa                                                                                <  2:1487304/100‑1 (MQ=255)
                       cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.g                                                                              >  1:1802406/1‑101 (MQ=255)
                        gAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.Gc                                                                             >  1:3982282/1‑101 (MQ=255)
                         aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.Gcc                                                                            <  1:3050159/101‑1 (MQ=255)
                          aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCg                                                                           <  1:741645/101‑1 (MQ=255)
                          aaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCg                                                                           >  2:3735689/1‑101 (MQ=255)
                            cccccTTGCGGATTGTCCCTTTCGCGGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCa                                                                         <  1:1325651/98‑1 (MQ=255)
                            ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCa                                                                         <  2:2897632/101‑1 (MQ=255)
                            ccaccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCa                                                                         <  2:194520/101‑1 (MQ=255)
                             caccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAg                                                                        >  2:621506/1‑101 (MQ=255)
                               ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTc                                                                      <  1:2520057/101‑1 (MQ=255)
                                cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCg                                                                     >  1:578040/1‑101 (MQ=255)
                                 ttGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGc                                                                    <  2:308217/101‑1 (MQ=255)
                                  tGCGGATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCt                                                                   <  1:1428999/101‑1 (MQ=255)
                                  tGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCt                                                                   <  1:3440399/101‑1 (MQ=255)
                                   ccGGATTGTCCCTTGCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTc                                                                  <  2:4190437/100‑1 (MQ=255)
                                      gATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTcccc                                                               <  1:2058529/101‑1 (MQ=255)
                                       aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTccccc                                                              <  1:3627139/101‑1 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCTA.GCCGCAGTCGCTCCCCCcg                                                            <  1:598947/101‑1 (MQ=255)
                                         tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                            <  1:3636560/101‑1 (MQ=255)
                                            cACTTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                         <  1:2981603/101‑1 (MQ=255)
                                             aCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                        >  1:1931550/1‑101 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                      >  2:712163/1‑101 (MQ=255)
                                               tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTACCCCCGCggggg                                                      >  2:2015694/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)
                                                    gggggCGACGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  1:1265649/100‑1 (MQ=255)
                                                    cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  2:2745645/101‑1 (MQ=255)
                                                     gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                 <  2:489962/93‑1 (MQ=2)
                                                     ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGGGTggg                                                >  2:1004670/1‑100 (MQ=255)
                                                      tgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGa                                                >  2:1717133/8‑101 (MQ=0)
                                                      gggGGAGGCGCCGCGCCCCCCGCGGGGGCGTGGATCGAAACTTCTGCATCGGGTGGGGGG.GCGGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:3506204/101‑1 (MQ=21)
                                                      gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:896881/101‑1 (MQ=255)
                                                      gggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                               <  2:3264323/101‑1 (MQ=255)
                                                        gCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  1:1096445/101‑1 (MQ=255)
                                                        gCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                             <  1:2359705/101‑1 (MQ=255)
                                                        ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                              <  2:3976626/96‑1 (MQ=2)
                                                          gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                           >  1:1869936/1‑101 (MQ=255)
                                                           aaGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  2:4071815/1‑101 (MQ=255)
                                                            tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  2:2212590/2‑101 (MQ=2)
                                                            tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  1:2140397/2‑101 (MQ=2)
                                                            aGCGTCGGCCCCCCCGGGGGGGCGTGGATCGAAACAGCTGCATGTAGTGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                         <  2:2459554/101‑1 (MQ=37)
                                                              cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                       >  1:311163/1‑101 (MQ=255)
                                                              cGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                       <  2:1576146/101‑1 (MQ=255)
                                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                      <  2:3788204/101‑1 (MQ=255)
                                                                gcGCCCCCCCCGGCGGGGGGTGGATCAACACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                     <  1:32959/100‑1 (MQ=25)
                                                                   ccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATg                                  <  1:217975/99‑1 (MQ=255)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  2:356006/1‑101 (MQ=255)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  2:2939734/1‑101 (MQ=255)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATgg                                 >  1:4068018/1‑101 (MQ=255)
                                                                     ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGa                                >  1:3098573/1‑101 (MQ=255)
                                                                           gcgGGGGCGTGGATCGAAACATCTGCAGGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                          <  2:1837398/101‑1 (MQ=255)
                                                                            cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgcg                         >  1:2722511/1‑101 (MQ=255)
                                                                              ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCa                       <  2:3908662/101‑1 (MQ=255)
                                                                                ggCGTGGATCGAAACATCTGCATGTAGTCGGGGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                     >  1:3951700/1‑101 (MQ=255)
                                                                                 gCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                    >  1:3368259/1‑101 (MQ=255)
                                                                                  cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACg                   >  1:2661199/1‑101 (MQ=255)
                                                                                   gtagATGGAAACATATGCATGTGGCCGGGGG.GGAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  <  1:1193797/98‑1 (MQ=25)
                                                                                   gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  >  1:2037351/1‑101 (MQ=255)
                                                                                   gTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACgg                  >  2:1547835/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 >  2:829063/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 >  2:229294/1‑101 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 <  2:1743662/101‑1 (MQ=255)
                                                                                    tGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                 >  2:1281276/1‑101 (MQ=255)
                                                                                       aTCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg              >  1:80252/1‑101 (MQ=255)
                                                                                        tCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGt             <  2:321850/101‑1 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg          >  2:1691194/1‑101 (MQ=255)
                                                                                           aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTg          >  2:1537956/1‑101 (MQ=255)
                                                                                             aCACCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtc        <  2:2702594/101‑1 (MQ=255)
                                                                                              cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcg       >  2:4056125/1‑101 (MQ=255)
                                                                                               agcTGCATGGAGTCGGCGG.GGAGTTCGA.CCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt      <  2:3044277/99‑1 (MQ=255)
                                                                                                tCTGGAGGTGGTGGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc     <  2:4341927/101‑1 (MQ=37)
                                                                                                tCTGCATGTGGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc     <  1:2672198/101‑1 (MQ=255)
                                                                                                   gCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCt  >  1:102695/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
GCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCT  >  NC_014034/1395570‑1395769

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