Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,139,406 Δ1 bp coding (1029/1644 nt) flgK → flagellar hook‑filament junction protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,139,4060T.100.0% 110.2 / NA 27coding (1029/1644 nt)flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (11/16);  total (11/16)

ACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGG  >  NC_000913/1139277‑1139536
                                                                                                                                 |                                                                                                                                  
aCTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGa                                                                                                                          >  1:545526/1‑139 (MQ=255)
aCTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGA‑TTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                         <  2:705596/139‑1 (MQ=255)
 cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                         <  2:487756/139‑1 (MQ=255)
     cgcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCTGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                                     <  2:545526/139‑1 (MQ=255)
     cgcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                                     <  1:405103/139‑1 (MQ=255)
                                cccAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTa                                                                                          <  2:488181/139‑1 (MQ=255)
                                 ccAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTAc                                                                                         >  1:466754/1‑139 (MQ=255)
                                     cacaAAGCCGGATTTGATGCTAACTGCGATGCCGGTGCAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCa                                                                                                                 <  2:506791/111‑1 (MQ=255)
                                     cacaAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCa                                                                                                                 >  1:506791/1‑111 (MQ=255)
                                         aaGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACa                                                                                 <  1:628591/139‑1 (MQ=255)
                                           gCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                               >  1:291957/1‑139 (MQ=255)
                                                        cTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCg                                                                  <  2:291957/139‑1 (MQ=255)
                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                                <  1:241627/139‑1 (MQ=255)
                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                                >  1:432371/1‑139 (MQ=255)
                                                             ggCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGa                                                             <  2:466754/139‑1 (MQ=255)
                                                                gATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGataa                                                          <  2:195672/139‑1 (MQ=255)
                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                         >  1:57345/1‑75 (MQ=255)
                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                         <  2:57345/75‑1 (MQ=255)
                                                                  tGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTg                                                                                                        <  2:471489/91‑1 (MQ=255)
                                                                  tGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTg                                                                                                        >  1:471489/1‑91 (MQ=255)
                                                                  tGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGataata                                                        >  2:51193/1‑139 (MQ=255)
                                                                    acGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATaataat                                                      <  2:205485/138‑1 (MQ=255)
                                                                                 cTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCa                                         >  2:273915/1‑139 (MQ=255)
                                                                                    tGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCAccc                                      <  1:51193/139‑1 (MQ=255)
                                                                                                   cgcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATa                       >  1:207722/1‑139 (MQ=255)
                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGa        >  2:748446/1‑139 (MQ=255)
                                                                                                                        aaaCAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCgg  <  1:631093/139‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                  
ACTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGG  >  NC_000913/1139277‑1139536

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