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% 147.6 / NA 41V200A (GTG→GCG) flkputative flagella assembly protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/23);  total (18/23)

CAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAA  >  NC_000913/2438415‑2438669
                                                                                                                                     |                                                                                                                         
cagcGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAg                                                                                                                      >  1:347579/1‑139 (MQ=255)
     cgAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGc                                                                                                                 <  2:154672/139‑1 (MQ=255)
        ggCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAg                                                                                                              >  1:927334/1‑139 (MQ=255)
           gaccgaccGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTa                                                                                                           >  2:489151/1‑139 (MQ=255)
                  cgTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAg                                                                                                    <  1:489151/139‑1 (MQ=255)
                     ccTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGa                                                                                                 <  1:835995/139‑1 (MQ=255)
                       tttATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                               <  2:347579/139‑1 (MQ=255)
                                  cggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCa                                                                                    >  2:311570/1‑139 (MQ=255)
                                   ggAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                   >  1:955061/1‑139 (MQ=255)
                                           aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCTGCAAACAGATCTGGCAATCGATGCTGGAACGTTCCGGGGCGTAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                           <  1:774798/139‑1 (MQ=255)
                                           aTGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                           <  1:311570/139‑1 (MQ=255)
                                                 tAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                  >  2:746239/1‑94 (MQ=255)
                                                 tAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATgg                                                                                                                  <  1:746239/94‑1 (MQ=255)
                                                                ccAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTa                                                                                        >  2:659223/1‑105 (MQ=255)
                                                                ccAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTa                                                                                        <  1:659223/105‑1 (MQ=255)
                                                                         cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                               <  1:954722/89‑1 (MQ=255)
                                                                         cggcggCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGaaa                                                                                               >  2:954722/1‑89 (MQ=255)
                                                                                           ccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                           <  2:7099/139‑1 (MQ=255)
                                                                                           ccAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                           <  2:927334/139‑1 (MQ=255)
                                                                                              gCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGa                        <  1:250116/139‑1 (MQ=255)
                                                                                                aaaCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                      <  2:955061/139‑1 (MQ=255)
                                                                                                  aCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTaa                                                                                       >  2:505775/1‑72 (MQ=255)
                                                                                                  aCAGATCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTaa                                                                                       <  1:505775/72‑1 (MQ=255)
                                                                                                      aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                   <  1:809915/72‑1 (MQ=255)
                                                                                                      aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCAt                                                                                   >  2:809915/1‑72 (MQ=255)
                                                                                                      aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACg                             <  1:739924/126‑1 (MQ=255)
                                                                                                      aTCTGGCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACg                             >  2:739924/1‑126 (MQ=255)
                                                                                                           gCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGtt                                                                                          >  1:225021/1‑60 (MQ=255)
                                                                                                           gCAATCGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGtt                                                                                          <  2:225021/60‑1 (MQ=255)
                                                                                                                cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                      >  1:776433/1‑123 (MQ=255)
                                                                                                                cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAAt                      <  2:776433/123‑1 (MQ=255)
                                                                                                                cGATGCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCt      >  1:882150/1‑139 (MQ=255)
                                                                                                                    gCTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTaaaa  <  2:882150/139‑1 (MQ=255)
                                                                                                                     cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACt                  >  1:95417/1‑122 (MQ=255)
                                                                                                                     cTGGAACTTTCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACt                  <  2:95417/122‑1 (MQ=255)
                                                                                                                              tCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                           >  1:472307/1‑104 (MQ=255)
                                                                                                                              tCCGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCt                           <  2:472307/104‑1 (MQ=255)
                                                                                                                                cGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                           <  1:320063/54‑1 (MQ=255)
                                                                                                                                cGGGGCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGAc                                                                           >  2:320063/1‑54 (MQ=255)
                                                                                                                                  gggCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCa                                                         <  2:128896/70‑1 (MQ=255)
                                                                                                                                  gggCGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCa                                                         >  1:128896/1‑70 (MQ=255)
                                                                                                                                     |                                                                                                                         
CAGCGCGAGGCGACCGACCGTCCTTTATTACCGGCGGAGCACAATGCGCTAAAACAGCTGGTGACCAAACTTGCGGCGGCAACGGGGGAACCCAGCAAACAGATCTGGCAATCGATGCTGGAACTTTCCGGGGTGAAAGATGGCGAGTTAATTCCAGCGAAACTGTTTAACCATCTGGTGACCTGGCTACAGGCGCGTCAGACGCTAAGCCAGCAAAATACGCCGACGCTGGAATCACTACAGATGACGCTAAAA  >  NC_000913/2438415‑2438669

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