Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,812 C→A 4.6% intergenic (+624/‑307) 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,8120C→A 4.6% 12.1 ‑0.3 ‑0.0204intergenic (+624/‑307)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (6/4):  ref base (98/96):  total (104/100)
Fisher's exact test strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.45e-01

Only 100 of 215 total aligned reads displayed.

CGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGC  >  NC_014034/1395712‑1395912
                                                                                                     |                                                                                                      
cgGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCATACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGcc                                                                                                        <  1:3572500/101‑1 (MQ=255)
 gggggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCa                                                                                                       >  1:3849061/1‑101 (MQ=255)
    ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                    >  2:1230011/1‑101 (MQ=255)
    ggCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTg                                                                                                    >  2:2328922/1‑101 (MQ=255)
      cGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGGG.GCGTGAGTTGCCAGTGAt                                                                                                  >  2:721981/1‑101 (MQ=255)
       gTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATa                                                                                                 <  2:1613234/101‑1 (MQ=255)
         ggATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGa                                                                                               >  1:631523/1‑101 (MQ=255)
          gATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.c                                                                                             >  2:806233/1‑101 (MQ=255)
           aTCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.Cg                                                                                            >  1:2725732/1‑101 (MQ=255)
            tCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                                                                           <  1:129554/101‑1 (MQ=255)
            tCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                                                                           <  1:1892411/101‑1 (MQ=255)
            tCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                                                                           >  2:3328829/1‑101 (MQ=255)
            tCGAAACCTGTTTGGGTCTGGGGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGc                                                                                           <  1:2756509/101‑1 (MQ=37)
             cGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGcc                                                                                          >  1:368324/1‑101 (MQ=255)
             cGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGcc                                                                                          >  1:1645150/1‑101 (MQ=255)
             cGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGcc                                                                                          <  2:3114734/101‑1 (MQ=255)
               aaaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGt                                                                                        >  2:1083631/1‑101 (MQ=255)
                gaCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTg                                                                                       >  2:1531239/2‑101 (MQ=255)
                 aCTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGTCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTgg                                                                                      <  2:1403057/101‑1 (MQ=255)
                  cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                     <  1:22412/101‑1 (MQ=255)
                  cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                     <  1:1757267/101‑1 (MQ=255)
                  cTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGa                                                                                     >  2:3422542/1‑101 (MQ=255)
                   tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                   <  2:2741296/101‑1 (MQ=255)
                   tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                   <  1:3200255/101‑1 (MQ=255)
                   tGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.c                                                                                   >  1:3450447/1‑101 (MQ=255)
                    ggTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.Cg                                                                                  <  1:1561851/101‑1 (MQ=255)
                      tATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTc                                                                                <  2:890403/101‑1 (MQ=255)
                       aTGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCg                                                                               <  1:1187743/101‑1 (MQ=255)
                       aTGGATCTGGCGCAAACGGAAAGGCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGCGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCg                                                                               >  2:298696/1‑101 (MQ=255)
                        tGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGc                                                                              <  2:965906/101‑1 (MQ=255)
                         ggATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCt                                                                             >  1:1877823/1‑101 (MQ=255)
                           aTCTGGCGCAAACGGAAAGTCTGTCGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcc                                                                           >  1:1244178/1‑101 (MQ=255)
                           aTCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcc                                                                           >  2:1937197/1‑101 (MQ=255)
                             cTGGCGCAAACGGAAAGTCTGTCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccc                                                                         <  1:1424778/101‑1 (MQ=255)
                               ggCGCAAACGGAAAGTCTGTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                       <  1:3389196/101‑1 (MQ=255)
                               ggCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTcccccc                                                                       >  1:170315/1‑101 (MQ=255)
                                 cgcAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCcgc                                                                     >  1:1110315/1‑101 (MQ=255)
                                   cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                   >  2:2193587/1‑101 (MQ=255)
                                   cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                   >  2:3023044/1‑101 (MQ=255)
                                   cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                   >  2:3767055/1‑101 (MQ=255)
                                   cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                   <  1:599329/101‑1 (MQ=255)
                                   cAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgg                                                                   >  1:2032776/1‑101 (MQ=255)
                                    aaaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggg                                                                  >  1:4081920/1‑101 (MQ=255)
                                     aaCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCgggg                                                                 >  2:3250835/1‑101 (MQ=255)
                                      aCGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCggggg                                                                <  2:2411430/101‑1 (MQ=255)
                                       cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                               <  1:2329090/101‑1 (MQ=255)
                                       cGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGc                                                               >  1:2588767/1‑101 (MQ=255)
                                        ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                              <  2:1783277/101‑1 (MQ=255)
                                        ggAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCg                                                              >  1:3215379/1‑101 (MQ=255)
                                         cggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCg                                                              >  2:2423160/11‑101 (MQ=0)
                                          ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                             >  1:1044080/10‑101 (MQ=0)
                                          ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                             >  2:3633167/10‑101 (MQ=0)
                                          ggacccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGGGGCGt                                                             >  2:1325824/10‑101 (MQ=0)
                                           aaGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTgg                                                           <  1:298741/101‑1 (MQ=255)
                                            aGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGa                                                          >  1:1807410/1‑101 (MQ=255)
                                            aGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCCTGGa                                                          >  1:2335903/1‑101 (MQ=255)
                                             cccgcccgtcgCTCCCCCCGCGGGGGCGTGGATCGAAA.CTGCGTGCGCGCCATG.AATTCATCGAGCGCCATCGAGCGTCGCTCCCCCCGCGGCGGCGTggc                                                          >  2:2853126/7‑100 (MQ=0)
                                                tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                      <  1:1994825/101‑1 (MQ=255)
                                                tGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                      <  1:1345908/101‑1 (MQ=255)
                                                 gtcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                     <  1:3329074/101‑1 (MQ=255)
                                                 gtcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                     >  1:3900193/1‑101 (MQ=255)
                                                   cgtcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                   <  2:14730/101‑1 (MQ=255)
                                                     tcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                 >  1:3410438/1‑101 (MQ=255)
                                                     tcgCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTg                                                 >  1:2397009/1‑101 (MQ=255)
                                                       gCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgcg                                               <  1:3508179/101‑1 (MQ=255)
                                                          ccccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccg                                            <  1:3422787/101‑1 (MQ=255)
                                                           cccccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGC.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGcgccgc                                           <  1:722057/101‑1 (MQ=255)
                                                             cccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCt                                         >  1:2128666/1‑101 (MQ=255)
                                                             cccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCt                                         >  2:1256389/1‑101 (MQ=255)
                                                              ccGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTa                                        >  1:928397/1‑101 (MQ=255)
                                                              ccGCGGGGACGTGGCTCGAAACCTGCG.GCGTGATTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTa                                        <  1:3654887/101‑1 (MQ=255)
                                                               cttccgatcttGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAt                                       <  1:3092116/91‑1 (MQ=255)
                                                                  gggggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGc                                    >  1:1475076/1‑101 (MQ=255)
                                                                    gggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCAt                                  >  2:479063/1‑101 (MQ=255)
                                                                     ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                 <  1:1465410/101‑1 (MQ=255)
                                                                     ggCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATc                                 >  1:1457696/1‑101 (MQ=255)
                                                                      gCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATcc                                >  2:3153011/1‑101 (MQ=255)
                                                                        gTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCgg                              >  2:92040/1‑101 (MQ=255)
                                                                         tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTAtagatcggaag                             >  2:3091886/1‑91 (MQ=255)
                                                                         tGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGt                             >  1:510830/1‑101 (MQ=255)
                                                                            aTCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGa                          >  2:2725884/1‑101 (MQ=255)
                                                                               gAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAgcg                       <  1:1603563/101‑1 (MQ=255)
                                                                                 aaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgcc                     >  2:403740/1‑101 (MQ=255)
                                                                                 aaCCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGcgcc                     <  2:2364187/101‑1 (MQ=255)
                                                                                    tccgatctgTGAGTTGCTCGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCAAGcgcctt                   <  2:2539471/93‑3 (MQ=25)
                                                                                     tGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCgcgc                 <  2:1645035/101‑1 (MQ=255)
                                                                                       cg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAg               <  1:1925214/101‑1 (MQ=255)
                                                                                       cg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAg               <  1:2204414/101‑1 (MQ=255)
                                                                                       cg.gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAg               <  1:536646/101‑1 (MQ=255)
                                                                                          gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTc             <  2:3977821/101‑1 (MQ=255)
                                                                                          gcgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGGCGCTCCCCCCGCGGGGGCGTGAATCGAAACTGCGCCGCCTTTCGCATCCGGTCGAGCGCCGCGCAGTc             <  2:3848802/101‑1 (MQ=255)
                                                                                           ggTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg            <  2:2811824/100‑1 (MQ=255)
                                                                                           cgTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCg            >  1:3001718/1‑101 (MQ=255)
                                                                                            gTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCAAGcgccttagatcgga           >  1:2539662/1‑91 (MQ=37)
                                                                                              gAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTc         <  2:2966121/101‑1 (MQ=255)
                                                                                               aGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcc        <  1:896238/101‑1 (MQ=255)
                                                                                                 ttGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTcccc      >  2:560565/1‑101 (MQ=255)
                                                                                                  tGCCAGAGATAGA.CCCCGTAGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTccccc     <  2:640224/101‑1 (MQ=255)
                                                                                                    ccAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcg   <  2:2588576/101‑1 (MQ=255)
                                                                                                    ccAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCcg   <  2:1802237/101‑1 (MQ=255)
                                                                                                     |                                                                                                      
CGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCG.GCGTGAGTTGCCAGTGATAGA.CGCCGTGGA.CGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGCGCCGCCTATCGCATCCGGTCGAGCGCCGCGCAGTCGCTCCCCCCGC  >  NC_014034/1395712‑1395912

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