Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,754 A→C 10.3% intergenic (+566/‑365) 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,7540A→C10.3% 72.1 ‑2.8 ‑0.0299intergenic (+566/‑365)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (10/21):  ref base (136/132):  total (146/153)
Fisher's exact test strand distribution p-value = 5.87e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

Only 100 of 309 total aligned reads displayed.

TGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1395654‑1395854
                                                                                                          |                                                                                                      
tGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..a                                                                                                        <  2:1743662/101‑1 (MQ=255)
  gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..aaa                                                                                                      <  1:472489/101‑1 (MQ=255)
   aTCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAg                                                                                                     >  1:80252/1‑101 (MQ=255)
   aTCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAg                                                                                                     >  2:3603562/1‑101 (MQ=255)
    tCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGt                                                                                                    >  1:1122675/1‑101 (MQ=255)
      gAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCt                                                                                                  <  2:1869812/101‑1 (MQ=255)
       aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g                                                                                                >  2:1691194/1‑101 (MQ=255)
          cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g.tcg                                                                                            >  2:4056125/1‑101 (MQ=255)
           ctcggCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                           <  2:4141045/97‑1 (MQ=255)
           agcTGCATGGAGTCGGCGGGGAGTTCGACCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g.tcgt                                                                                           <  2:3044277/99‑1 (MQ=255)
            tCTGCATGTGGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                                          <  1:2672198/101‑1 (MQ=255)
            tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                                          <  1:3172250/101‑1 (MQ=255)
               tgaTGTAGGGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGAGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCt                                                                                       <  1:1000128/99‑1 (MQ=255)
                  tGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccc                                                                                    <  1:3039365/101‑1 (MQ=255)
                   gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                                   >  2:1301508/1‑101 (MQ=255)
                    aattcgaagagCAGTACGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTAGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccccc                                                                                  <  2:2789949/91‑1 (MQ=37)
                     aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                                 >  1:2032090/1‑101 (MQ=255)
                     aGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccccc                                                                                 >  1:2850836/1‑101 (MQ=255)
                      gTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTccccccc                                                                                >  2:1499411/1‑100 (MQ=255)
                      gTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcg                                                                                >  2:819203/1‑101 (MQ=255)
                        cggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgcg                                                                              >  2:3784859/1‑101 (MQ=255)
                             ggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGc                                                                         <  2:2528726/101‑1 (MQ=255)
                             ggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGc                                                                         >  2:2836204/1‑101 (MQ=255)
                                aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCTTg                                                                      >  1:335673/1‑101 (MQ=255)
                                aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTg                                                                      >  1:3816219/1‑101 (MQ=255)
                                  ttCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                    >  2:4399679/1‑101 (MQ=255)
                                   tggAGCCGCAGTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  1:2309088/99‑1 (MQ=255)
                                    tgcgcgttggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGa                                                                    <  2:880848/92‑1 (MQ=2)
                                    cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  >  2:3196943/1‑101 (MQ=255)
                                    cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCCCGGGGGCGTGGATc                                                                  >  1:1668730/1‑101 (MQ=255)
                                     gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   >  1:1433172/9‑101 (MQ=0)
                                     gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  1:1048791/93‑1 (MQ=0)
                                      agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  2:488347/96‑1 (MQ=2)
                                      agattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   >  1:521557/6‑101 (MQ=2)
                                      aGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                >  2:3000738/1‑101 (MQ=255)
                                      aGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                >  2:2588555/1‑101 (MQ=255)
                                       gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCAAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 <  2:833335/95‑1 (MQ=0)
                                       gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               >  1:354525/1‑101 (MQ=255)
                                         ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                                                >  1:2736164/3‑101 (MQ=2)
                                          tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                               <  2:2010199/100‑1 (MQ=2)
                                           cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                           >  1:1203010/1‑101 (MQ=255)
                                           cAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                           <  2:1943643/101‑1 (MQ=255)
                                           cAGCCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                                           <  1:3981769/101‑1 (MQ=255)
                                            ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                                            <  2:3694717/100‑1 (MQ=2)
                                            aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                                             >  1:579377/1‑101 (MQ=2)
                                            aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                                          <  1:1383700/101‑1 (MQ=255)
                                             gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                                         >  1:3116456/1‑101 (MQ=255)
                                              tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcac                                                           <  2:1783854/101‑3 (MQ=2)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                                          >  2:1160083/1‑98 (MQ=2)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.ATGACCAGTT..CCGCAATCGG..CCGGGCTCGATCGGCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                                         >  1:762208/1‑99 (MQ=2)
                                               cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                                         <  2:70453/101‑3 (MQ=2)
                                                gccccccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                                         <  1:2289510/98‑5 (MQ=0)
                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAACtatg                                                        >  2:1589147/1‑98 (MQ=2)
                                                gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATATG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                                      >  2:3221327/1‑101 (MQ=255)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                        >  1:3061356/1‑96 (MQ=2)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                        >  1:1755366/1‑96 (MQ=2)
                                                 cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                                        <  1:2566513/101‑6 (MQ=2)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                                       >  2:3439305/1‑95 (MQ=2)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                                      >  2:3734982/1‑96 (MQ=2)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGt                                                    >  2:3096836/1‑101 (MQ=255)
                                                  tCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGt                                                    >  2:1103746/1‑101 (MQ=255)
                                                   ccccccGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                                   <  2:1445556/101‑1 (MQ=255)
                                                    cccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                                                     >  1:608148/1‑93 (MQ=2)
                                                      cccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcgcc                                                   <  2:2623561/101‑11 (MQ=2)
                                                       ccGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGtt                                               <  1:4072075/101‑1 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGc                                             >  1:710457/1‑101 (MQ=255)
                                                         gcgGGGGCGTGGATCGAAACT.GGTATG.GATCCG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGc                                             <  2:4089459/101‑1 (MQ=255)
                                                          cgGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTGGcc                                            >  2:2943048/1‑101 (MQ=255)
                                                           gggggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                                           >  2:4292067/1‑101 (MQ=255)
                                                           gggggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                                           >  1:1904105/1‑101 (MQ=255)
                                                           gggggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                                           >  1:3205556/1‑101 (MQ=255)
                                                                cGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                                      >  1:3453326/1‑101 (MQ=255)
                                                                   ggATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                                   <  2:577972/101‑1 (MQ=255)
                                                                     aTCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                                 >  1:4153850/1‑101 (MQ=255)
                                                                     aTCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                                 >  1:1654518/1‑101 (MQ=255)
                                                                     aTCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                                 <  2:1914804/101‑1 (MQ=255)
                                                                      tCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc                                <  1:974084/101‑1 (MQ=255)
                                                                        gAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCg                              >  1:1273450/1‑101 (MQ=255)
                                                                         aaaCT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTTGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCgg                             >  1:4401889/1‑100 (MQ=255)
                                                                          aaCT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTg                            <  1:4280754/101‑1 (MQ=255)
                                                                          aaCT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTg                            >  2:3314314/1‑101 (MQ=255)
                                                                               ggTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACg                        <  1:1538061/101‑1 (MQ=255)
                                                                                gTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCCGTGATAGACGCCGTGGACGt                       >  1:824162/1‑101 (MQ=255)
                                                                                   tG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGc                    >  1:2244579/1‑101 (MQ=255)
                                                                                      gATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTc                  >  2:2392933/1‑101 (MQ=255)
                                                                                      gATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTc                  <  2:3368022/101‑1 (MQ=255)
                                                                                        tCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTccc                >  2:3944024/1‑101 (MQ=255)
                                                                                        tCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTccc                <  1:1301605/101‑1 (MQ=255)
                                                                                        tCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTccc                <  2:1273357/101‑1 (MQ=255)
                                                                                         cTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccc               <  1:2939955/101‑1 (MQ=255)
                                                                                           g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcccccc             <  1:4121378/101‑1 (MQ=255)
                                                                                           g..gCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGGGAGTTGCCAGTGATAGACGCCGTGGGCGTCGCTcccccc             >  2:2102470/1‑101 (MQ=255)
                                                                                              gcgcAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCcg            <  2:2711784/101‑1 (MQ=255)
                                                                                                gcAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCcgcg          >  1:2075912/1‑101 (MQ=255)
                                                                                                   aaCGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCgggg       <  2:3116232/101‑1 (MQ=255)
                                                                                                     cGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGc     <  1:1864366/101‑1 (MQ=255)
                                                                                                      gg..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCg    >  2:1232612/1‑101 (MQ=255)
                                                                                                      gg..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCg    >  1:2151269/1‑101 (MQ=255)
                                                                                                       g..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGt   >  2:2690181/1‑101 (MQ=255)
                                                                                                          aaaGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTg  <  2:545206/101‑1 (MQ=255)
                                                                                                          |                                                                                                      
TGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCCCCCCGCGGGGGCGTG  >  NC_014034/1395654‑1395854

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