Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,761 +A 10.6% intergenic (+573/‑358) 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,7611.→A10.6% 115.9 ‑1.4 ‑0.8274intergenic (+573/‑358)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (11/18):  ref base (122/123):  total (133/141)
Fisher's exact test strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.72e-01

Only 100 of 306 total aligned reads displayed.

AAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1395661‑1395861
                                                                                                           |                                                                                                    
aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                                       <  2:4403571/101‑1 (MQ=255)
aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                                       >  1:3585406/1‑101 (MQ=255)
aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g                                                                                                       <  1:89249/101‑1 (MQ=255)
 aaCATCTGCATGTGGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.t                                                                                                     <  1:2106180/101‑1 (MQ=255)
  aCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tc                                                                                                    >  1:1928921/1‑101 (MQ=255)
   cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcg                                                                                                   >  1:2817565/1‑101 (MQ=255)
   cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcg                                                                                                   <  1:4909450/101‑1 (MQ=255)
    ctCTGCATGTGGCGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                                  <  2:1548907/100‑1 (MQ=255)
      cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                                                >  1:3700067/1‑101 (MQ=255)
      cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                                                >  1:2398007/1‑101 (MQ=255)
       tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                                               >  2:2260291/1‑101 (MQ=255)
            gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                                          >  1:3087414/1‑101 (MQ=255)
            gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                                          <  2:4599327/101‑1 (MQ=255)
              cggcggcggGCAGTTCGACCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACT.GGTTTGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                                        <  1:798152/98‑1 (MQ=38)
               gTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcg                                                                                       <  1:3210529/101‑1 (MQ=255)
                  ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                                                    >  2:2121547/1‑101 (MQ=255)
                     gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACA.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCggggg                                                                                 >  1:4274480/1‑101 (MQ=255)
                          gTTCGAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                                            <  2:4180378/101‑1 (MQ=255)
                            tCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          >  1:2637204/1‑101 (MQ=255)
                            ccGAGCCGCCGTCGCCCCCCCCGCGGGGGCGTTGATCGAAACT.GGTATGGATCTG..GCCCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                          <  1:2138346/100‑1 (MQ=38)
                              gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                        <  1:2328698/101‑1 (MQ=255)
                              gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                        <  1:4511956/101‑1 (MQ=255)
                              cagattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGa                                                                           <  2:815732/95‑1 (MQ=2)
                               cgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                         <  2:5018250/94‑1 (MQ=0)
                                gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                        <  1:3977772/95‑1 (MQ=0)
                                gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                         >  2:1027640/5‑101 (MQ=2)
                                  cGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  1:69988/101‑1 (MQ=255)
                                   tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCgca                                                                      >  2:2672616/2‑99 (MQ=2)
                                   tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                                      >  2:3575117/2‑101 (MQ=2)
                                   ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                                                     <  2:4923984/98‑1 (MQ=2)
                                    cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCt                                                                  >  1:4532405/1‑101 (MQ=255)
                                     ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                                   >  2:3711699/2‑101 (MQ=2)
                                     aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCTCCGCGGGGGCGTGGATCGAAAc                                                                    <  2:515778/101‑1 (MQ=2)
                                     aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                                    <  2:4756481/101‑1 (MQ=2)
                                     aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                                                 <  2:1488763/101‑1 (MQ=255)
                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcac                                                                  <  1:2535894/101‑3 (MQ=2)
                                       tCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcg                                                               >  2:123139/1‑101 (MQ=255)
                                        cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                                                 <  1:1849072/101‑4 (MQ=2)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                                                <  1:3638760/101‑5 (MQ=2)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                               >  2:3416623/1‑98 (MQ=2)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                                               >  1:4871950/1‑98 (MQ=2)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggt                                                             <  2:4565843/101‑2 (MQ=255)
                                         gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                                             >  2:4992675/1‑101 (MQ=255)
                                          cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                               <  2:4893434/101‑6 (MQ=2)
                                          cTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggcg                                                            >  1:2267527/1‑101 (MQ=255)
                                           tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgtatgag                                                             <  1:2066794/101‑8 (MQ=2)
                                           tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                                             <  1:3037247/101‑6 (MQ=2)
                                           cccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                                              <  2:2086059/100‑7 (MQ=0)
                                           cccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                                              <  2:3894399/100‑7 (MQ=0)
                                            gatctcgcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                                          <  1:5017625/96‑1 (MQ=255)
                                            ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGGGGATCGAAACCTGCGGCGTg                                                          >  1:2194961/1‑101 (MQ=255)
                                             cccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                                                            >  2:3771739/1‑93 (MQ=2)
                                               cccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTTAGt                                                       <  2:4274163/101‑1 (MQ=255)
                                                ccGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                                                      <  1:133561/101‑1 (MQ=255)
                                                  gcgGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGc                                                    <  2:3127920/101‑1 (MQ=255)
                                                     ggggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                                                 >  2:1964466/1‑101 (MQ=255)
                                                     ggggCGTGGATCGAAACT.GGTACGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAg                                                 >  2:2219794/1‑101 (MQ=255)
                                                      gggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAgg                                                >  2:2848348/1‑100 (MQ=255)
                                                       ggCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTg                                               <  2:3359435/101‑1 (MQ=255)
                                                         cGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                                             <  2:2696841/101‑1 (MQ=255)
                                                           tGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                                           <  2:4009932/101‑1 (MQ=255)
                                                           tGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                                           <  1:274270/101‑1 (MQ=255)
                                                            ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                                          <  1:2685758/101‑1 (MQ=255)
                                                            ggATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                                          >  1:4601131/1‑101 (MQ=255)
                                                              aTCGAAAAT.GGAATGGATCAG..GCGCAAACGG..TAAGTCT.G.TCGTCGGTCCCCCCGCGGGGGCGTGGATCGAAGCACGCGACGGGAGGTCCCAGCGATAGagg                                        >  2:2443276/1‑99 (MQ=16)
                                                               tCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc                                       <  1:4627997/101‑1 (MQ=255)
                                                                 gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg                                     >  2:2298444/1‑101 (MQ=255)
                                                                 gAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg                                     >  1:3327564/1‑101 (MQ=255)
                                                                  aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCgg                                    >  1:1093676/1‑100 (MQ=255)
                                                                  aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt                                    >  2:1840929/1‑101 (MQ=255)
                                                                  aaaCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGt                                    >  2:3326507/1‑101 (MQ=255)
                                                                    aCT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg                                  <  2:2958186/101‑1 (MQ=255)
                                                                     cT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa                                 >  2:4311946/1‑101 (MQ=255)
                                                                      t.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCAGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc                                >  1:4332930/1‑101 (MQ=255)
                                                                             ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGTCGTCGCt                          >  2:4501529/1‑101 (MQ=255)
                                                                             ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCt                          >  1:3807805/1‑101 (MQ=255)
                                                                             ggATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCCAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCt                          <  2:4362772/101‑1 (MQ=255)
                                                                               aTCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcc                        <  1:3425773/101‑1 (MQ=255)
                                                                                  tG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTccccc                     <  2:3862515/101‑1 (MQ=255)
                                                                                   g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                    >  2:4673375/1‑101 (MQ=255)
                                                                                   g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                    <  1:929284/101‑1 (MQ=255)
                                                                                   g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc                    <  1:4949833/101‑1 (MQ=255)
                                                                                         cAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCgg                >  2:3678384/1‑101 (MQ=255)
                                                                                         cAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCgg                >  1:1500376/1‑101 (MQ=255)
                                                                                           aaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTAAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCgggg              <  1:4029958/101‑1 (MQ=255)
                                                                                            aCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCggggg             >  2:4205555/1‑101 (MQ=255)
                                                                                             cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc            >  1:1230956/1‑101 (MQ=255)
                                                                                             cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc            <  1:1008604/101‑1 (MQ=255)
                                                                                             cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc            >  1:1614343/1‑101 (MQ=255)
                                                                                             cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc            >  1:4783394/1‑101 (MQ=255)
                                                                                             cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc            <  1:3482243/101‑1 (MQ=255)
                                                                                                  aaaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTg         <  1:2678950/101‑1 (MQ=255)
                                                                                                   aaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTgg        >  1:3237075/1‑101 (MQ=255)
                                                                                                      tCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATc     >  2:647316/1‑101 (MQ=255)
                                                                                                       cT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCg    <  1:4212407/101‑1 (MQ=255)
                                                                                                        t.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGa   >  1:2872901/1‑101 (MQ=255)
                                                                                                          g.tcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATCGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGTTCGaa  >  2:2015318/1‑101 (MQ=255)
                                                                                                          g.tcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa  >  1:126455/1‑101 (MQ=255)
                                                                                                          g.tcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa  >  1:168761/1‑101 (MQ=255)
                                                                                                          g.tcgtcGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGaa  >  1:3863081/1‑101 (MQ=255)
                                                                                                           |                                                                                                    
AAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTGGATCGAA  >  NC_014034/1395661‑1395861

Base quality scores:  ATCG  < 3 ≤  ATCG  < 21 ≤  ATCG  < 29 ≤  ATCG  < 34 ≤  ATCG  < 41 ≤  ATCG