Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,538 G→C 6.4% intergenic (+350/‑581) 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,5380G→C 6.4% 34.8 ‑1.8 ‑0.0304intergenic (+350/‑581)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/7):  ref base (134/150):  total (147/157)
Fisher's exact test strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.03e-01

Only 100 of 328 total aligned reads displayed.

CTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGC  >  NC_014034/1395438‑1395637
                                                                                                    |                                                                                                     
cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                                                                       <  1:963569/101‑1 (MQ=255)
cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCCTGGATCGAAACATCTCgg                                                                                                       >  2:796291/1‑101 (MQ=255)
 tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                                                                                      >  1:1401946/1‑101 (MQ=255)
     cccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCg                                                                                                  <  2:1716733/101‑1 (MQ=255)
       cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                                                                                >  2:4364560/1‑101 (MQ=255)
        gcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCg                                                                                               <  2:1991481/101‑1 (MQ=255)
         cgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGc                                                                                              >  1:3435219/1‑101 (MQ=255)
          gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                                                             >  2:3386569/1‑101 (MQ=255)
          gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                                                             <  1:2495025/101‑1 (MQ=255)
          gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAAATGTCGGTGACATCggt                                                                                             >  2:1901307/1‑99 (MQ=38)
            gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                                                                                           <  1:564219/101‑1 (MQ=255)
            gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                                                                                           >  2:4653729/1‑101 (MQ=255)
               cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Ccc                                                                                       >  1:1158353/1‑101 (MQ=255)
                gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCTGTGCCGCCGCTCG.TCCa                                                                                      <  2:3877787/101‑1 (MQ=255)
                gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                      <  2:1191325/101‑1 (MQ=255)
                 tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                     >  2:600199/1‑101 (MQ=255)
                 tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                     <  2:2982669/101‑1 (MQ=255)
                  ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGc                                                                                    >  1:4643952/1‑101 (MQ=255)
                      cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.c                                                                               >  1:3863634/1‑101 (MQ=255)
                       gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.GCCAGCCGG.cc                                                                              >  2:2575544/1‑101 (MQ=255)
                         aaCATCGCCGCTTTGGCGGTCGACAGGCTGGGAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTg                                                                            >  1:4133367/1‑101 (MQ=255)
                           cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCt                                                                          <  2:2693202/101‑1 (MQ=255)
                            aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                         <  2:1849123/101‑1 (MQ=255)
                            aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                         <  2:3750842/101‑1 (MQ=255)
                             tCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcg                                                                        <  2:2102500/101‑1 (MQ=255)
                              cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgt                                                                       <  1:3152431/101‑1 (MQ=255)
                              cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgt                                                                       <  1:4364852/101‑1 (MQ=255)
                                ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtcg                                                                     >  1:3599679/1‑101 (MQ=255)
                                 cgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGc                                                                    >  1:424543/1‑101 (MQ=255)
                                    tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcc                                                                 <  2:1550479/101‑1 (MQ=255)
                                       ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTccccc                                                              <  2:1286611/101‑1 (MQ=255)
                                        gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                             >  2:3068333/1‑101 (MQ=255)
                                        gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                             >  2:1489786/1‑101 (MQ=255)
                                        gCGGTCCACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                             <  2:714088/101‑1 (MQ=255)
                                              gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCcggg                                                       >  1:2991581/1‑101 (MQ=255)
                                                    gccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  1:89064/91‑1 (MQ=0)
                                                    cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  2:3562191/101‑1 (MQ=255)
                                                    aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTgg                                                 <  1:26870/91‑1 (MQ=2)
                                                    aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTgg                                                 >  1:1515942/11‑101 (MQ=2)
                                                     tGGTAACGTGTCCCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGa                                                <  2:924568/101‑1 (MQ=255)
                                                      gctaccGTGTCGCCCCCCACGCGGGGGCGGGTGCCGAACCACCTCGGTGCCGACGCTTG.TCCAGCCGG.CCTGCTCGCCGCTCCCCCCGCGGGGGCGTGGAt                                               <  1:1098424/96‑1 (MQ=19)
                                                       gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGGTc                                              >  2:1001599/8‑101 (MQ=2)
                                                       ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                              <  1:3056539/94‑1 (MQ=2)
                                                       ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                              >  2:363785/8‑101 (MQ=2)
                                                         aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                            >  2:363133/1‑101 (MQ=255)
                                                          tgcggtCGCGCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           <  1:2159723/97‑1 (MQ=2)
                                                          ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           >  2:2338999/5‑101 (MQ=2)
                                                          ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           <  2:1627459/97‑1 (MQ=2)
                                                          aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                           <  1:130389/101‑1 (MQ=255)
                                                           gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          <  1:2817406/98‑1 (MQ=2)
                                                           cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          >  1:4471838/4‑101 (MQ=2)
                                                           cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                          <  2:662342/101‑1 (MQ=255)
                                                            gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                         <  1:3751065/101‑1 (MQ=255)
                                                            gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGCCGCTCG.TCCAGCCGG.CCTGCTCGCCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                         >  2:2171356/1‑101 (MQ=255)
                                                            acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                         >  1:4591596/3‑101 (MQ=2)
                                                             tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                        <  2:4581392/101‑1 (MQ=255)
                                                              gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                       <  2:3430560/101‑2 (MQ=2)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccac                                      <  2:413716/101‑1 (MQ=255)
                                                               tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGCGGGCGTGGATCGAAacccc                                      >  1:1932749/1‑99 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                     <  1:20762/101‑1 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                     >  1:3234017/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGgaaatgcg                                     >  1:511801/1‑96 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcg                                     >  2:2811199/1‑97 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                     <  2:1332531/101‑4 (MQ=2)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                     >  1:2272647/1‑97 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                    >  1:3128882/1‑101 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                    <  1:2237231/101‑1 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                    <  1:760688/101‑6 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgc                                    >  1:851367/1‑96 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCtt                                   <  2:1789683/101‑1 (MQ=255)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtg                                   <  1:365563/101‑7 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcg                                   <  1:4672490/101‑6 (MQ=2)
                                                                  cTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcg                                   >  1:3462210/1‑96 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                  >  2:4549592/1‑94 (MQ=2)
                                                                   cccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                  <  1:2901129/100‑8 (MQ=2)
                                                                    ccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCATGGATCGAAACCACCTTGc                                 <  1:1540702/101‑1 (MQ=255)
                                                                    ccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcg                                 >  1:651747/1‑93 (MQ=2)
                                                                     cccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                                <  2:4284612/101‑1 (MQ=255)
                                                                      ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgggt                               >  2:3311105/1‑91 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                               >  1:2418490/1‑91 (MQ=2)
                                                                      ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                               <  1:1625567/101‑11 (MQ=2)
                                                                        ccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGGCGCTCG.CCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCTTGCGGAt                             <  1:167807/101‑1 (MQ=255)
                                                                             ggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCa                        <  2:4089418/101‑1 (MQ=255)
                                                                              gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                       <  1:2538023/101‑1 (MQ=255)
                                                                              gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                       >  2:2450173/1‑101 (MQ=255)
                                                                              gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                       >  2:4613261/1‑101 (MQ=255)
                                                                               ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                      <  2:822148/101‑1 (MQ=255)
                                                                                 cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                    <  2:3233772/101‑1 (MQ=255)
                                                                                  gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                   <  2:4423923/101‑1 (MQ=255)
                                                                                  gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                   >  2:3158274/1‑101 (MQ=255)
                                                                                         gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc            >  1:4621023/1‑101 (MQ=255)
                                                                                         gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCCGGGc            >  2:2870648/1‑101 (MQ=255)
                                                                                          aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg           <  1:1828130/101‑1 (MQ=255)
                                                                                          aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg           <  2:4232/101‑1 (MQ=255)
                                                                                          aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg           <  1:3368984/101‑1 (MQ=255)
                                                                                              aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAgg       >  2:2270624/1‑100 (MQ=255)
                                                                                              aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGc       <  2:2350607/101‑1 (MQ=255)
                                                                                              aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGc       >  2:1052279/1‑101 (MQ=255)
                                                                                                 tcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTc    <  1:3524789/101‑1 (MQ=255)
                                                                                                   ggTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCAGATTGTCACTTTCGCGGGGCGAAGCGTCGc  >  1:4232931/1‑101 (MQ=255)
                                                                                                    |                                                                                                     
CTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGC  >  NC_014034/1395438‑1395637

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