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% 132.8 / NA 32coding (1029/1644 nt)flgKflagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (15/17);  total (15/17)

GGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGT  >  NC_000913/1139280‑1139528
                                                                                                                              |                                                                                                                          
ggCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCg                                                                                                               >  2:982385/1‑139 (MQ=255)
      ggCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCa                                                                                                         >  1:124919/1‑139 (MQ=255)
                 aGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGAt                                                                                              <  1:187348/139‑1 (MQ=255)
                             cccAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGt                                                                                                             >  1:282679/1‑112 (MQ=255)
                             cccAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGt                                                                                                             <  2:282679/112‑1 (MQ=255)
                                aaCACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTg                                                                               >  1:585175/1‑139 (MQ=255)
                                 acacAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTgg                                                                              >  1:1120781/1‑139 (MQ=255)
                                       aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGAt                                                                                                                 <  1:512272/98‑1 (MQ=255)
                                       aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGAt                                                                                                                 >  2:512272/1‑98 (MQ=255)
                                       aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGt                                                                                                             >  1:730712/1‑102 (MQ=255)
                                       aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAg                                                                        >  2:444339/1‑139 (MQ=255)
                                       aGCCGGATATGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATTGGt                                                                                                             <  2:730712/102‑1 (MQ=255)
                                        gCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       <  1:1028552/139‑1 (MQ=255)
                                        gCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       >  2:416500/1‑139 (MQ=255)
                                        gCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                       <  2:271233/139‑1 (MQ=255)
                                                       aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGc                                                                                                           >  1:488108/1‑88 (MQ=255)
                                                       aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGc                                                                                                           <  2:488108/88‑1 (MQ=255)
                                                       aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                        <  1:416500/139‑1 (MQ=255)
                                                          ggCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGa                                                     <  2:1120781/139‑1 (MQ=255)
                                                                     tGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGt                                          <  1:349963/139‑1 (MQ=255)
                                                                      gAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACaaaa                                                               >  1:845146/1‑117 (MQ=255)
                                                                      gAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACaaaa                                                               <  2:845146/117‑1 (MQ=255)
                                                                                              cccGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCaa                 <  1:958906/139‑1 (MQ=255)
                                                                                               ccGCGGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAAt                <  1:872518/139‑1 (MQ=255)
                                                                                                cgcgGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATa               >  1:795991/1‑139 (MQ=255)
                                                                                                 gcgGTTCTGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATAc              <  2:585175/139‑1 (MQ=255)
                                                                                                        tGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                           >  1:287585/1‑71 (MQ=255)
                                                                                                        tGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                           <  2:287585/71‑1 (MQ=255)
                                                                                                        tGCAAAACACGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACtttt       >  1:588037/1‑139 (MQ=255)
                                                                                                             aaCACGAAAAACAAATGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGt  <  1:759621/139‑1 (MQ=255)
                                                                                                               cacGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcg                                                                                       <  2:565920/52‑1 (MQ=255)
                                                                                                               cacGAAAAACAAAGGGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcg                                                                                       >  1:565920/1‑52 (MQ=255)
                                                                                                                              |                                                                                                                          
GGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGT  >  NC_000913/1139280‑1139528

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