Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,661,263 Δ1 bp 7.3% coding (708/2424 nt) ynfF → S‑ and N‑oxide reductase, A subunit, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,661,2630G. 7.3% 184.9 / 7.1 55coding (708/2424 nt)ynfFS‑ and N‑oxide reductase, A subunit, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (20/31);  new base . (2/2);  total (22/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TCAGTCGCTATGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCACCGATGGCGCA  >  NC_000913/1661118‑1661407
                                                                                                                                                 |                                                                                                                                                
tcAGTCGCTATGGCCGCTACAGTACCGCCCAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAAGACGAAACTGGTCGGTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAgc                                                                                                                                               >  5:19168/1‑149 (MQ=255)
  aGTCGCTATGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAgcg                                                                                                                                              <  5:228983/148‑1 (MQ=255)
    tCGCTATGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAgcggcg                                                                                                                                           <  6:142761/149‑1 (MQ=255)
       cTATGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCggtg                                                                                                                                        <  4:283851/149‑1 (MQ=255)
          tGGCAGCTACAGTACCGCACAGTTCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCGGTATGATATCTCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCggtggt                                                                                                                                      <  6:148862/148‑1 (MQ=255)
          tGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGcggcgg                                                                                                                                          >  1:302749/1‑144 (MQ=255)
          tGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGcggcgg                                                                                                                                          <  2:302749/144‑1 (MQ=255)
             cAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCa                                                                                                                                  <  3:362577/149‑1 (MQ=255)
              aGCTACAGTACCACACAGATCAGTGCCGCAATGAGTTATACGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAgcgg                                                                                                                                             <  1:88313/137‑1 (MQ=255)
              aGCTACAGTACCACACAGATCAGTGCCGCAATGAGTTATACGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAgcgg                                                                                                                                             >  2:88313/1‑137 (MQ=255)
                    aGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAATTAACCCGTCAGAAACGCGGATGAGCGGCGGTGGTGTCACTtacta                                                                                                                           <  1:10316/149‑1 (MQ=255)
                        ccGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTc                                                                                                                       <  2:2644/149‑1 (MQ=255)
                              aGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCaa                                                                                                                 <  5:290455/149‑1 (MQ=255)
                                       ccGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATTATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGaa                                                                                                        <  1:13805/149‑1 (MQ=255)
                                        cGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAAc                                                                                                       <  2:238739/149‑1 (MQ=255)
                                        cGCAATCAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAAc                                                                                                       <  2:24407/149‑1 (MQ=255)
                                          cAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAACACGAAACTGGTCGGTATGTTCGGAAATAACCCGGCAGCAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCACGCCCGCGAACGt                                                                                                     >  8:188505/1‑149 (MQ=255)
                                           aaTGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGtt                                                                                                    >  6:71739/1‑149 (MQ=255)
                                             tGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCACTACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCa                                                                                                  <  4:33274/149‑1 (MQ=255)
                                                                cAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCa                                                                               <  6:324472/149‑1 (MQ=255)
                                                                 aaTGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCAt                                                                              <  4:190217/149‑1 (MQ=255)
                                                                 aaTGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCAt                                                                              >  8:325580/1‑149 (MQ=255)
                                                                  atgatgGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATc                                                                             >  6:5168/1‑149 (MQ=255)
                                                                   tgatgGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGAT‑AGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGa                                                                           >  4:331904/1‑149 (MQ=255)
                                                                   tgatgGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGAT‑AGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGa                                                                           >  4:337255/1‑149 (MQ=255)
                                                                     atgGCAACAGCCCGGATGATAACGCCAATCCGTAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGAt                                                                          <  2:140152/149‑1 (MQ=255)
                                                                        gCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATcc                                                                        <  7:325580/148‑1 (MQ=255)
                                                                         cAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCAc                                                                      <  7:266272/149‑1 (MQ=255)
                                                                          aaCAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACg                                                                     >  1:153616/1‑149 (MQ=255)
                                                                          aaCAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACg                                                                     >  6:168731/1‑149 (MQ=255)
                                                                               cccGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTtata                                                                <  6:66025/149‑1 (MQ=255)
                                                                                ccGGATGATATCGACAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAATCGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATaa                                                               <  3:220618/149‑1 (MQ=255)
                                                                                  ggATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAAc                                                              >  6:279471/1‑148 (MQ=255)
                                                                                   gatgatATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTTGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCTTCATCGATCCACGTTATAACGa                                                            <  5:110609/149‑1 (MQ=255)
                                                                                   gatgatATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATc                                                                                   >  5:9883/1‑126 (MQ=255)
                                                                                   gatgatATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATc                                                                                   <  6:9883/126‑1 (MQ=255)
                                                                                        tatCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGACCGTTCAAACGCACGCATGATCGTCATCGAGCCACGTTATAACGACActg                                                       >  6:284918/1‑149 (MQ=255)
                                                                                          tCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACActgct                                                     >  7:91512/1‑149 (MQ=255)
                                                                                             ccAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTta                                                                  >  2:350632/1‑133 (MQ=255)
                                                                                             ccAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTta                                                                  <  1:350632/133‑1 (MQ=255)
                                                                                               aaTACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCgg                                                <  7:80591/149‑1 (MQ=255)
                                                                                                    gAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGt                                            >  7:78907/1‑148 (MQ=255)
                                                                                                     aaaCTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGa                                          <  5:70376/149‑1 (MQ=255)
                                                                                                     aaaCTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGa                                          <  7:361083/149‑1 (MQ=255)
                                                                                                       aCTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGAT‑AGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGa                                       <  3:337255/149‑1 (MQ=255)
                                                                                                       aCTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGAT‑AGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGa                                       <  3:331904/149‑1 (MQ=255)
                                                                                                       aCTGGTCGTTATGTTCGGAAAAAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCCAACGTTCCAACTCACGCATGATCGTCATCCAACCACGTTATACCGACACTGCTGCCGGGCgag                                           >  8:150633/1‑144 (MQ=255)
                                                                                                           gTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCTTCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGAt                                      <  5:284918/147‑1 (MQ=255)
                                                                                                                         aaaTAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTc                      <  5:71739/149‑1 (MQ=255)
                                                                                                                              aCCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCt                 <  5:5168/149‑1 (MQ=255)
                                                                                                                              aCCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCt                 <  8:78907/149‑1 (MQ=255)
                                                                                                                                  ggCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCa             >  7:40825/1‑149 (MQ=255)
                                                                                                                                   gCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCAc            >  4:157002/1‑149 (MQ=255)
                                                                                                                                           gcgGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCACCGATGgcg    >  4:270959/1‑149 (MQ=255)
                                                                                                                                             ggATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGCTCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCACCGATGGCGCa  >  8:138641/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCAGTCGCTATGGCAGCTACAGTACCGCACAGATCAGTGCCGCAATGAGTTATATGTTCGGTGCCAATGATGGCAACAGCCCGGATGATATCGCCAATACGAAACTGGTCGTTATGTTCGGAAATAACCCGGCAGAAACGCGGATGAGCGGCGGTGGTGTCACTTACTACGTCGAGCAAGCCCGCGAACGTTCAAACGCACGCATGATCGTCATCGATCCACGTTATAACGACACTGCTGCCGGGCGTGAAGATGAATGGCTGCCCATTCGCCCTGGCACCGATGGCGCA  >  NC_000913/1661118‑1661407

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: ‑