Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,535 T→G 16.4% intergenic (+347/‑584) 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,5350T→G16.4% 145.8 ‑2.0 ‑0.2410intergenic (+347/‑584)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (33/22):  ref base (135/145):  total (207/203)
Fisher's exact test strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.94e-01

Only 100 of 478 total aligned reads displayed.

GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395434‑1395635
                                                                                                     |                                                                                                        
 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                                                                          <  2:5773985/101‑1 (MQ=255)
 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                                                                          >  2:1560552/1‑101 (MQ=255)
  cGCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgcttt                                                                                                        >  2:213881/1‑99 (MQ=2)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                                                                                       >  2:4703055/1‑98 (MQ=2)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                        <  2:5235566/101‑1 (MQ=255)
   gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTGGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                        <  1:774801/101‑1 (MQ=255)
    cTCCCCCCGGGGGGGCGTGGATCGAAACATCGCCGCTTTGGGGGGCGACAGGCTGGTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACAtctccg                                                                                                       <  2:4125728/101‑3 (MQ=37)
    cTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttca                                                                                                      >  1:985534/1‑97 (MQ=2)
    cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg.g                                                                                                      >  2:228094/1‑101 (MQ=255)
     tCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaa                                                                                                     >  1:3481832/1‑96 (MQ=2)
     tCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaa                                                                                                     >  1:5565866/1‑96 (MQ=2)
     tCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaa                                                                                                     >  2:2439313/1‑96 (MQ=2)
      ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                                                                    >  2:4824836/1‑95 (MQ=0)
              gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCggt                                                                                            >  2:3572195/1‑99 (MQ=255)
               gcggCGTGGATCGAACCATCGCCACTATGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTc                                                                                           <  2:2649431/99‑1 (MQ=255)
                 ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.t                                                                                        <  2:1772832/101‑1 (MQ=255)
                   cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGTATGCCGTcttctgct                                                                                        >  1:5293820/1‑94 (MQ=37)
                    gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCa                                                                                     <  2:3027110/101‑1 (MQ=255)
                     tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAg                                                                                    <  1:2684305/101‑1 (MQ=255)
                      ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGc                                                                                   <  2:4554082/101‑1 (MQ=255)
                       gATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGcc                                                                                  >  1:3367475/1‑101 (MQ=255)
                         tCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCgg                                                                                >  1:440305/1‑101 (MQ=255)
                           gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.cc                                                                             >  1:1663205/1‑101 (MQ=255)
                            aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            >  1:2930073/1‑101 (MQ=255)
                            aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            <  2:4106099/101‑1 (MQ=255)
                            aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            >  1:5465414/1‑101 (MQ=255)
                            aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            >  1:3841906/1‑101 (MQ=255)
                             aactgcgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCGCCCCCCGCGGGGGCGTGGATCGCAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                          <  2:1277447/96‑1 (MQ=255)
                             aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                          <  1:3062437/101‑1 (MQ=255)
                             aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                          <  1:1326460/101‑1 (MQ=255)
                               cATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCt                                                                        >  1:1446338/1‑101 (MQ=255)
                                aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                       >  1:2438454/1‑101 (MQ=255)
                                       cTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGGCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTc                                                                >  1:3637435/1‑101 (MQ=255)
                                          tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccc                                                             <  1:423811/101‑1 (MQ=255)
                                              ggTCGACAGGCTGGTAACGTGGCGCTCCCCCCGCGGGGGCGTGGATCGAGACCTCTAG.GGGTTGTCGCTCG.CCCCGCCGG.CCT.GCTCCTCCCTcccccccc                                                         >  2:1983678/1‑99 (MQ=21)
                                                tCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGTCTCG.GAGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCCCgg                                                       >  2:2855760/1‑101 (MQ=255)
                                                 cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggg                                                      >  1:353716/1‑101 (MQ=255)
                                                  gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                     <  1:5153307/101‑1 (MQ=255)
                                                   aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                    >  2:376208/1‑101 (MQ=255)
                                                     aGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCg                                                  >  2:5664604/1‑101 (MQ=255)
                                                      ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGt                                                 <  1:4003006/101‑1 (MQ=255)
                                                       gCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTg                                                <  1:2206726/101‑1 (MQ=255)
                                                        cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCCCGGGGGCGTgg                                               >  2:3437753/1‑101 (MQ=255)
                                                        cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTAGCTCG.TCCAGCCGG.CCT.GATCGTCGCTCCCCCCCCGGGGGCGTgg                                               >  1:4866181/1‑101 (MQ=255)
                                                          ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                             >  2:1060318/1‑101 (MQ=255)
                                                           gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGacc                                            <  1:2125864/94‑3 (MQ=0)
                                                           gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                            <  2:2932105/94‑1 (MQ=2)
                                                           gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                            <  2:187755/94‑1 (MQ=2)
                                                           gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                            <  2:464786/101‑1 (MQ=255)
                                                            tgccgcgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGAT.GATCAAATCCGCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                            <  2:827328/95‑1 (MQ=2)
                                                            tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           >  1:611848/1‑101 (MQ=255)
                                                            tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           <  2:2337974/101‑1 (MQ=255)
                                                            tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.ACT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           <  1:236097/101‑1 (MQ=255)
                                                              aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  1:690549/101‑1 (MQ=255)
                                                               cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                        >  2:5455354/4‑101 (MQ=2)
                                                               cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                        <  2:5019911/101‑1 (MQ=255)
                                                                acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                       >  1:3215625/3‑101 (MQ=2)
                                                                acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAac                                       >  2:4506727/3‑101 (MQ=2)
                                                                 tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      >  2:4372877/1‑101 (MQ=255)
                                                                 ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      <  2:149606/100‑1 (MQ=2)
                                                                 cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                      >  2:1414322/2‑100 (MQ=2)
                                                                  gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                     <  1:156012/101‑1 (MQ=255)
                                                                  gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacca                                     >  1:847506/1‑101 (MQ=255)
                                                                  gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacta                                     <  2:4742461/101‑3 (MQ=2)
                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                    <  2:3970681/101‑3 (MQ=2)
                                                                   tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                    >  1:3457500/1‑99 (MQ=2)
                                                                     tcgcccccaGCGGGGGCGTGGATTGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                  <  2:5098867/98‑6 (MQ=0)
                                                                     gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                  <  1:471534/101‑6 (MQ=2)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtg                                 >  1:5062185/1‑95 (MQ=2)
                                                                      cTCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtg                                 >  1:2856248/1‑95 (MQ=2)
                                                                       tCCCCCCGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                <  1:4599151/101‑8 (MQ=2)
                                                                        ccccccGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCttcc                               >  2:4918852/1‑99 (MQ=255)
                                                                        ccccccGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                               >  1:257802/1‑101 (MQ=255)
                                                                        ccccccGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                               <  2:5563063/101‑8 (MQ=2)
                                                                        ccccccGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                               >  1:110422/1‑93 (MQ=2)
                                                                         cccccGCGGGGGCGTGGATCGAAACAAAAAG.CCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcgccg                              >  2:766770/1‑93 (MQ=0)
                                                                         cccccGCGCGGGCTTGGGTCGAAGCGTCTCG.ATGCCGTCGCCCT.TCCAGCCGG.CCT.GCTGGTCGCGCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                              <  2:5293530/101‑1 (MQ=25)
                                                                          cccgGCGGGGGCGTGGATCGAAACAATGCG.GCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgcg                             <  2:1691699/101‑11 (MQ=0)
                                                                          ccccGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCgg                             <  1:1291006/101‑1 (MQ=255)
                                                                          ccccGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                             >  1:3744870/1‑91 (MQ=2)
                                                                          ccccGCGGGGGCGTGGATCGAAACA.CCCG.GCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                             >  2:1168536/1‑91 (MQ=2)
                                                                             cgcgGGGGCGTGGATCGAAGCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAtt                          <  1:2288666/101‑1 (MQ=255)
                                                                              gcgGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                         <  2:637539/101‑1 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  1:4679000/1‑101 (MQ=255)
                                                                                gggggCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  2:2557926/1‑101 (MQ=255)
                                                                                    gCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                   <  1:1006253/101‑1 (MQ=255)
                                                                                     cGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                  <  1:1939080/101‑1 (MQ=255)
                                                                                      gTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                 <  2:2047717/101‑1 (MQ=255)
                                                                                       tGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                >  1:2165979/1‑101 (MQ=255)
                                                                                         gATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg              >  1:4135716/1‑101 (MQ=255)
                                                                                         gATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg              >  1:2467323/1‑101 (MQ=255)
                                                                                             gAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc          <  1:5008258/101‑1 (MQ=255)
                                                                                             caaaCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc          <  2:321067/100‑1 (MQ=255)
                                                                                              aaaCATTTCG.CTGCCCTGGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCCCTTGCGGATTGTCACTTTCGCGGGGCg         <  2:2792827/101‑1 (MQ=37)
                                                                                              aaaCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg         <  2:431107/101‑1 (MQ=255)
                                                                                               aaCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa        <  1:2738048/101‑1 (MQ=255)
                                                                                                aCATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGaa       <  2:92620/101‑1 (MQ=255)
                                                                                                    ctcG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGt   <  1:1060393/101‑1 (MQ=255)
                                                                                                     tcG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTc  <  1:5676786/101‑1 (MQ=255)
                                                                                                     tcG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTc  >  1:1846979/1‑101 (MQ=255)
                                                                                                     |                                                                                                        
GTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCG.GTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395434‑1395635

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 39 ≤  ATCG