Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 250,396 A→G 14.8% pseudogene (325/756 nt) lafU → flagellar system protein, promoterless fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913250,3960AG14.8% 69.4 / 5.1 27pseudogene (325/756 nt)lafUflagellar system protein, promoterless fragment
Reads supporting (aligned to +/- strand):  ref base A (13/10);  new base G (2/2);  total (15/12)
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.87e-01

CGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGAT  >  NC_000913/250145‑250659
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                       
tgCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCGGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGGGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAgc                                                                                                                                                                                                                                                     <  1:481057/271‑1 (MQ=255)
     taaaaaaacccaggCCCGCCGCTTCCGACAAACCCTGCCACCGGGGGGGGGCCGGAAGGGACGGGGAAAAAATCGCGGGATGGGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGa                                                                                                                                                                                                                                                 <  2:466818/258‑1 (MQ=255)
                                                                                 gTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAattattatt                                                                                                                                                                    >  1:293691/1‑272 (MQ=255)
                                                                                                  aCGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCGATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAAc                                                                                                                                                                                                                           >  1:339384/1‑200 (MQ=255)
                                                                                                  aCGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGActcgctcg                                                                                                                                                                                     >  1:162103/1‑238 (MQ=255)
                                                                                                   cGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCGATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACg                                                                                                                                                                                                                          <  2:339384/200‑1 (MQ=255)
                                                                                                   cGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGa                                                                                                                                                                                    <  2:162103/238‑1 (MQ=255)
                                                                                                    ggCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAAcg                                                                                                                                                                                                                                                            >  2:248937/1‑165 (MQ=255)
                                                                                                    ggCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACcgatgcga                                                                                                                                                 >  1:230958/1‑272 (MQ=255)
                                                                                                     gcgcTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAAcgc                                                                                                                                                                                                                                                           <  1:248937/165‑1 (MQ=255)
                                                                                                                                                          cAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAg                                                                                                                                                                                                                                                      >  2:90517/1‑117 (MQ=255)
                                                                                                                                                           aaCACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAgc                                                                                                                                                                                                                                                     <  1:90517/117‑1 (MQ=255)
                                                                                                                                                                      ccGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCtggt                                                                                                                                                                                                                  >  1:58997/1‑141 (MQ=255)
                                                                                                                                                                       cGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCtggtg                                                                                                                                                                                                                 <  2:58997/141‑1 (MQ=255)
                                                                                                                                                                          gATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGGACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATCCcga                                                                                                                                                      >  1:447739/1‑197 (MQ=255)
                                                                                                                                                                                            aaTCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCGATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCttt                                                                                                                                                                                                                                >  2:354390/1‑105 (MQ=255)
                                                                                                                                                                                             ttttggggtgggagatTGTTCCCCAGGGTTTACGCGTGCTGTTTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATg                                                        <  2:230958/262‑1 (MQ=255)
                                                                                                                                                                                             aTCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCGATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTa                                                                                                                                                                                                                               <  1:354390/105‑1 (MQ=255)
                                                                                                                                                                                                        gagatTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATaa                                             >  1:291203/1‑272 (MQ=255)
                                                                                                                                                                                                          gatTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAattatta                                                                                                                                                                      >  1:430542/1‑149 (MQ=255)
                                                                                                                                                                                                           atTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAattattat                                                                                                                                                                     <  2:430542/149‑1 (MQ=255)
                                                                                                                                                                                                                  cccAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTAcaacaa                                                                                                                     >  2:329314/1‑190 (MQ=255)
                                                                                                                                                                                                                   ccAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACaacaac                                                                                                                    <  1:329314/190‑1 (MQ=255)
                                                                                                                                                                                                                           tACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTg                                                                                                                  <  1:256034/184‑1 (MQ=255)
                                                                                                                                                                                                                                    tGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAAt                                                         >  1:383337/1‑232 (MQ=255)
                                                                                                                                                                                                                                                   aGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAACACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGAt  >  1:374029/1‑272 (MQ=255)
                                                                                                                                                                                                                                                   aGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGAt  >  1:330019/1‑272 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                       
CGCCGTTAAATAAAATCAGCCCGTCTCATCCGCCAAAACCTGCCACCGTGGCGGTGCCGGAAGAGACGGAGAAAAAAGCGCGTGATGTGAATGAGAAAACGGCGCTGCTGAAGAAGAAATCCGCCACCGAGCTTGGTGAGCTGGCAACCAGTATCAACACCATTGCCCGTGATGCGCATATGGAAGCGAATCTGGAGATGGAGATTGTTCCCCAGGGATTACGCGTGCTGATTAAAGACGACCAGAACCGCAATATGTTTGAACGCGGCAGCGCGAAGATTATGCCGTTCTTTAAAACGCTGCTGGTGGAGCTCGCGCCAGTGTTCGACTCGCTCGATAATAAAATTATTATTACCGGGCATACCGATGCGATGGCCTACAAAAACAATATCTACAACAACTGGAACCTTTCGGGTGACCGCGCGCTTTCGGCTCGTCGGGTGCTGGAAGAGGCCGGAATGCCGGAAGATAAAGTGATGCAGGTAAGCGCAATGGCGGACCAGATGCTGCTGGAT  >  NC_000913/250145‑250659

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