Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,536 C→A 8.8% intergenic (+348/‑583) 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,5360C→A 8.8% 72.3 ‑2.4 ‑0.0354intergenic (+348/‑583)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/22):  ref base (148/175):  total (157/197)
Fisher's exact test strand distribution p-value = 8.86e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.66e-01

Only 100 of 385 total aligned reads displayed.

CGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCG  >  NC_014034/1395436‑1395636
                                                                                                    |                                                                                                      
cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACAtctc                                                                                                        <  1:1893708/101‑1 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                       >  1:34010/1‑101 (MQ=255)
 gCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCg                                                                                                       <  1:2899007/101‑1 (MQ=255)
  cTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCgg                                                                                                      >  1:3615483/1‑101 (MQ=255)
               ggCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.t                                                                                        <  2:3676538/101‑1 (MQ=255)
                  gTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCa                                                                                     <  2:1526942/101‑1 (MQ=255)
                    ggATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGc                                                                                   >  2:3615592/1‑101 (MQ=255)
                     gATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGcc                                                                                  >  2:2400467/1‑101 (MQ=255)
                        cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGGCGTCGCTCG.TCCAGCCGG.c                                                                              >  2:1765393/1‑101 (MQ=255)
                        cGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.c                                                                              >  2:680381/1‑101 (MQ=255)
                              aTCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.CCCAGCCGG.CCTGCtc                                                                        >  1:1076904/1‑101 (MQ=255)
                                cgccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgt                                                                      >  2:2645470/1‑101 (MQ=255)
                                  ccgcTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCtcgtcg                                                                    <  2:2879717/101‑1 (MQ=255)
                                      tttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcc                                                                >  2:3948998/1‑101 (MQ=255)
                                       ttGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCCCTccc                                                               >  1:4972627/1‑101 (MQ=255)
                                        tGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccc                                                              >  1:1677328/1‑101 (MQ=255)
                                          gCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTcccccc                                                            <  2:581116/101‑1 (MQ=255)
                                               cGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCggg                                                       <  2:34000/101‑1 (MQ=255)
                                                     gCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCCCCCCCCCCGGGGGCggg                                                 >  2:2161716/1‑99 (MQ=255)
                                                      gccctcgcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGCCTATG.GCTGGCCGC.CCCGCCCCTCGCTCCCCCCGCGGGGGCGTgg                                                >  2:2014018/11‑101 (MQ=2)
                                                      cccctcccccgtCGCCCCCCCCGCGGTGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:3265946/91‑1 (MQ=0)
                                                      cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg                                                <  2:3051784/101‑1 (MQ=255)
                                                      cTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTgg                                                >  1:4234010/1‑101 (MQ=255)
                                                      agcggctgcggtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:3526510/91‑1 (MQ=0)
                                                       gcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCACCCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGa                                               <  2:1661797/92‑1 (MQ=0)
                                                       gcggctgcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGa                                               <  1:477855/92‑1 (MQ=0)
                                                        ggccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGAt                                              >  1:4833911/9‑101 (MQ=2)
                                                        ggTAACGTGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  1:1795421/101‑1 (MQ=255)
                                                         gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  1:1988851/94‑1 (MQ=2)
                                                         gccggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  1:2483963/94‑1 (MQ=2)
                                                         gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGTATc                                             <  2:4745741/101‑1 (MQ=255)
                                                         gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  2:3891191/101‑1 (MQ=255)
                                                         gTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGGCGCTCCCCCCGCGGGGGCGGGGATc                                             >  2:2830023/1‑101 (MQ=255)
                                                          tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                            <  1:2836702/101‑1 (MQ=255)
                                                          tAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGGGGGTCg                                            >  2:929179/1‑101 (MQ=255)
                                                           cggtggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                           <  1:490431/96‑1 (MQ=2)
                                                           aaCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGGGGATcaa                                           >  2:2212785/1‑99 (MQ=255)
                                                            gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          <  1:3695411/97‑1 (MQ=2)
                                                            gcacgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                          >  1:4421915/5‑101 (MQ=2)
                                                            aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          >  1:433048/1‑101 (MQ=255)
                                                            aCGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          <  1:2293623/101‑1 (MQ=255)
                                                             gcggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         <  1:3846625/98‑1 (MQ=2)
                                                             cGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         <  2:3189479/101‑1 (MQ=255)
                                                              tggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                        <  1:3772024/99‑1 (MQ=2)
                                                              gtgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGGGGGTCGAAac                                        >  1:895680/1‑101 (MQ=255)
                                                               tgtCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacc                                       <  2:1931878/101‑1 (MQ=255)
                                                               ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                       >  1:4161174/2‑100 (MQ=2)
                                                               ggtCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAact                                       <  1:223518/100‑2 (MQ=2)
                                                                gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                      >  2:4192343/1‑99 (MQ=2)
                                                                gtCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                      >  1:4009300/1‑99 (MQ=2)
                                                                gtCGCCCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactg                                      <  1:4639615/101‑3 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgc                                     <  1:5053765/101‑4 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacccc                                     <  1:2411224/101‑3 (MQ=2)
                                                                 tCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                     >  2:310388/1‑98 (MQ=2)
                                                                 gccccccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGCCCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgc                                     <  2:4448557/96‑4 (MQ=2)
                                                                 gcGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactat                                     <  2:1057343/100‑4 (MQ=2)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                    <  2:5031607/101‑1 (MQ=255)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccacc                                    >  2:2714598/1‑101 (MQ=255)
                                                                  cGCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatg                                    >  2:2544048/1‑97 (MQ=2)
                                                                   gccccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgt                                   <  1:4462868/98‑6 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgc                                   <  1:3928218/101‑5 (MQ=2)
                                                                   gCTCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgc                                   >  1:2141749/1‑96 (MQ=2)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                 <  1:3358851/101‑1 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTg                                 <  1:2974364/101‑1 (MQ=255)
                                                                     tCCCCCCGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                 >  1:1176760/1‑94 (MQ=2)
                                                                     cccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgc                                 <  1:3062797/100‑8 (MQ=2)
                                                                     cccccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctg                                 <  2:1410700/100‑8 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGc                                <  1:1542028/101‑1 (MQ=255)
                                                                      ccccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcg                                <  1:1300738/101‑9 (MQ=2)
                                                                      ccccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggc                                >  1:780972/1‑94 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCg                               >  2:1718024/1‑101 (MQ=255)
                                                                       cccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgtgcgc                               >  1:2126314/1‑92 (MQ=2)
                                                                       cccccGCGGGGGCGTGGATCGAAACAATGCGGCGCGGGC.CTAT.GGCTGACGGACCCGCCCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactgcgcgcgc                               >  2:3498382/1‑92 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcgt                              >  1:2193081/1‑92 (MQ=2)
                                                                        ccccGCGGGGGCGTGGATCGAAACA.CCCGGCGCCGCAGCGAG.TGCGTCTGGAACAGGATGTCGCTCCCCCCGCGGGGGCGTGGATCGAAactatgctggct                              <  2:1590595/101‑11 (MQ=2)
                                                                         cccGCGGGGGCGTGGATCGAAACAAAAAGCCGCCGACAAACGCGGCAACCGG.CCCG.GCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAaccccgcggcgtc                             <  2:3830578/101‑11 (MQ=2)
                                                                            gcgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTg                          <  1:3502072/101‑1 (MQ=255)
                                                                             cgGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGt                         <  2:1091778/101‑1 (MQ=255)
                                                                              gggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTc                        >  2:1499608/1‑101 (MQ=255)
                                                                               ggggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCa                       <  2:4410284/101‑1 (MQ=255)
                                                                                gggCGTGCATGGAAACCTCTCGTTGCCGTCGCTCG.GCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCAc                      <  2:1209283/101‑1 (MQ=37)
                                                                                 ggCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACCACCTTGCGGATTGTCACt                     <  1:4543985/101‑1 (MQ=255)
                                                                                  gCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACtt                    <  2:454096/101‑1 (MQ=255)
                                                                                   cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTTTCACttt                   >  2:4627674/1‑101 (MQ=255)
                                                                                   cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                   >  1:4026601/1‑101 (MQ=255)
                                                                                   cGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                   <  1:1018856/101‑1 (MQ=255)
                                                                                   cGTGGATCAAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACttt                   <  1:1507009/101‑1 (MQ=255)
                                                                                    gTGGATCGAAACATCTCGGTGCCGTCGCTCN.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                  <  1:602455/101‑1 (MQ=255)
                                                                                    gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                  >  2:4332590/1‑101 (MQ=255)
                                                                                    gTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTc                  <  2:3254458/101‑1 (MQ=255)
                                                                                     tGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTcg                 <  2:2289565/101‑1 (MQ=255)
                                                                                      ggATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCAGCTTGCGGATTGTCCCTTTcgc                >  2:2106032/1‑101 (MQ=255)
                                                                                          cGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgggg            >  1:1613284/1‑101 (MQ=255)
                                                                                          cGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCgggg            <  2:2484409/101‑1 (MQ=255)
                                                                                           gAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGc           <  1:1718154/101‑1 (MQ=255)
                                                                                             aaCATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGa         <  2:2982933/101‑1 (MQ=255)
                                                                                               cATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAg       <  1:640608/101‑1 (MQ=255)
                                                                                                aTCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGc      <  2:4972297/101‑1 (MQ=255)
                                                                                                 tctcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGGGAAGCg     >  1:4566164/1‑101 (MQ=255)
                                                                                                  ctcGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGt    >  2:4029655/1‑101 (MQ=255)
                                                                                                    |                                                                                                      
CGCTCCCCCCGCGGGGGCGTGGATCGAAACATCGCCGCTTTGGCGGTCGACAGGCTGGTAACGTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTCGGTGCCGTCGCTCG.TCCAGCCGG.CCTGCTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCG  >  NC_014034/1395436‑1395636

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