Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,064,023 T→A 7.5% C119S (TGC→AGC)  ubiA → 4‑hydroxybenzoate octaprenyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,064,0230TA 7.5% 187.6 / 5.5 67C119S (TGC→AGC) ubiA4‑hydroxybenzoate octaprenyltransferase
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (2/3);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

TGCTGCATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCTTTCACAGCCGCT  >  NC_002947/6063876‑6064169
                                                                                                                                                   |                                                                                                                                                  
tgctgcATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCtgtg                                                                                                                                                   >  6:215857/1‑149 (MQ=255)
  ctgcATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGc                                                                                                                                                  <  3:180150/148‑1 (MQ=255)
  ctgcATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCa                                                                                                                                                 <  2:385739/149‑1 (MQ=255)
   tgcATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCATGCTGGTGCTGTGCAc                                                                                                                                                >  8:82001/1‑149 (MQ=255)
         aaCGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGTTGGTGCACAGCACCAGCa                                                                                                                                           <  4:208276/148‑1 (MQ=255)
         aaCGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCACAGCACCAGCa                                                                                                                                           <  3:316767/148‑1 (MQ=255)
          aCGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGc                                                                                                                                         >  6:273678/1‑149 (MQ=255)
            gACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCACAGCACCAGCa                                                                                                                                           <  3:88504/145‑1 (MQ=255)
            gACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCACAGCACCAGCa                                                                                                                                           >  4:88504/1‑145 (MQ=255)
                    cGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCt                                                                                                                               >  7:177882/1‑149 (MQ=255)
                      aCCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTgg                                                                                                                             <  3:137753/149‑1 (MQ=255)
                          cAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGt                                                                                                                         >  7:350133/1‑149 (MQ=255)
                          cAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGt                                                                                                                         <  5:59598/149‑1 (MQ=255)
                          cAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGt                                                                                                                         >  2:318921/1‑149 (MQ=255)
                           aaGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTc                                                                                                                        <  5:241905/149‑1 (MQ=255)
                             ggTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGt                                                                                                                      >  3:380103/1‑149 (MQ=255)
                             ggTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCAGACAGTGCGGCTGttgg                                                                                                                      <  1:247512/149‑4 (MQ=255)
                                    ggCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGc                                                                                                               >  4:281419/1‑149 (MQ=255)
                                    ggCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGc                                                                                                               >  5:54915/1‑149 (MQ=255)
                                    ggCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGc                                                                                                               >  6:139860/1‑149 (MQ=255)
                                    ggCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGc                                                                                                               <  1:237928/149‑1 (MQ=255)
                                      ccACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAg                                                                                                                                          <  2:77447/120‑1 (MQ=255)
                                      ccACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAg                                                                                                                                          >  1:77447/1‑120 (MQ=255)
                                        aCGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTg                                                                                                           <  6:40278/149‑1 (MQ=255)
                                              aGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTTCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGGCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTg                                                                                                     <  8:333817/149‑1 (MQ=255)
                                                cgcACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGc                                                                                                   <  4:173461/149‑1 (MQ=255)
                                                    cagcagACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTc                                                                                               >  1:171358/1‑149 (MQ=255)
                                                      gcagACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTg                                                                                             <  5:139860/149‑1 (MQ=255)
                                                      gcagACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTCTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTg                                                                                             <  1:318921/149‑1 (MQ=255)
                                                       cagACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCTCCAACAGCCGCTCTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGGGTTCTGc                                                                                            >  4:358658/1‑149 (MQ=255)
                                                          aCCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTAc                                                                                         >  5:313915/1‑149 (MQ=255)
                                                             gCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCg                                                                                      <  7:90567/149‑1 (MQ=255)
                                                               cccTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGtt                                                                                    <  5:64181/149‑1 (MQ=255)
                                                                 cTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCa                                                                                  <  6:304743/149‑1 (MQ=255)
                                                                 cTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCa                                                                                  <  8:177882/149‑1 (MQ=255)
                                                                 cTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCa                                                                                  >  1:56102/1‑149 (MQ=255)
                                                                    ggCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCAt                                                                                 <  3:281419/147‑1 (MQ=255)
                                                                    ggCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGGGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCAt                                                                                 <  8:219120/147‑1 (MQ=255)
                                                                        aGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTAccc                                                                                       >  8:281664/1‑137 (MQ=255)
                                                                        aGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTAccc                                                                                       <  7:281664/137‑1 (MQ=255)
                                                                        aGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAgcg                                                                           <  4:380103/149‑1 (MQ=255)
                                                                                ggTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACct                                                                   <  5:273678/149‑1 (MQ=255)
                                                                                     ggCCACGCGAGGTGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACctact                                                                >  2:379978/1‑147 (MQ=255)
                                                                                     ggCCACGCGAGGTGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACctact                                                                <  1:379978/147‑1 (MQ=255)
                                                                                     ggCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCtactac                                                              <  2:171358/149‑1 (MQ=255)
                                                                                     ggCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCtactac                                                              <  1:57795/149‑1 (MQ=255)
                                                                                       gcACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTAccc                                                            >  6:283187/2‑149 (MQ=255)
                                                                                              aGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAggtg                                                     <  6:1944/149‑1 (MQ=255)
                                                                                                  gcTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGc                                                 >  5:99931/1‑149 (MQ=255)
                                                                                                    tGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAg                                                                             >  8:43930/1‑119 (MQ=255)
                                                                                                    tGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAg                                                                             <  7:43930/119‑1 (MQ=255)
                                                                                                      acacTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTg                                               <  6:99931/147‑1 (MQ=255)
                                                                                                            ttCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGggcggc                                       >  7:333417/1‑149 (MQ=255)
                                                                                                                 gATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTAtt                                  >  7:16090/1‑149 (MQ=255)
                                                                                                                    tCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCg                                                                                                  >  4:78772/1‑82 (MQ=255)
                                                                                                                    tCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCg                                                                                                  <  3:78772/82‑1 (MQ=255)
                                                                                                                             ggTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTc                      >  7:308544/1‑149 (MQ=255)
                                                                                                                                 aGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCAt                  <  6:54915/149‑1 (MQ=255)
                                                                                                                                  tgcggctgttgGTGCACAGCACCAGCAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATg                 >  7:353451/6‑149 (MQ=255)
                                                                                                                                   tttCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATgg                >  7:300119/1‑149 (MQ=255)
                                                                                                                                   tttCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATgg                >  3:362465/1‑149 (MQ=255)
                                                                                                                                   tttCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATg                 <  6:313915/148‑1 (MQ=255)
                                                                                                                                     tCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGc               <  3:373844/148‑1 (MQ=255)
                                                                                                                                       ctgctgGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCttt            >  3:266732/1‑149 (MQ=255)
                                                                                                                                       ctgctgGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCttt            <  5:283187/149‑1 (MQ=255)
                                                                                                                                       ctgctgGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCtt             <  8:350133/148‑1 (MQ=255)
                                                                                                                                          ctgGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCTTTca          <  7:106720/148‑1 (MQ=255)
                                                                                                                                           tgGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCTTTcc          >  7:352264/1‑146 (MQ=255)
                                                                                                                                                 tgtgCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCTTTCACAGCCGCt  >  4:177078/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TGCTGCATCAACGACTTTGCCGACCGCAAGGTTGACGGCCACGTCAAGCGCACAGCAGACCGCCCCTTGGCCAGCGGCCGGGTGCGGCCACGCGAGGCGCTGACACTGTTCGCGATTCTGGTCGGGGTGAGTTTCCTGCTGGTGCTGTGCACCAACAGCCGCACTGTCTGGCTGTCGTTCGGGGCGGTGGCCTTGGCGTTCTGCTACCCGTTCATGAAGCGCTACACCTACTACCCGCAGGTGGTGCTGGGGGCGGCGTATTCGTGGGGCATTCCCATGGCTTTCACAGCCGCT  >  NC_002947/6063876‑6064169

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: