Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,535 T→C 23.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→C23.4% 128.6 ‑1.9 ‑1.9302intergenic (+347/‑584)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (29/20):  ref base (89/101):  total (147/155)
Fisher's exact test strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.31e-02

Only 100 of 333 total aligned reads displayed.

TCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395435‑1395635
                                                                                                    |                                                                                                       
tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtct                                                                                                         <  2:1977338/101‑1 (MQ=255)
 cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                        <  1:2853850/101‑1 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttc                                                                                                      >  1:1979495/1‑98 (MQ=2)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                       <  1:459324/101‑1 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                       >  1:831056/1‑101 (MQ=255)
  gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                       >  1:3511728/1‑101 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGt                                                                                                     >  1:2359869/1‑101 (MQ=255)
     ccccccGCGGGGGTGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTg                                                                                                    <  2:1047259/101‑1 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                                                                   <  1:2748894/101‑7 (MQ=0)
     ccccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaat                                                                                                   <  1:1767511/101‑7 (MQ=0)
     ccccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGTCCGACAGGCCGTTATAGGGGCGGTCGCCCCGCGGGGCCGGGGTTCGAAACTACTTGGTg                                                                                                    >  2:833080/1‑101 (MQ=16)
      cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGc                                                                                                   <  1:2472404/101‑1 (MQ=255)
      cccccGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGGGc                                                                                                   >  2:2952045/1‑101 (MQ=255)
       ccccGCGGGGGCGTGGATCGAAACTACACCGGGGCTCTGGTC.ACGGTTTGGGACGGAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAgctttcaatgg                                                                                                 <  1:3321997/101‑9 (MQ=2)
          cgcgGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTc                                                                                               <  1:3219436/101‑1 (MQ=255)
             gggggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGcc                                                                                            >  2:1988290/1‑100 (MQ=255)
                 gCGTTGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCTCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.Tc                                                                                       <  2:554174/101‑1 (MQ=255)
                  cGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTTGATCGAAACATCTCGGTGCCGTCGCTCG.Tcc                                                                                      <  1:3754861/101‑1 (MQ=255)
                   gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                     >  1:1936922/1‑101 (MQ=255)
                   gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                     <  2:2359693/101‑1 (MQ=255)
                    tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                    <  1:1341215/101‑1 (MQ=255)
                       aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                 >  2:3402862/1‑101 (MQ=255)
                       aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                 >  2:3624506/1‑101 (MQ=255)
                          gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                             >  2:346787/1‑101 (MQ=255)
                          gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                             >  2:3085897/1‑101 (MQ=255)
                           aaaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCt                                                                            <  2:587678/101‑1 (MQ=255)
                            aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.g                                                                          <  1:2257267/101‑1 (MQ=255)
                             aCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.Gc                                                                         >  2:630643/1‑101 (MQ=255)
                               aTCGCCGCTTTGGCGGTCGACAGCCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCtc                                                                       <  2:544165/101‑1 (MQ=255)
                                       tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcc                                                               >  1:1740061/1‑101 (MQ=255)
                                       tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcc                                                               >  2:1557093/1‑101 (MQ=255)
                                         tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccc                                                             >  1:2098719/1‑101 (MQ=255)
                                         tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCACCCGG.CCT.GCTCGTCGCTcccc                                                             <  1:2411508/101‑1 (MQ=255)
                                          ggCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTccccc                                                            >  1:4171262/1‑101 (MQ=255)
                                           gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTcccccc                                                           <  2:1162012/101‑1 (MQ=255)
                                             ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCcgc                                                         >  1:2643189/1‑101 (MQ=255)
                                                 gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                     <  2:2151676/101‑1 (MQ=255)
                                                 gACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCgggg                                                     <  1:2018102/101‑1 (MQ=255)
                                                  aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCggggg                                                    >  2:2081737/1‑101 (MQ=255)
                                                   cAGGCTGGTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGc                                                   <  1:579254/101‑1 (MQ=255)
                                                       ttgcgcgctggtCGCTCCCCCCGCGGGGGCGTGGATCGAAAC.TTGATGACCAGATCCTCA.ATCGGCCGG.GCTCGATCGTCGCTCCCCCCGCGGGGGCGTgg                                               <  2:1605716/91‑1 (MQ=0)
                                                       agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTgg                                               <  2:2285179/91‑1 (MQ=0)
                                                       agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTgg                                               >  2:1021061/11‑101 (MQ=0)
                                                       aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTgg                                               >  1:641202/11‑101 (MQ=2)
                                                       aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTgg                                               >  2:459927/11‑101 (MQ=2)
                                                       aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTgg                                               >  1:1035369/11‑101 (MQ=2)
                                                        ccctcgcccgtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGa                                              <  1:3575894/92‑1 (MQ=2)
                                                        ccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGa                                              <  1:1796616/92‑1 (MQ=2)
                                                        aggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGa                                              >  1:3435224/10‑101 (MQ=0)
                                                         ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                             <  1:1837026/101‑1 (MQ=255)
                                                         cctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGAt                                             <  2:404372/93‑1 (MQ=2)
                                                           tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                           >  1:440169/1‑101 (MQ=255)
                                                           gctgcggtCGCCCCCCCCGCGGGGGCGTGGATCAAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGCGCGTGGATCg                                           <  1:3658676/95‑1 (MQ=0)
                                                            cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          <  1:3729463/96‑1 (MQ=2)
                                                            cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          <  2:1713020/96‑1 (MQ=2)
                                                            aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                          >  1:1096869/1‑101 (MQ=255)
                                                             ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  1:3373446/97‑1 (MQ=2)
                                                             gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGGGGGGGTCgga                                         >  1:2962384/5‑99 (MQ=2)
                                                             cgttgcccccccccccGGGGGGGGGGGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  1:3937285/91‑1 (MQ=37)
                                                             aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                         <  1:2511906/101‑1 (MQ=255)
                                                              gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                        >  2:3552237/4‑101 (MQ=2)
                                                               tggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                       <  1:1526138/99‑1 (MQ=2)
                                                               gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                       <  1:3001814/101‑1 (MQ=255)
                                                               acgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCC.GCCCGTCGCTCCCCCCGCGGGGGCGTGGCTCGAAac                                       >  2:721106/3‑101 (MQ=2)
                                                                ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      <  2:961627/100‑1 (MQ=2)
                                                                ggtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAAAACCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                      <  1:1378143/100‑1 (MQ=2)
                                                                 gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacta                                     <  2:3126514/101‑3 (MQ=2)
                                                                 ggcGCGCCCCCCGGGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                     <  2:1806895/99‑2 (MQ=0)
                                                                  tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccac                                    <  1:962466/101‑1 (MQ=255)
                                                                   cGCGCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                   <  2:2554911/101‑5 (MQ=2)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCt                                  >  2:3322624/1‑101 (MQ=255)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                  >  2:899049/1‑97 (MQ=2)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCtt                                 >  1:3155118/1‑101 (MQ=255)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCtt                                 <  2:4126783/101‑1 (MQ=255)
                                                                     cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgcg                                 >  2:3965954/1‑95 (MQ=2)
                                                                       ccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                               <  2:1096775/101‑1 (MQ=255)
                                                                       ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                               >  1:1261818/1‑94 (MQ=2)
                                                                       ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCC.G.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                               >  2:3931076/1‑94 (MQ=2)
                                                                        cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                              <  1:277709/101‑10 (MQ=2)
                                                                        cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                              <  2:2618532/101‑10 (MQ=2)
                                                                         ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACA.GGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                             <  2:2518090/101‑11 (MQ=2)
                                                                          cccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCACCTTGCggg                            >  2:276987/1‑100 (MQ=255)
                                                                             gcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                         <  1:3892067/101‑1 (MQ=255)
                                                                              cgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                        <  2:1447839/101‑1 (MQ=255)
                                                                               gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTggc                       >  2:1908889/1‑99 (MQ=255)
                                                                               gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  2:1310057/1‑101 (MQ=255)
                                                                               gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                       >  2:2415076/1‑101 (MQ=255)
                                                                                ggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCa                      >  1:3309431/1‑101 (MQ=255)
                                                                                     gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                 <  1:117110/101‑1 (MQ=255)
                                                                                      tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                >  2:2885313/1‑101 (MQ=255)
                                                                                      tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                <  2:4171000/101‑1 (MQ=255)
                                                                                      tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                <  1:3086127/101‑1 (MQ=255)
                                                                                       ggATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcag               >  1:3017568/1‑99 (MQ=255)
                                                                                        gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg              >  1:4031179/1‑101 (MQ=255)
                                                                                          tCGAAGCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCggg            <  2:388774/101‑1 (MQ=255)
                                                                                             aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcgagatc         >  2:177818/1‑96 (MQ=255)
                                                                                              aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa        >  2:3812846/1‑101 (MQ=255)
                                                                                              aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa        <  2:2478036/101‑1 (MQ=255)
                                                                                                cATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAg      >  2:3193785/1‑101 (MQ=255)
                                                                                                  tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCg    >  1:1151060/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
TCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCT.GCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTC  >  NC_014034/1395435‑1395635

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