Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,534 C→T 10.7% intergenic (+346/‑585) 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,5340C→T10.7% 124.7 ‑0.4 0.0509intergenic (+346/‑585)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (28/25):  ref base (210/214):  total (253/256)
Fisher's exact test strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 576 total aligned reads displayed.

GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGT  >  NC_014034/1395434‑1395634
                                                                                                    |                                                                                                      
gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtc                                                                                                        >  1:3642701/1‑101 (MQ=255)
  cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                      >  2:19991/1‑101 (MQ=255)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                                                                                    <  1:3516848/101‑4 (MQ=2)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                     <  2:6248931/101‑1 (MQ=255)
    cgccccgcGCGGGGGCGTGGATCGCACCATCGCCGCTTTGGCCGTCGACAGGCTCGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                                                                    <  2:65015/99‑1 (MQ=37)
    cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                                                                   >  1:6738877/1‑97 (MQ=2)
       cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGc                                                                                                 <  1:3370757/101‑1 (MQ=255)
          ccGCGGGGGCGTGGATCGAATCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGt                                                                                              <  1:5604294/101‑1 (MQ=255)
          ccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGt                                                                                              >  1:1891619/1‑101 (MQ=255)
              gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCt                                                                                          <  1:1727795/101‑1 (MQ=255)
               ggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTc                                                                                         >  2:1677370/1‑101 (MQ=255)
                gggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCg                                                                                        >  1:3928541/1‑101 (MQ=255)
                    gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                   <  2:219770/101‑1 (MQ=255)
                     tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                  >  2:4308110/1‑101 (MQ=255)
                      ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGc                                                                                 >  2:3809545/1‑101 (MQ=255)
                          cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGCAACATCTCGGTGCCGTCGCTCG.TCCAGCCTG.c                                                                            >  2:7067430/1‑101 (MQ=255)
                           gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                           <  1:5350350/101‑1 (MQ=255)
                           gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                           >  1:4575201/1‑101 (MQ=255)
                             aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTg                                                                         <  1:2240420/101‑1 (MQ=255)
                             aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTg                                                                         >  1:6098403/1‑101 (MQ=255)
                              aCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGc                                                                        <  2:4771977/101‑1 (MQ=255)
                                aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                      >  2:192184/1‑101 (MQ=255)
                                aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                      >  1:2880381/1‑101 (MQ=255)
                                 tCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcg                                                                     <  1:644849/101‑1 (MQ=255)
                                   tctgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtc                                                                   <  1:6156252/98‑1 (MQ=255)
                                      gcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGcc                                                                >  2:3440429/1‑100 (MQ=255)
                                      gcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCt                                                                <  2:6306192/101‑1 (MQ=255)
                                       cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTc                                                               >  1:6594888/1‑101 (MQ=255)
                                       cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTc                                                               >  1:2404489/1‑101 (MQ=255)
                                        tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcc                                                              >  2:3252295/1‑101 (MQ=255)
                                           ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTccccc                                                           >  2:2457311/1‑101 (MQ=255)
                                            gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                          <  2:5440726/101‑1 (MQ=255)
                                              ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcgc                                                        <  1:1780708/101‑1 (MQ=255)
                                               gTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TGCCGCCGG.CCTGCTCGTCGCTCCCCCcccg                                                       >  2:6163043/1‑101 (MQ=255)
                                                 cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCggg                                                     >  2:1822824/1‑101 (MQ=255)
                                                     aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCg                                                 >  1:5749150/1‑101 (MQ=255)
                                                     aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCg                                                 >  2:4053940/1‑101 (MQ=255)
                                                      ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGt                                                <  1:4187910/101‑1 (MQ=255)
                                                       gCTGGCAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGCTGCCGTCGCTCG.TCCAGCCGG.CCCGCGCGTCGCCACCCCCCCGGCGGCCTg                                               >  2:6723102/1‑101 (MQ=25)
                                                        agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTgg                                              <  2:2388445/91‑1 (MQ=0)
                                                        aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTgg                                              <  2:1053672/91‑1 (MQ=2)
                                                          ggccggtggtCGCGCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  2:2786868/93‑1 (MQ=2)
                                                          ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  1:4069645/101‑1 (MQ=255)
                                                          cctcgcccgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  2:1453383/93‑1 (MQ=2)
                                                          cctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                            <  1:1042213/93‑1 (MQ=2)
                                                           gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  2:3367885/94‑1 (MQ=2)
                                                           ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  1:2614646/94‑1 (MQ=2)
                                                           ctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATc                                           <  1:1563365/94‑1 (MQ=2)
                                                            tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                          <  1:3809706/101‑1 (MQ=255)
                                                             cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                         >  2:576294/6‑101 (MQ=2)
                                                              gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  1:6974482/97‑1 (MQ=2)
                                                              gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        >  2:2930367/5‑101 (MQ=2)
                                                              gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  1:5190702/97‑1 (MQ=2)
                                                              aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                        <  1:3450370/101‑1 (MQ=255)
                                                                gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      <  1:4346602/101‑1 (MQ=255)
                                                                gtgtCCCGCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      <  2:2703251/101‑1 (MQ=255)
                                                                ccgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      <  2:3045167/99‑1 (MQ=2)
                                                                ccgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      <  1:6267312/99‑1 (MQ=2)
                                                                acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      <  1:6010375/99‑1 (MQ=2)
                                                                acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                      >  2:5278950/3‑101 (MQ=2)
                                                                 tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                     <  1:5678268/101‑1 (MQ=255)
                                                                  gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAaccc                                    >  2:6615877/1‑100 (MQ=255)
                                                                  gtCGCTCCCCCCGCGGGGGCGTGGAGCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                    <  2:886510/101‑1 (MQ=255)
                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccac                                   <  1:5571122/101‑1 (MQ=255)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                  >  1:2771360/1‑101 (MQ=255)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                  <  2:6071020/101‑4 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                  >  2:1570366/1‑97 (MQ=2)
                                                                    cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                  <  1:4870860/101‑5 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATccgaactgagt                                 >  2:5012884/1‑91 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgg                                 >  1:1432997/1‑96 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                 >  1:1129728/1‑97 (MQ=2)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgc                                 >  2:4196589/1‑96 (MQ=2)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCtt                                <  2:749328/101‑1 (MQ=255)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcg                                >  2:5370035/1‑96 (MQ=2)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TTCAGGCGG.CCCGCTCGTCGATCCCCCCGCGGGGGCGGGGATCGAAACCACCTTg                               >  2:4585460/1‑101 (MQ=37)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                               >  1:233886/1‑94 (MQ=2)
                                                                        ccccccGCGGGGTCGTCGGCCAAAACATCTCGGTGCCGTCGCTCG.TCCAGTCGG.CCCGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                              <  1:7067804/101‑1 (MQ=25)
                                                                        ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                              >  2:5162555/1‑93 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                             <  1:3011482/101‑10 (MQ=2)
                                                                          ccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCgg                            >  1:6331641/1‑101 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCgg                            <  2:1745152/101‑1 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                            <  1:1944764/101‑11 (MQ=2)
                                                                             cgcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAtt                         <  1:1224457/101‑1 (MQ=255)
                                                                               cgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                       <  2:5927532/101‑1 (MQ=255)
                                                                               cgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                       >  1:2548852/1‑101 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                      <  2:2064902/101‑1 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                      <  2:483657/101‑1 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                      <  1:3155529/101‑1 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                      <  1:721727/101‑1 (MQ=255)
                                                                                    gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                  <  1:2073807/101‑1 (MQ=255)
                                                                                      gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                >  2:6921959/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                >  2:1125884/1‑101 (MQ=255)
                                                                                      gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                <  1:1606054/101‑1 (MQ=255)
                                                                                         gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg             >  1:465423/1‑101 (MQ=255)
                                                                                         gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg             >  1:854342/1‑101 (MQ=255)
                                                                                           tCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCggg           >  1:4445581/1‑101 (MQ=255)
                                                                                             gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc         >  1:916201/1‑101 (MQ=255)
                                                                                               aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa       <  2:5636286/101‑1 (MQ=255)
                                                                                                 cATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAg     >  2:4464570/1‑101 (MQ=255)
                                                                                                   tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCg   <  2:2629236/101‑1 (MQ=255)
                                                                                                    |                                                                                                      
GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGT  >  NC_014034/1395434‑1395634

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