Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,438,548 T→C V200A (GTG→GCG)  flk → putative flagella assembly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,438,5480TC100.0% 111.5 / NA 32V200A (GTG→GCG) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/17);  total (15/17)

GCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACA  >  NC_000913/2438417‑2438671
                                                                                                                                   |                                                                                                                           
gcgcGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGAt                                                                                                                      >  1:1029851/1‑139 (MQ=255)
  gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                    >  2:76899/1‑139 (MQ=255)
  gcgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                    >  1:278985/1‑139 (MQ=255)
          accgaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTaa                                                                                                            >  2:133753/1‑139 (MQ=255)
             gaccgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTc                                                                                                         >  1:488315/1‑139 (MQ=255)
                    cTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaa                                                                                                  <  2:488315/139‑1 (MQ=255)
                                    gCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTg                                                                                  <  1:76899/139‑1 (MQ=255)
                                      acaATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGt                                                                                <  1:468000/139‑1 (MQ=255)
                                      acaATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGt                                                                                <  1:48865/139‑1 (MQ=255)
                                                aaaaCAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTa                                                                      <  1:133753/139‑1 (MQ=255)
                                                 aaaCAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                    <  2:701888/92‑1 (MQ=255)
                                                 aaaCAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                    >  1:701888/1‑92 (MQ=255)
                                                        tGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                   >  1:285595/1‑102 (MQ=255)
                                                        tGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                   <  2:285595/102‑1 (MQ=255)
                                                           tGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCa                                                           <  2:525925/139‑1 (MQ=255)
                                                               cAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACg                                                       >  2:582293/1‑139 (MQ=255)
                                                                        ggcggcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                   <  2:943984/70‑1 (MQ=255)
                                                                        ggcggcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                   >  1:119971/1‑70 (MQ=255)
                                                                        ggcggcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                   <  2:119971/70‑1 (MQ=255)
                                                                        ggcggcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                   >  1:943984/1‑70 (MQ=255)
                                                                            gcAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCaaa                                          >  1:197837/1‑139 (MQ=255)
                                                                                  ggggAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGc                                    <  1:582293/139‑1 (MQ=255)
                                                                                              aaaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                        <  2:763361/139‑1 (MQ=255)
                                                                                                        ggCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGAt              <  2:982737/139‑1 (MQ=255)
                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCa                                                                                      >  2:514969/1‑61 (MQ=255)
                                                                                                              cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCa                                                                                      <  1:514969/61‑1 (MQ=255)
                                                                                                                       aaCTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAaaca  <  1:72487/136‑1 (MQ=255)
                                                                                                                       aaCTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAaaca  >  2:72487/1‑136 (MQ=255)
                                                                                                                           ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTg                                                                                  >  1:842289/1‑52 (MQ=255)
                                                                                                                           ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTg                                                                                  <  2:842289/52‑1 (MQ=255)
                                                                                                                           ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                                >  1:868744/1‑70 (MQ=255)
                                                                                                                           ttCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGgcg                                                                <  2:868744/70‑1 (MQ=255)
                                                                                                                                   |                                                                                                                           
GCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAACA  >  NC_000913/2438417‑2438671

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