Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,400,872 C→A 9.0% intergenic (+35/‑134) ybgF → / → queE TolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,400,8720CA 9.0% 260.5 / 14.2 89intergenic (+35/‑134)ybgF/queETolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base C (37/44);  new base A (4/4);  total (41/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.42e-02

CCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACA  >  NC_002947/1400729‑1401015
                                                                                                                                                |                                                                                                                                                
ccGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTc                                                                                                                                             >  7:165144/1‑149 (MQ=255)
ccGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTc                                                                                                                                             <  1:221109/149‑1 (MQ=255)
   tATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGTCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTcgcg                                                                                                                                          >  7:405713/1‑149 (MQ=255)
     tGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCgg                                                                                                                                         <  1:11754/148‑1 (MQ=255)
     tGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCgg                                                                                                                                         <  8:392645/148‑1 (MQ=255)
       ggCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtt                                                                                                                                      >  5:18340/1‑149 (MQ=255)
        gCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGttt                                                                                                                                     >  7:277150/1‑149 (MQ=255)
        gCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGttt                                                                                                                                     <  1:6553/149‑1 (MQ=255)
          cATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGACTGGCGCGGGtttct                                                                                                                                   <  8:206579/149‑3 (MQ=255)
            tACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTc                                                                                                                                >  1:128822/1‑149 (MQ=255)
                 aCAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGTTTCTTTCgag                                                                                                                             >  1:251089/1‑145 (MQ=38)
                 aCAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAAAACCCGCGACTGGCGCGGGTTTCTTTCgag                                                                                                                             <  2:251089/147‑3 (MQ=38)
                    agctAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAAt                                                                                                                        <  8:281329/146‑1 (MQ=255)
                    aGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAAt                                                                                                                        <  4:245864/149‑1 (MQ=255)
                     ggTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAAt                                                                                                                        >  7:6395/1‑148 (MQ=255)
                     ggTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGACAGTCGCGGGTTTTTTCGTTAGAATc                                                                                                                       <  1:234336/149‑1 (MQ=255)
                         aaaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTg                                                                                                                   >  4:247473/1‑149 (MQ=255)
                          aaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGc                                                                                                                  <  5:42729/149‑1 (MQ=255)
                          aaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGc                                                                                                                  >  3:304077/1‑149 (MQ=255)
                                  cTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCtttttt                                                                                                          >  3:74133/1‑149 (MQ=255)
                                    tCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCg                                                                                                        >  1:115870/1‑149 (MQ=255)
                                     cagcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCTCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGTGTTTTTTCGTTAGAATCACTGCCCTTTTTTc                                                                                                         <  4:78592/147‑1 (MQ=255)
                                      agcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTa                                                                                                      <  4:420259/149‑1 (MQ=255)
                                      agcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTa                                                                                                      <  4:413691/149‑1 (MQ=255)
                                          ggTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAaca                                                                                                  >  4:292838/1‑149 (MQ=255)
                                           gTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAacac                                                                                                 >  3:280407/1‑149 (MQ=255)
                                            tCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACg                                                                                                <  8:19201/149‑1 (MQ=255)
                                                gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCGAGGCGTTGTCGCGGTCCCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTGGTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgctgc                                                                                            >  5:140919/1‑149 (MQ=255)
                                                 aCCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgctgct                                                                                           >  3:150479/1‑149 (MQ=255)
                                                          cccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAAc                                                                                  >  5:316863/1‑149 (MQ=255)
                                                           gcGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACg                                                                                 <  2:128822/148‑1 (MQ=255)
                                                           ccGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACg                                                                                 <  3:292838/149‑1 (MQ=255)
                                                            gggCAACTCTGCCGCTCAACTGGCCCAGCGTGACGTTCTGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTc                                                                                                                                <  2:265793/100‑1 (MQ=255)
                                                            cGGCACCTCTGCCGCTCAACTGGCCCAGCGTTTACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtttttta                                                                                                                                >  1:265793/1‑101 (MQ=255)
                                                            cGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGc                                                                                >  1:345522/1‑149 (MQ=255)
                                                                aCCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATg                                                                            >  1:73434/1‑149 (MQ=255)
                                                                 ccTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGt                                                                           <  8:6395/149‑1 (MQ=255)
                                                                       ccGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGa                                                                        <  3:247473/146‑1 (MQ=255)
                                                                         gCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcg                                                                   <  2:253173/149‑1 (MQ=255)
                                                                          cTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                  >  4:267247/1‑149 (MQ=255)
                                                                          cTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                  <  6:140919/149‑1 (MQ=255)
                                                                             aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGa                                                                >  6:423034/1‑148 (MQ=255)
                                                                               cTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAg                                                             >  3:141505/1‑149 (MQ=255)
                                                                                  gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                          <  4:304077/149‑1 (MQ=255)
                                                                                   cccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGATGCGGATAACGCTATGTCGAGTCCGCGACAGTGcc                                                         >  8:289188/1‑149 (MQ=255)
                                                                                    ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCACTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                        >  6:249797/1‑149 (MQ=255)
                                                                                           tGACCTGCAGAAGGTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCTGTCGCGGGTTTTTTCGTTAGAATCACTGCGCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTGCGCGAGAGTGCCTGACGGa                                                  <  4:172698/149‑1 (MQ=255)
                                                                                             aCCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAg                                                 <  6:18340/148‑1 (MQ=255)
                                                                                               cTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCg                                              <  4:150479/149‑1 (MQ=255)
                                                                                               cTGCAGAAGCTCTAGGCGTTGTCGCTGAACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCg                                              <  8:165144/149‑1 (MQ=255)
                                                                                                tGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                  >  5:273727/1‑128 (MQ=255)
                                                                                                tGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                  <  6:273727/128‑1 (MQ=255)
                                                                                                 gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGCGGCGGa                                            >  6:53656/1‑149 (MQ=255)
                                                                                                 gCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGa                                            >  5:287097/1‑149 (MQ=255)
                                                                                                   aGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTTTAATAACTGCCCTTTTTTAGTAAACAAGCAGCTGCCGATAACGCTATGTAGAGTACGCGACACTGCCTGACGGAGGCGGa                                            >  2:414235/1‑147 (MQ=255)
                                                                                                   aGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTACGCGACAGTGCCTGACGGAGGCGGa                                            >  5:271113/1‑147 (MQ=255)
                                                                                                     aaGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAAGGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGc                                        >  8:225939/1‑149 (MQ=255)
                                                                                                        ctctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTg                                     <  2:414639/149‑1 (MQ=255)
                                                                                                           tAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                   >  1:154009/1‑148 (MQ=255)
                                                                                                           tAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGt                                    <  3:330806/147‑1 (MQ=255)
                                                                                                            aGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTGCGCGACAGTGCCTGACGGAGGCGGACAGACTGtt                                   >  7:334691/1‑147 (MQ=255)
                                                                                                             ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                   <  2:141396/146‑1 (MQ=255)
                                                                                                             ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                   <  4:74133/146‑1 (MQ=255)
                                                                                                             ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAg                                <  8:405713/149‑1 (MQ=255)
                                                                                                             ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAg                                >  7:122420/1‑149 (MQ=255)
                                                                                                              gCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGc                               <  4:280407/149‑1 (MQ=255)
                                                                                                               cGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCt                              <  4:141505/149‑1 (MQ=255)
                                                                                                                        tGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGccc                     <  7:65403/149‑1 (MQ=255)
                                                                                                                          tACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGccc                     <  5:423034/147‑1 (MQ=255)
                                                                                                                             cGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                <  6:316863/149‑1 (MQ=255)
                                                                                                                             cGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                >  7:39688/1‑149 (MQ=255)
                                                                                                                              gCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCg               <  8:277150/149‑1 (MQ=255)
                                                                                                                               cTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGa              <  3:267247/149‑1 (MQ=255)
                                                                                                                                 acgaaaaaaCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAg             >  6:45737/7‑148 (MQ=255)
                                                                                                                                  gaaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCa           <  1:241594/149‑1 (MQ=255)
                                                                                                                                   aaagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCAt          <  2:115870/149‑1 (MQ=255)
                                                                                                                                    aagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATg         >  7:349137/1‑149 (MQ=255)
                                                                                                                                    aagaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAAAGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACACCATGTTTAGATGTACCTCCCCTCGAGAGCAt          >  6:300631/1‑148 (MQ=255)
                                                                                                                                     agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                          >  2:365094/1‑131 (MQ=255)
                                                                                                                                     agaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                          <  1:365094/131‑1 (MQ=255)
                                                                                                                                      gaaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCa       <  7:125688/149‑1 (MQ=255)
                                                                                                                                       aaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa      <  2:73434/149‑1 (MQ=255)
                                                                                                                                       aaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa      <  8:349137/149‑1 (MQ=255)
                                                                                                                                       aaaCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCaa      <  2:318895/149‑1 (MQ=255)
                                                                                                                                         aCCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGa    <  1:414235/149‑1 (MQ=255)
                                                                                                                                         aCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                >  1:201730/1‑137 (MQ=255)
                                                                                                                                         aCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGc                <  2:201730/137‑1 (MQ=255)
                                                                                                                                         aCCCGCGACTGGCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAg     >  8:279400/1‑148 (MQ=255)
                                                                                                                                            cgcgCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGaca  >  8:113704/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTTTTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACA  >  NC_002947/1400729‑1401015

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