Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,698 Δ1 bp 8.0% intergenic (+510/‑421) 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,6980A→. 8.0% 97.5 ‑1.5 ‑0.0447intergenic (+510/‑421)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (13/23):  ref base (198/213):  total (211/236)
Fisher's exact test strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.96e-01

Only 100 of 478 total aligned reads displayed.

CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798
                                                                                                      |                                                                                                       
  cccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGt                                                                                                       <  2:437371/100‑1 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  2:2877296/1‑101 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  2:2435596/1‑101 (MQ=255)
    cTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCg                                                                                                     >  2:1080628/1‑101 (MQ=255)
       gCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTc                                                                                                  >  1:1017147/1‑101 (MQ=255)
           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                              <  1:8596316/101‑1 (MQ=255)
           aTTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTccccc                                                                                              >  1:2764870/1‑101 (MQ=255)
             tGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                                                            <  2:1236883/101‑1 (MQ=255)
              gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                                                           <  2:2113328/101‑1 (MQ=255)
                cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCgg                                                                                         >  1:1731193/1‑101 (MQ=255)
                cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GACGCAGTCGCTCCCCCCGCgg                                                                                         >  1:5642593/1‑101 (MQ=255)
                cACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACATCTGCATGCAGTCGGCGG.GCAGCTCTA.GCCGCAGTCGCTCCCGCCGCgg                                                                                         >  1:2979398/1‑101 (MQ=38)
                 ccTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                                                        <  2:873687/100‑1 (MQ=255)
                    ttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.CCCGCTGGA.GCCGCAGGCGCTTCCCCCGCGGCGGc                                                                                     >  2:7347798/1‑101 (MQ=25)
                        ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  2:5734907/10‑101 (MQ=2)
                        ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  1:1189167/10‑101 (MQ=2)
                        ggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTAGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTg                                                                                  >  1:6473577/10‑101 (MQ=2)
                         gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                                                 <  2:1471612/93‑1 (MQ=2)
                         ggggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGa                                                                                >  2:1935834/1‑101 (MQ=255)
                          gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                               <  1:3738050/101‑1 (MQ=255)
                           ggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:8476112/101‑1 (MQ=255)
                            gCGAAGCGCCGCTCCCCCCGCGGGGGCGTGGATCAAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                             <  2:8566409/101‑1 (MQ=255)
                            ccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATc                                                                              <  1:8118958/96‑1 (MQ=2)
                              gctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCGGCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                            >  2:3243457/4‑101 (MQ=2)
                              gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                           >  1:3855902/1‑101 (MQ=255)
                                  cGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.g                                                                      >  2:5217229/1‑101 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGTAACT.gg                                                                     >  2:8493761/1‑101 (MQ=255)
                                   gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.gg                                                                     >  2:4695163/1‑101 (MQ=255)
                                    ccGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGt                                                                    <  2:7727175/100‑1 (MQ=255)
                                     cGCCCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATTTAGTCGGCGT.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTa                                                                   <  1:1574047/101‑1 (MQ=37)
                                      gccccccccGCGGGGGGGGGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTAt                                                                  <  1:2322296/98‑1 (MQ=255)
                                       cccccccgGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 <  2:1565832/99‑1 (MQ=255)
                                       cTCCCCCCGGGGGGGCGTGGAGCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 <  2:504086/101‑1 (MQ=255)
                                       cTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATg                                                                 <  1:61707/101‑1 (MQ=255)
                                         ccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGa                                                               <  1:6704232/101‑1 (MQ=255)
                                         ccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGa                                                               >  1:4099766/1‑101 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  2:179784/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  2:6674320/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  1:437804/101‑1 (MQ=255)
                                            cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                            <  1:4520389/101‑1 (MQ=255)
                                             ccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTg                                                           <  1:2033584/101‑1 (MQ=255)
                                                  ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCa                                                      <  2:6494688/101‑1 (MQ=255)
                                                    ggCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCaaa                                                    <  1:8744175/101‑1 (MQ=255)
                                                      cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGGGTGGATCGAAACG.GGTATGAGTTCGGCGCaaagg                                                  >  2:8966096/1‑99 (MQ=37)
                                                      cGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACg                                                  <  1:2565719/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTGGTCGGCGG.GCGGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        <  1:9137598/101‑1 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  1:7403700/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCt                                        >  1:5123917/1‑101 (MQ=255)
                                                              gAAACATCTGCATGTAGTCGGCGG.GCAGTTCAA.GCCGCAGTCGCACCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTcg                                        >  1:8971841/1‑100 (MQ=255)
                                                               aaaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTg                                       >  2:4446939/1‑101 (MQ=255)
                                                                aaCATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgt                                      <  1:1080646/101‑1 (MQ=255)
                                                                 aCATCGGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtc                                     <  1:6592999/101‑1 (MQ=255)
                                                                  cATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcg                                    >  2:5587481/1‑101 (MQ=255)
                                                                   cttggCATGTAGTCGGCGG.GCAGTTGGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  2:80515/97‑1 (MQ=255)
                                                                   ctCTGCATGTGGCCGGCGG.GCGGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTgtcgt                                   <  2:2240799/100‑1 (MQ=255)
                                                                   aTCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACT.GGTATGGGTCTGGCGCAAACGG..AAAGTCTggcgg                                   >  2:3725528/1‑100 (MQ=255)
                                                                     cTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGtcgtcg                                 >  2:3537640/1‑101 (MQ=255)
                                                                         aTGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcc                             <  2:7953520/101‑1 (MQ=255)
                                                                         aTGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcc                             <  1:5941947/101‑1 (MQ=255)
                                                                         aTGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCGGTGGTCGCTcc                             >  2:4936558/1‑101 (MQ=255)
                                                                          tGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTccc                            >  2:7513387/1‑101 (MQ=255)
                                                                           gTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccc                           >  2:6874540/1‑101 (MQ=255)
                                                                             ggTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTcccccc                         <  2:6050490/100‑1 (MQ=255)
                                                                              gTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCcg                        <  1:7699651/101‑1 (MQ=255)
                                                                                   ggg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCgggg                   <  2:1017129/100‑1 (MQ=255)
                                                                                   ggg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCgggg                   <  1:6257551/100‑1 (MQ=255)
                                                                                    gg.gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCggggg                  >  1:1899628/1‑101 (MQ=255)
                                                                                       gCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCg                >  1:5351869/1‑101 (MQ=255)
                                                                                         aGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTg              >  1:3554439/1‑101 (MQ=255)
                                                                                         aGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCggg              >  1:7216967/1‑99 (MQ=255)
                                                                                          gtgcGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             <  2:6172622/98‑1 (MQ=255)
                                                                                          gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             >  2:2686860/1‑101 (MQ=255)
                                                                                          gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             >  1:6840928/1‑101 (MQ=255)
                                                                                          gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             >  2:8333411/1‑101 (MQ=255)
                                                                                          gTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTgg             >  2:5496538/1‑101 (MQ=255)
                                                                                           ttCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGa            >  1:3214051/1‑101 (MQ=255)
                                                                                           ttCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGa            <  1:7513711/101‑1 (MQ=255)
                                                                                            tCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGAt           <  1:5044304/101‑1 (MQ=255)
                                                                                              gA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:3855793/101‑1 (MQ=255)
                                                                                              gA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:2776266/101‑1 (MQ=255)
                                                                                              caggattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTgg             <  2:9119184/94‑1 (MQ=2)
                                                                                               gggCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:5131107/99‑1 (MQ=255)
                                                                                               aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGa            >  1:287577/7‑101 (MQ=2)
                                                                                               aagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGa            >  1:4784935/7‑101 (MQ=2)
                                                                                                agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGAt           <  1:7020256/96‑1 (MQ=2)
                                                                                                 gcagCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  1:1556136/98‑1 (MQ=255)
                                                                                                 gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       >  1:1401975/1‑101 (MQ=255)
                                                                                                  ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa      <  2:5329308/101‑1 (MQ=255)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  1:2129890/98‑1 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         >  1:5813000/4‑101 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         >  1:7143125/4‑101 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         <  2:7579598/98‑1 (MQ=2)
                                                                                                  attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCg         >  2:7186577/4‑101 (MQ=2)
                                                                                                   gtcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGa        <  2:2533321/99‑1 (MQ=2)
                                                                                                    tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       <  1:6364818/100‑1 (MQ=2)
                                                                                                    tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa       >  2:6326193/2‑101 (MQ=2)
                                                                                                     cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGGGGATCGaaa      >  2:3148402/1‑101 (MQ=2)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCGCCGATCTGGGTAAGGCCGGTGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc     >  1:7880537/1‑101 (MQ=2)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg  >  2:5907301/1‑101 (MQ=255)
                                                                                                      aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGGGGGTCGCAAcccc  >  2:4223438/1‑99 (MQ=255)
                                                                                                      |                                                                                                       
CCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGG.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTGGCGCAAACGG..AAAGTCTGTCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTG  >  NC_014034/1395598‑1395798

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