Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,741 C→A 15.7% intergenic (+553/‑378) 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,7410C→A15.7% 118.2 ‑0.7 ‑0.1292intergenic (+553/‑378)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/26):  ref base (122/124):  total (142/150)
Fisher's exact test strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.71e-01

Only 100 of 301 total aligned reads displayed.

CCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCC  >  NC_014034/1395641‑1395839
                                                                                                      |                                                                                                         
ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATc                                                                                                           <  2:1935384/101‑1 (MQ=255)
 cccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCt                                                                                                          <  1:2611943/101‑1 (MQ=255)
    gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..gc                                                                                                     >  1:1235156/1‑101 (MQ=255)
    gcgGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..gc                                                                                                     <  2:325383/101‑1 (MQ=255)
     ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..gcg                                                                                                    <  1:3230763/100‑1 (MQ=255)
     cgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..gcg                                                                                                    >  2:1161170/1‑101 (MQ=255)
      gggggGCTGGACCGAAATATCTGCGTGTAGTCGGGGGGCTGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..gcgc                                                                                                   <  2:3835088/101‑1 (MQ=25)
       ggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCa                                                                                                  >  2:3915629/1‑101 (MQ=255)
         ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCaaa                                                                                                >  2:521171/1‑101 (MQ=255)
          gCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAAc                                                                                               <  2:2212242/101‑1 (MQ=255)
             tGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..a                                                                                          >  1:4838796/1‑101 (MQ=255)
              ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..aa                                                                                         >  2:170361/1‑101 (MQ=255)
              ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..aa                                                                                         >  1:4133036/1‑101 (MQ=255)
               gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..aaa                                                                                        >  1:138248/1‑101 (MQ=255)
                aTCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAg                                                                                       <  2:1310250/101‑1 (MQ=255)
                  cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTc                                                                                     >  2:2767760/1‑101 (MQ=255)
                  cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGCATG.GATCTG..GCGCAAACGG..AAAGTc                                                                                     >  2:291157/1‑101 (MQ=255)
                    aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g                                                                                  <  1:89249/101‑1 (MQ=255)
                    aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g                                                                                  >  1:3585406/1‑101 (MQ=255)
                    aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g                                                                                  <  1:220889/101‑1 (MQ=255)
                     aaCATCTGCATGTGGTGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g..t                                                                               <  1:2106180/101‑1 (MQ=255)
                      aCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g..tc                                                                              >  1:1928921/1‑101 (MQ=255)
                        ctCTGCATGTGGCGGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.g..tcgt                                                                            <  2:1548907/100‑1 (MQ=255)
                          cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..tcgtcg                                                                          >  1:3700067/1‑101 (MQ=255)
                           tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGc                                                                         >  2:2857642/1‑101 (MQ=255)
                              aTGTAGTCGGCGGGAAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTcc                                                                      <  1:3124245/101‑1 (MQ=255)
                               tGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTccc                                                                     >  2:228840/1‑101 (MQ=255)
                                gTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTcccc                                                                    <  1:1161287/101‑1 (MQ=255)
                                  cggcggcggGCAGTTCGACCCGCTGTCGCTGCCCCCGCGGGGGCGTGGATCTAAACT.GGTTTG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTcccccc                                                                  <  1:798152/98‑1 (MQ=38)
                                      ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCgg                                                              >  2:2121547/1‑101 (MQ=255)
                                         gggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCggggg                                                           >  1:4889940/1‑101 (MQ=255)
                                          ggCGGTCCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGc                                                          <  1:3888888/101‑1 (MQ=255)
                                            cAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGt                                                        <  1:2690571/101‑1 (MQ=255)
                                                 cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                   <  2:721739/101‑1 (MQ=255)
                                                  gcgcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                    <  1:98905/93‑1 (MQ=0)
                                                    gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..GGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                   >  2:506593/5‑101 (MQ=2)
                                                    gCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                <  2:4364384/101‑1 (MQ=255)
                                                     ccGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  1:4527538/101‑1 (MQ=255)
                                                     attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                  <  1:60664/98‑1 (MQ=2)
                                                       tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCgca                                                >  2:2672616/2‑99 (MQ=2)
                                                       tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                                >  2:3575117/2‑101 (MQ=2)
                                                       gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                             >  2:4111666/1‑101 (MQ=255)
                                                       ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                               <  2:4923984/98‑1 (MQ=2)
                                                        tggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                              >  2:2308317/3‑101 (MQ=2)
                                                         ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                             >  2:3711699/2‑101 (MQ=2)
                                                         aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                              <  2:4756481/101‑1 (MQ=2)
                                                         aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTg                                           <  2:1488763/101‑1 (MQ=255)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                             >  1:1683291/1‑100 (MQ=2)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                          >  1:1314846/1‑101 (MQ=255)
                                                          gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                          >  2:4007226/1‑101 (MQ=255)
                                                          gTCGCTCCCCCCGCGGGGGAGTGGATCGAAACT.GGTATG.GATCTG..GCGTAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                          <  1:2830236/101‑1 (MQ=255)
                                                           tCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcac                                            <  1:2535894/101‑3 (MQ=2)
                                                           gcGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.ATGACCAAAT..CCGCAATCGG..CCGGGCT.GGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCta                                           <  1:4264405/100‑2 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGTGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                          <  2:4539446/101‑3 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacc                                           <  1:1849072/101‑4 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                          <  1:4949381/101‑3 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtat                                          <  2:4289950/101‑3 (MQ=2)
                                                            cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGGCT.G..TCCTCGCTCCCCCCGCGGGGGCGGGGCTCGAAAGCTgtgg                                        >  2:4184789/1‑101 (MQ=25)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                          >  1:2790962/1‑97 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                          >  2:4984721/1‑97 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                         >  2:3416623/1‑98 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                         <  1:155601/101‑4 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                         <  1:1057341/101‑4 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                         <  1:3994964/101‑4 (MQ=2)
                                                             gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                       >  1:3488700/1‑101 (MQ=255)
                                                              cTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccg                                         <  2:4893434/101‑6 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccgg                                        >  2:4618857/1‑95 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcgtatgag                                       <  1:2066794/101‑8 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                       >  2:3106645/1‑96 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                       <  1:3037247/101‑6 (MQ=2)
                                                               tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT..GACCAGATCCGCAATCGG..CCGGGCTCGA.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                       >  1:2439499/1‑96 (MQ=2)
                                                                gatctcgcgGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                    <  1:5017625/96‑1 (MQ=255)
                                                                ccccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                       >  1:2159312/1‑94 (MQ=0)
                                                                ccccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                       >  2:1887599/1‑94 (MQ=0)
                                                                ccccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                       <  2:4338879/101‑8 (MQ=0)
                                                                ccccccGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                    >  1:3330049/1‑101 (MQ=255)
                                                                 cccccGCGGGGGCGTGGATCGAAACTCGGCAAG.GACATC..GCCGATCTGGGTAAGGCCG.G..TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggcg                                      >  2:3771739/1‑93 (MQ=2)
                                                                     cgcgGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAgatc                               >  2:5017295/1‑98 (MQ=255)
                                                                       cgGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGcc                             <  1:2679592/101‑1 (MQ=255)
                                                                        gggggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCa                            >  1:2319675/1‑101 (MQ=255)
                                                                           ggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTg                         <  2:3359435/101‑1 (MQ=255)
                                                                           ggCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCGGCGTGAGTTGCCAGTg                         >  1:3621578/1‑101 (MQ=255)
                                                                             cGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                       <  1:4326667/101‑1 (MQ=255)
                                                                             cGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                       <  2:2696841/101‑1 (MQ=255)
                                                                             cGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                       <  1:2384191/101‑1 (MQ=255)
                                                                             cGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                       >  2:2051248/1‑101 (MQ=255)
                                                                              gTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                      >  2:1974341/1‑101 (MQ=255)
                                                                               tGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                     >  1:728176/1‑101 (MQ=255)
                                                                                ggATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                    >  1:4601131/1‑101 (MQ=255)
                                                                                 gATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGACAGAc                   <  2:4273559/101‑1 (MQ=255)
                                                                                  aTCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACg                  <  1:444478/101‑1 (MQ=255)
                                                                                   tCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGc                 <  1:4627997/101‑1 (MQ=255)
                                                                                     gAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGCGATAGACGACg               >  2:4598871/1‑101 (MQ=255)
                                                                                        aCT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg            >  1:2971400/1‑101 (MQ=255)
                                                                                          t.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc          <  1:3282583/101‑1 (MQ=255)
                                                                                             gTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGt        <  2:3621307/101‑1 (MQ=255)
                                                                                                 g.gATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGTCGTCGCt    >  2:4501529/1‑101 (MQ=255)
                                                                                                 g.gATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCt    >  1:3807805/1‑101 (MQ=255)
                                                                                                    aTCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcc  <  1:3425773/101‑1 (MQ=255)
                                                                                                    aTCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTcc  >  2:2225371/1‑101 (MQ=255)
                                                                                                      |                                                                                                         
CCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATG.GATCTG..GCGCAAACGG..AAAGTCT.G..TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGCTCC  >  NC_014034/1395641‑1395839

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