Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,395,736 T→A 15.4% intergenic (+548/‑383) 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,7360T→A15.4% 136.2 ‑0.3 ‑1.5300intergenic (+548/‑383)RCAP_rcc01306/RCAP_rcc01307Cas2 family CRISPR‑associated protein/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (20/26):  ref base (119/134):  total (140/160)
Fisher's exact test strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.28e-02

Only 100 of 308 total aligned reads displayed.

CTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGC  >  NC_014034/1395637‑1395836
                                                                                                    |                                                                                                           
 tCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                           >  2:2805090/1‑101 (MQ=255)
 cccccccGGGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                           <  1:1760071/100‑1 (MQ=255)
 cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                           <  1:2108345/100‑1 (MQ=255)
 cccccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATgg                                                                                                           <  1:2853699/100‑1 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                         >  2:1041889/1‑101 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                         <  1:3621762/101‑1 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                         <  1:299709/101‑1 (MQ=255)
   cccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGAt                                                                                                         <  1:2087672/101‑1 (MQ=255)
    ccccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATc                                                                                                        <  2:1935384/101‑1 (MQ=255)
     cccGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCt                                                                                                       <  1:2805292/101‑1 (MQ=255)
        gcgGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gc                                                                                                 >  1:1235156/1‑101 (MQ=255)
        gcgGGGGCGGGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gc                                                                                                 <  2:325383/101‑1 (MQ=255)
         ggggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gcg                                                                                                <  1:3230763/100‑1 (MQ=255)
          gggggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..gcgc                                                                                               <  1:5029480/101‑1 (MQ=255)
            gggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaa                                                                                             <  2:351360/101‑1 (MQ=255)
             ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaaa                                                                                            >  1:728828/1‑101 (MQ=255)
             ggCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCaaa                                                                                            >  2:521171/1‑101 (MQ=255)
              gCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAAc                                                                                           <  2:2212242/101‑1 (MQ=255)
                gTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGAGCGAAACT.GGTATGGATCT.G..GCGCAAACgg                                                                                         >  2:479927/1‑101 (MQ=255)
                  ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aa                                                                                     >  2:170361/1‑101 (MQ=255)
                  ggATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aa                                                                                     <  1:4171172/101‑1 (MQ=255)
                   gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                    >  1:138248/1‑101 (MQ=255)
                   gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                    >  2:3603901/1‑101 (MQ=255)
                   gATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..aaa                                                                                    >  2:3665253/1‑101 (MQ=255)
                      cGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGCATGGATCT.G..GCGCAAACGG..AAAGTc                                                                                 >  2:291157/1‑101 (MQ=255)
                        aaaCATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g                                                                              >  1:2734176/1‑101 (MQ=255)
                           cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcg                                                                          >  2:4155464/1‑101 (MQ=255)
                           cATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcg                                                                          <  1:4909450/101‑1 (MQ=255)
                             tCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.g.tcgtc                                                                        >  2:2678765/1‑101 (MQ=255)
                              cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                       >  1:2398007/1‑101 (MQ=255)
                              cTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.tcgtcg                                                                       >  1:3700067/1‑101 (MQ=255)
                               tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                      >  2:2857642/1‑101 (MQ=255)
                               tGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGc                                                                      >  2:2260291/1‑101 (MQ=255)
                                  aTGTAGTCGGCGGGAAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcc                                                                   <  1:3124245/101‑1 (MQ=255)
                                   ggTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AACGTCT.G.TCGTCGCTccc                                                                  <  1:4413511/100‑1 (MQ=255)
                                    gTAGTCGGGGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                 <  1:1161287/101‑1 (MQ=255)
                                    gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                 <  2:4599327/101‑1 (MQ=255)
                                    gTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                 >  1:4734729/1‑101 (MQ=255)
                                    gTAGTCGGCGGGAAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTcccc                                                                 <  2:4215985/101‑1 (MQ=255)
                                        tgggggggCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCcgc                                                             <  2:4502216/96‑1 (MQ=255)
                                          ggcggGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                           >  2:2121547/1‑101 (MQ=255)
                                          ggcggGCAGTGCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCgg                                                           <  2:3864247/101‑1 (MQ=255)
                                                 aGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCNCCCCCCGCGGGGGCGTg                                                    >  2:1196052/1‑101 (MQ=255)
                                                  gTTGGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTgg                                                   <  1:3915921/101‑1 (MQ=255)
                                                    tCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 >  1:2637204/1‑101 (MQ=255)
                                                    tCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 <  1:3604172/101‑1 (MQ=255)
                                                     cGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATc                                                <  2:721739/101‑1 (MQ=255)
                                                      gAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               <  1:2328698/101‑1 (MQ=255)
                                                        gcgctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               <  1:851139/95‑1 (MQ=0)
                                                        gattcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.GGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                >  2:506593/5‑101 (MQ=2)
                                                        atctcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATc                                                <  1:4612303/97‑1 (MQ=2)
                                                         attcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACATCG..GCGATCTGGGT.AAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCg                                               >  1:2002408/4‑101 (MQ=2)
                                                          ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              <  1:2681320/99‑1 (MQ=2)
                                                          ttcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              <  2:4752545/99‑1 (MQ=2)
                                                           tcAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGaa                                             >  2:3778276/2‑101 (MQ=2)
                                                           gCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                          >  2:273323/1‑101 (MQ=255)
                                                           ctggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGaaa                                            <  2:4923984/98‑1 (MQ=2)
                                                            tggTCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                           <  2:1387470/99‑1 (MQ=2)
                                                             ggTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcc                                          >  2:2289756/2‑101 (MQ=2)
                                                             aGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAc                                           <  1:955020/101‑1 (MQ=2)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAca                                          >  1:1683291/1‑100 (MQ=2)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCt                                         <  2:4361923/101‑1 (MQ=2)
                                                              gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgc                                       >  1:1314846/1‑101 (MQ=255)
                                                              gTCGCCCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAGAca                                          <  1:1896147/101‑2 (MQ=2)
                                                               tCGCCCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcg                                      <  1:4286753/101‑1 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                     >  2:4903140/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                     >  1:452418/1‑101 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcgg                                     <  1:4564855/101‑1 (MQ=255)
                                                                cGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCCgcgg                                     >  1:2222896/1‑101 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcaccc                                       >  2:4984721/1‑97 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                      <  1:155601/101‑4 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                      >  1:4871950/1‑98 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatg                                      <  1:1057341/101‑4 (MQ=2)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggt                                    <  2:4565843/101‑2 (MQ=255)
                                                                 gCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggc                                    <  1:618059/101‑1 (MQ=255)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACT.TG.AT.GACCA.GATCCGCAATCGG..CCGGGCTCGATCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCtatgag                                    >  2:3106645/1‑96 (MQ=2)
                                                                   tCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTgcggcgg                                  >  2:1599849/1‑100 (MQ=255)
                                                                    ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                    >  2:1887599/1‑94 (MQ=0)
                                                                    ccccccGCGGGGGCGTGGATCGAAACTCGGCAAGGACAT.C..GCCGATCTGGGTAAGGCCG.G.TGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAcacccggc                                    >  1:2159312/1‑94 (MQ=0)
                                                                    ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                 <  2:1530746/101‑1 (MQ=255)
                                                                    ccccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTg                                 >  1:3330049/1‑101 (MQ=255)
                                                                     cccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGa                                <  2:3727125/101‑1 (MQ=255)
                                                                      ccccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAg                               <  1:3810210/101‑1 (MQ=255)
                                                                       cccGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTTAGt                              <  2:4274163/101‑1 (MQ=255)
                                                                         cgcgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAgatc                            >  2:5017295/1‑98 (MQ=255)
                                                                         cgcgGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTg                            <  1:3435277/101‑1 (MQ=255)
                                                                               ggCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTg                      <  2:3359435/101‑1 (MQ=255)
                                                                                 cGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGAt                    >  2:2051248/1‑101 (MQ=255)
                                                                                  gTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATa                   >  2:1974341/1‑101 (MQ=255)
                                                                                   tGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTTCGGCGTGAGTTGCCAGTGATAg                  <  2:3948830/101‑1 (MQ=255)
                                                                                   tGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAg                  <  2:4009932/101‑1 (MQ=255)
                                                                                   gggATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCT.CCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGa                 <  1:3491890/100‑1 (MQ=255)
                                                                                      aTCGAAAAT.GGAATGGATCA.G..GCGCAAACGG..TAAGTCT.G.TCGTCGGTCCCCCCGCGGGGGCGTGGATCGAAGCACGCGACGGGAGGTCCCAGCGATAGagg               >  2:2443276/1‑99 (MQ=16)
                                                                                          aaaCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCGGTGATAGACGCCGt           >  1:4973808/1‑101 (MQ=255)
                                                                                            aCT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGCGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTgg         >  1:2971400/1‑101 (MQ=255)
                                                                                             cT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGa        >  1:4036796/1‑101 (MQ=255)
                                                                                              t.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGAc       <  1:3282583/101‑1 (MQ=255)
                                                                                                  tATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc    <  2:4305053/101‑1 (MQ=255)
                                                                                                  tATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTc    >  1:2558902/1‑101 (MQ=255)
                                                                                                   aTGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACgctg   >  1:2485892/1‑98 (MQ=255)
                                                                                                    |                                                                                                           
CTCCCCCCGCGGGGGCGTGGATCGAAACATCTGCATGTAGTCGGCGGGCAGTTCGAGCCGCAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACT.GGTATGGATCT.G..GCGCAAACGG..AAAGTCT.G.TCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTGCGGCGTGAGTTGCCAGTGATAGACGCCGTGGACGTCGC  >  NC_014034/1395637‑1395836

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