Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,668 T→C 6.1% intergenic (+480/‑451) 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,6680T→C 6.1% 33.5 ‑0.2 ‑0.1309intergenic (+480/‑451)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (9/10):  ref base (151/138):  total (160/149)
Fisher's exact test strand distribution p-value = 8.14e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 8.24e-01

Only 100 of 344 total aligned reads displayed.

TCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGC  >  NC_014034/1395568‑1395768
                                                                                                    |                                                                                                        
tcgCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCt                                                                                                          >  2:3146135/1‑101 (MQ=255)
 cgCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTg                                                                                                         >  1:4072822/1‑101 (MQ=255)
 cgCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTg                                                                                                         >  1:5727880/1‑101 (MQ=255)
   cTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCa                                                                                                       >  1:4901833/1‑101 (MQ=255)
   cTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCa                                                                                                       <  1:1946553/101‑1 (MQ=255)
    tCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCAt                                                                                                      <  1:3441969/101‑1 (MQ=255)
     ccccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATg                                                                                                     >  2:1788375/1‑101 (MQ=255)
      cccccGCGGGGGCGTGGATAGAAACCACCTTGCGGATTGGCACTTTCGCGGGGGGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGt                                                                                                    <  1:3093963/101‑1 (MQ=255)
       ccccGCGGGGGCGTGGATCTAAACCACCTTGCGTATTGTCACTTTCGCTGGGCGAAGCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                   <  1:844480/101‑1 (MQ=38)
       ccccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTa                                                                                                   <  1:3637709/101‑1 (MQ=255)
         ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                 <  2:1483754/101‑1 (MQ=255)
         ccGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGt                                                                                                 >  2:136214/1‑101 (MQ=255)
           gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                               >  2:3193757/1‑101 (MQ=255)
           gcgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcg                                                                                               >  2:5631264/1‑101 (MQ=255)
            cgGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcgg                                                                                              >  2:3198446/1‑101 (MQ=255)
             gggggCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTcggc                                                                                             >  1:4751729/1‑101 (MQ=255)
                 gCGTGGATCGAATCCACCTTGCGGATTGTCACGTTCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.Gc                                                                                      <  2:5559917/101‑1 (MQ=255)
                  cGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCa                                                                                     <  1:2558068/101‑1 (MQ=255)
                   gTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAg                                                                                    >  1:5275284/1‑101 (MQ=255)
                    tGGATCGAACCCACCTTGCGGATTGTCCCTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                   <  1:1840623/101‑1 (MQ=255)
                    tGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGt                                                                                   >  1:1478976/1‑101 (MQ=255)
                     ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGtt                                                                                  <  2:575428/101‑1 (MQ=255)
                     ggATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACGTCGGAATGTAGTCGGCG..G.GCAGtt                                                                                  >  2:5400493/1‑101 (MQ=38)
                      gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGGTCGAAACATCTGCATGTAGTCGGCG..G.GCAGGTc                                                                                 >  2:5043884/1‑101 (MQ=255)
                      gATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTc                                                                                 >  1:2943090/1‑101 (MQ=255)
                       aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                                <  2:2467189/101‑1 (MQ=255)
                       aTCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCg                                                                                <  2:2566338/101‑1 (MQ=255)
                         cGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.g                                                                             >  1:4450568/1‑101 (MQ=255)
                           aaaCCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.Gcc                                                                           >  1:5571891/1‑101 (MQ=255)
                            aaCCACCTTGAGGATTGTCACTTTCGCGGGGCGAAGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCg                                                                          >  1:3332393/1‑101 (MQ=255)
                               ccccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAg                                                                       <  1:2407563/99‑1 (MQ=255)
                                gccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTGCGA.GCCGCAGt                                                                      <  2:2179038/100‑1 (MQ=255)
                                cccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGCCGGCG..G.GCAGTTCGA.GCCGCAGt                                                                      <  2:3377254/100‑1 (MQ=255)
                                 ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTc                                                                     <  2:1151501/101‑1 (MQ=255)
                                 ccTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTc                                                                     <  2:1874184/101‑1 (MQ=255)
                                   ttGGGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGc                                                                   <  1:5427231/101‑1 (MQ=255)
                                      cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                <  1:4120155/101‑1 (MQ=255)
                                      cGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcc                                                                >  2:4567470/1‑101 (MQ=255)
                                         aTTGTCCCTTCCGCGGGGCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTccccc                                                             <  2:3329068/101‑1 (MQ=255)
                                          ttGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTcccccc                                                            >  1:2191126/1‑101 (MQ=255)
                                           tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcg                                                           <  2:419156/101‑1 (MQ=255)
                                           tGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCACCCCCcg                                                           >  1:978275/1‑101 (MQ=255)
                                            gTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgc                                                          >  2:5311451/1‑101 (MQ=255)
                                            gTCACTTTCGC.GGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                         >  2:398376/1‑101 (MQ=255)
                                             tCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCcgcg                                                         <  1:3976555/101‑1 (MQ=255)
                                               acttcggcgGGGCCCAGCGTCCCTCCGCCCGCGGGGGCGAGGATCCAAACCTCGGCATGTAGCCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggg                                                       <  2:2099077/95‑1 (MQ=255)
                                                 tttCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCggggg                                                     >  1:5188986/1‑101 (MQ=255)
                                                     tggtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGt                                                  <  2:64120/91‑1 (MQ=2)
                                                     gcgGGGCGAGGCGTCGCGCCCCCCGCGGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                 <  1:3193946/101‑1 (MQ=255)
                                                     gcgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTg                                                 <  1:2892873/101‑1 (MQ=255)
                                                      cgGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.TCCGCAGTCGCTCCCCCCGCGGGGGCGTgg                                                <  2:3527274/101‑1 (MQ=255)
                                                       gtgccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTgg                                                <  1:268315/93‑1 (MQ=2)
                                                        gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                              >  2:4215510/1‑101 (MQ=255)
                                                        gggCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                              >  2:906044/1‑101 (MQ=255)
                                                        gggCGAAGCGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  2:3592893/101‑1 (MQ=255)
                                                         ggCGAGGCGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATc                                             <  1:1019576/101‑1 (MQ=255)
                                                         gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                              >  1:3495146/7‑101 (MQ=2)
                                                         gccgctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGAt                                              <  2:2835973/95‑1 (MQ=2)
                                                           cGAAGCGTCGCTCCCCCCGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                           <  1:5695492/101‑1 (MQ=255)
                                                            gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGGTCGaa                                          >  1:1210895/1‑101 (MQ=255)
                                                            gAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          >  2:5449150/1‑101 (MQ=255)
                                                             ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATGGCTATGGCGGCAGTGTCA.GGTTTTCCA.GCGGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          <  1:3548070/99‑1 (MQ=2)
                                                             ctgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                          <  2:4929980/99‑1 (MQ=2)
                                                              tgCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACCCCGAACCTGTCGGCG..GCGAACACCGATACCGC.GTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                         >  2:2704099/2‑101 (MQ=2)
                                                                 gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTgg                                     >  2:3816160/1‑101 (MQ=255)
                                                                  tCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGt                                    >  1:4182080/1‑101 (MQ=255)
                                                                   cGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTa                                   >  2:4349504/1‑101 (MQ=255)
                                                                    gccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                  <  2:2224101/98‑1 (MQ=255)
                                                                    gCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                  >  1:1979828/1‑101 (MQ=255)
                                                                    cccccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTAt                                  <  2:2904973/98‑1 (MQ=255)
                                                                         ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGGATGGATc                             >  2:4017561/1‑101 (MQ=255)
                                                                          cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                            <  1:1788474/101‑1 (MQ=255)
                                                                          cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                            >  2:474148/1‑101 (MQ=255)
                                                                          cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                            <  1:4508461/101‑1 (MQ=255)
                                                                          cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCt                            <  1:5369011/101‑1 (MQ=255)
                                                                             gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGgc                         >  1:5070293/1‑101 (MQ=255)
                                                                                  ggCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCaaa                    >  2:1243562/1‑101 (MQ=255)
                                                                                   gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                   >  1:4829436/1‑101 (MQ=255)
                                                                                   gCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAAc                   <  1:5631556/101‑1 (MQ=255)
                                                                                      tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                >  2:1574342/1‑101 (MQ=255)
                                                                                      tGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGGAGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGa                >  2:5002243/1‑101 (MQ=255)
                                                                                        gATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGaaa              >  2:3893119/1‑101 (MQ=255)
                                                                                         aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg             >  2:195141/1‑101 (MQ=255)
                                                                                         aTCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAg             >  1:624109/1‑101 (MQ=255)
                                                                                           cGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc           <  2:413936/101‑1 (MQ=255)
                                                                                           cGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTc           >  2:734116/1‑101 (MQ=255)
                                                                                            gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt          >  2:664709/1‑101 (MQ=255)
                                                                                            gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt          >  1:201600/1‑101 (MQ=255)
                                                                                            gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt          >  2:5222899/1‑101 (MQ=255)
                                                                                            gAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCt          <  1:5218408/101‑1 (MQ=255)
                                                                                              aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt        >  1:3056047/1‑101 (MQ=255)
                                                                                              aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt        >  2:4846339/1‑101 (MQ=255)
                                                                                              aaCATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt        >  1:4311345/1‑101 (MQ=255)
                                                                                              aaCATCTGCATGTAGGCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgt        <  1:1067295/101‑1 (MQ=255)
                                                                                                 aTCTGCAGGTAGTCGGCG..G.GCAGTCCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgt     <  2:4097254/101‑1 (MQ=255)
                                                                                                  tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc    <  2:3016019/101‑1 (MQ=255)
                                                                                                  tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc    >  1:1564862/1‑101 (MQ=255)
                                                                                                  tCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc    <  2:4495295/101‑1 (MQ=255)
                                                                                                  tCTGCATGTAGGCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTgtcgtc    <  1:3792757/101‑1 (MQ=255)
                                                                                                   cTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGgcgtcg   >  2:2000575/1‑101 (MQ=255)
                                                                                                    |                                                                                                        
TCGCTCCCCCCGCGGGGGCGTGGATCGAAACCACCTTGCGGATTGTCACTTTCGCGGGGCGAAGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCG..G.GCAGTTCGA.GCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTGGTATGGATCTGGCGCAAACGGAAAGTCTGTCGTCGC  >  NC_014034/1395568‑1395768

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