Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,549 C→A 26.6% intergenic (+361/‑570) 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,5490C→A26.6% 233.8 ‑1.8 ‑0.0308intergenic (+361/‑570)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (45/37):  ref base (103/123):  total (148/160)
Fisher's exact test strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.72e-01

Only 100 of 322 total aligned reads displayed.

GGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGG  >  NC_014034/1395449‑1395649
                                                                                                    |                                                                                                      
      tGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAg                                                                                                 <  2:2590120/101‑1 (MQ=255)
       ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGc                                                                                                >  2:1836906/1‑101 (MQ=255)
         aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                              <  2:1655313/101‑1 (MQ=255)
         aTCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCg                                                                                              >  2:3624506/1‑101 (MQ=255)
            gAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                                          >  2:346787/1‑101 (MQ=255)
            atctcATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.cc                                                                                          <  2:912609/97‑1 (MQ=255)
              aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTg                                                                                        <  1:2257267/101‑1 (MQ=255)
              aaCATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTg                                                                                        <  2:2167831/101‑1 (MQ=255)
                 aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtc                                                                                     <  2:1567323/101‑1 (MQ=255)
                    gccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtc                                                                                  <  2:2066604/101‑1 (MQ=255)
                          ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTccc                                                                            <  2:4207246/101‑1 (MQ=255)
                             gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                                         <  2:1162012/101‑1 (MQ=255)
                              cGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcg                                                                        >  1:3543589/1‑101 (MQ=255)
                               ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcgc                                                                       >  2:4149808/1‑101 (MQ=255)
                               ggTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcgc                                                                       >  1:2643189/1‑101 (MQ=255)
                                gTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCcccg                                                                      <  1:2416821/101‑1 (MQ=255)
                                  cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGGTCGTCGCTTCCCCCGCggg                                                                    >  2:3205056/1‑101 (MQ=255)
                                    aCAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCggggg                                                                  <  2:814997/101‑1 (MQ=255)
                                     cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGc                                                                 >  1:1794387/1‑101 (MQ=255)
                                     cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGc                                                                 <  2:2623965/101‑1 (MQ=255)
                                     cAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGc                                                                 >  2:1083280/1‑101 (MQ=255)
                                       ggCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGt                                                               <  1:3403107/101‑1 (MQ=255)
                                         gccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTgg                                                             <  1:599133/91‑1 (MQ=0)
                                         agcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACGAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTgg                                                             >  2:3093939/11‑101 (MQ=2)
                                         agcggctgcggtCCCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTgg                                                             <  2:574877/91‑1 (MQ=0)
                                         aaggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTgg                                                             >  1:1035369/11‑101 (MQ=2)
                                          ccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGa                                                            <  2:212854/92‑1 (MQ=2)
                                           ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  2:1984778/93‑1 (MQ=2)
                                           ggTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                           <  1:1837026/101‑1 (MQ=255)
                                            gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  1:1592516/94‑1 (MQ=2)
                                            gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  2:549361/101‑1 (MQ=255)
                                            gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                          <  2:3426737/101‑1 (MQ=255)
                                             tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                         >  2:508899/1‑101 (MQ=255)
                                              cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  1:3729463/96‑1 (MQ=2)
                                              cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        <  2:1713020/96‑1 (MQ=2)
                                              cgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  1:1951001/6‑101 (MQ=2)
                                              aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                        >  1:1096869/1‑101 (MQ=255)
                                               ggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  2:3191313/97‑1 (MQ=2)
                                               cgttgcccccccccccGGGGGGGGGGGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                       <  1:3937285/91‑1 (MQ=37)
                                                gtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  2:3552237/4‑101 (MQ=2)
                                                cacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                      >  2:2396321/4‑101 (MQ=2)
                                                 tggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                                     >  1:3003947/3‑101 (MQ=2)
                                                  ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                                    <  2:961627/100‑1 (MQ=2)
                                                  ggtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAAAACCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                                    <  1:1378143/100‑1 (MQ=2)
                                                  cgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                                    <  1:444287/100‑2 (MQ=2)
                                                   gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                                   >  1:418309/1‑99 (MQ=2)
                                                   ggcGCGCCCCCCGGGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccc                                                   <  2:1806895/99‑2 (MQ=0)
                                                    tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                                  <  2:58969/101‑4 (MQ=2)
                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                                 <  1:329804/101‑1 (MQ=255)
                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                                 <  2:1341585/101‑4 (MQ=2)
                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccg                                                 >  1:543445/1‑98 (MQ=2)
                                                     cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                                 >  1:3801984/1‑97 (MQ=2)
                                                      gCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgg                                                >  2:3956885/1‑96 (MQ=2)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCtt                                               >  1:3155118/1‑101 (MQ=255)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCCCGGGGGCGGGGGTCGAAaccgcgtg                                               >  1:764575/1‑96 (MQ=0)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgct                                               <  2:3393013/101‑7 (MQ=2)
                                                       cTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgcg                                               >  2:1362953/1‑95 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                              >  1:2974060/1‑101 (MQ=255)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                              >  2:664447/1‑101 (MQ=255)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              <  1:1225729/101‑8 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              >  1:3319607/1‑94 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              >  1:2133854/1‑94 (MQ=2)
                                                        tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              <  2:1035278/101‑8 (MQ=2)
                                                        cccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                              <  1:3625092/100‑8 (MQ=2)
                                                         ccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctgg                                             >  1:896318/1‑93 (MQ=2)
                                                          cccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggc                                            <  1:277709/101‑10 (MQ=2)
                                                           ccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgggcgcg                                           >  2:4119558/1‑91 (MQ=2)
                                                           ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                                           <  2:2518090/101‑11 (MQ=2)
                                                             ccGCGGGGGCGTGGATCGAAACATCTCGGCGCCGTCGCTCG.TCCCGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGAt                                         <  1:2483203/101‑1 (MQ=255)
                                                               gcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                                       <  1:3892067/101‑1 (MQ=255)
                                                               gcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                                       >  2:2951296/1‑101 (MQ=255)
                                                                 gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTGGCGGATTGTc                                     >  1:965997/1‑101 (MQ=255)
                                                                   gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                                   <  1:3033485/101‑1 (MQ=255)
                                                                   gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                                   <  1:3678978/101‑1 (MQ=255)
                                                                   gggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                                   <  1:2726450/101‑1 (MQ=255)
                                                                    ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACt                                  <  2:886684/101‑1 (MQ=255)
                                                                     gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                                 >  2:3068781/1‑101 (MQ=255)
                                                                      cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                >  1:1354653/1‑101 (MQ=255)
                                                                      cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                <  2:391545/101‑1 (MQ=255)
                                                                      cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                                <  1:1988416/101‑1 (MQ=255)
                                                                       gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCGCTTTc                               >  2:2778443/1‑101 (MQ=255)
                                                                       gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                               <  1:117110/101‑1 (MQ=255)
                                                                       gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                               >  2:4001114/1‑101 (MQ=255)
                                                                        tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                              <  1:3086127/101‑1 (MQ=255)
                                                                          gATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcg                            >  1:4031179/1‑101 (MQ=255)
                                                                             cGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgggg                         >  1:174884/1‑101 (MQ=255)
                                                                              gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc                        >  2:3873211/1‑101 (MQ=255)
                                                                               aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcgcgagatc                       >  2:177818/1‑96 (MQ=255)
                                                                               aaaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCg                       >  2:2330706/1‑101 (MQ=255)
                                                                                aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa                      >  2:3812846/1‑101 (MQ=255)
                                                                                aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa                      <  2:2478036/101‑1 (MQ=255)
                                                                                  cATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAg                    >  1:2890949/1‑101 (MQ=255)
                                                                                    tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCg                  >  1:1151060/1‑101 (MQ=255)
                                                                                      tcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTc                <  2:3309190/101‑1 (MQ=255)
                                                                                        ggTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGc              <  2:2418301/101‑1 (MQ=255)
                                                                                          tGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTc            <  1:336675/101‑1 (MQ=255)
                                                                                           gCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTcc           >  1:3451549/1‑101 (MQ=255)
                                                                                            ccGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTccc          <  1:2081875/101‑1 (MQ=255)
                                                                                                cGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCcg      >  2:1494075/1‑101 (MQ=255)
                                                                                                   tCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCgg   <  1:2403334/101‑1 (MQ=255)
                                                                                                    |                                                                                                      
GGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGG  >  NC_014034/1395449‑1395649

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