Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 585,935 A→G intergenic (‑302/+212) ompT ← / ← envY outer membrane protease VII (outer membrane protein 3b)/DNA‑binding transcriptional activator EnvY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913585,9350AG97.0% 108.7 / ‑3.1 33intergenic (‑302/+212)ompT/envYouter membrane protease VII (outer membrane protein 3b)/DNA‑binding transcriptional activator EnvY
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (11/21);  minor base . (0/1);  total (11/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

TACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGG  >  NC_000913/585806‑586054
                                                                                                                                 |                                                                                                                       
tACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCaa                                                                                                                      >  2:368777/1‑133 (MQ=255)
tACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCaa                                                                                                                      <  1:496518/133‑1 (MQ=255)
 aCTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCaaa                                                                                                                     >  1:474618/1‑133 (MQ=255)
       tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTgag                                                                                                               <  1:18479/133‑1 (MQ=255)
       tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTgag                                                                                                               >  2:443945/1‑133 (MQ=255)
       tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTgag                                                                                                               <  1:47858/133‑1 (MQ=255)
       tgGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCAATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTgag                                                                                                               >  1:441400/1‑133 (MQ=255)
              tttATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATa                                                                                                        <  2:69878/133‑1 (MQ=255)
                        tgtgGTGCTTTTATCTATTTAAGTGCAATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCa                                                                                              <  2:441400/133‑1 (MQ=255)
                                ttttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTg                                                                                      <  2:64924/133‑1 (MQ=255)
                                  ttATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCa                                                                                    <  1:414387/133‑1 (MQ=255)
                                       tATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCt                                                                               >  1:263229/1‑133 (MQ=255)
                                          ttAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACGTAAGGCCGCGCAATAGCTTGCAAGGCTTTa                                                                            >  1:249128/1‑133 (MQ=255)
                                           tAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACa                                                                          <  1:1172/133‑1 (MQ=255)
                                           tAAGTGCCATGCCTTTAGAGGCATATAAGCAAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTAc                                                                           <  2:270940/133‑1 (MQ=255)
                                            aaGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACa                                                                          <  1:44356/133‑1 (MQ=255)
                                            aaGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACa                                                                          <  1:68156/133‑1 (MQ=255)
                                              gTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATg                                                                        <  1:443945/133‑1 (MQ=255)
                                                   aagCCTTTAGAGGCACATAAGCGAAATTTGCACGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTg                                                                   <  2:249128/131‑1 (MQ=255)
                                                   aTGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTg                                                                   <  2:474618/133‑1 (MQ=255)
                                                         ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAg                                                             <  2:406743/133‑1 (MQ=255)
                                                         ttAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAg                                                             <  1:11249/133‑1 (MQ=255)
                                                                 atatAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt                                                     >  1:477/1‑133 (MQ=255)
                                                                 atatAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgt                                                     >  1:332984/1‑133 (MQ=255)
                                                                   ataAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACgtgt                                                   <  2:263229/133‑1 (MQ=255)
                                                                                gCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATg                                      <  1:480725/133‑1 (MQ=255)
                                                                                 cATGAGGTTTATCCTCAATTACTATGTTTTTTACAACAAAAATGAGGGGCAAAACTGAGACACATAAGCCGTTGCAATGGCTTACAAGGCTTAACTTGTTTTAAGGAAGTGGGACGTGTGAGCCCAGAGAttg                                     >  2:83359/1‑131 (MQ=255)
                                                                                  aTGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGgc                                    <  1:368777/133‑1 (MQ=255)
                                                                                                 aaTTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGa                     >  2:383254/1‑133 (MQ=255)
                                                                                                 aaTTACTATGTTTTTTAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGa                     >  2:469815/1‑133 (MQ=255)
                                                                                                      cTATGTTTTTTAGTAAAACAAAGAGGGGCACAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTaa                <  1:200430/133‑1 (MQ=255)
                                                                                                             ttttAGTACAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCg         >  2:327267/1‑133 (MQ=255)
                                                                                                                    aCAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGTGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAgg  <  2:332984/133‑1 (MQ=255)
                                                                                                                    aCAAAAAAGAGGGGCAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAgg  <  1:64155/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                       
TACTGTGTGGTTTTTTTATGCTTTTGTGGTGCTTTTATCTATTTAAGTGCCATGCCTTTAGAGGCATATAAGCGAAAATAGCATGAGGTTTATCCTCAATTACTATGTTTTTTAGTACAAAAAAGAGGGACAAAACTGAGACACATAAGGCCTCGCAATGGCTTGCAAGGCTTTACATGTTTTGAGGTAGTGGGACGTGTGAGCGCAGAGATGGCGCGGTAAGTTGTTGACTTAAAATGTCGTTCTAGG  >  NC_000913/585806‑586054

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