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_000913	1,139,406	0	T	.	100.0%	80.9 / NA	20	coding (1029/1644 nt)	flgK	flagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (6/14);  total (6/14)

CTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTT  >  NC_000913/1139278‑1139522                                                                                                                                 |                                                                                                                     cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                           >  2:569750/1‑139 (MQ=255)     cgcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                       <  1:569750/139‑1 (MQ=255)       cTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGc                                                                                                     <  2:569272/139‑1 (MQ=255)         ggCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCa                                                                                                   <  1:409797/139‑1 (MQ=255)           cATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACg                                                                                                 <  2:775777/139‑1 (MQ=255)                                       caAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                     <  1:734525/139‑1 (MQ=255)                                          aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAg                                                                  >  1:498202/1‑139 (MQ=255)                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAGGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcgcg                                                                               <  2:580199/110‑1 (MQ=255)                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcgcg                                                                               >  1:580199/1‑110 (MQ=255)                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                  <  2:464274/139‑1 (MQ=255)                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                  <  2:498202/139‑1 (MQ=255)                                                               cGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGata                                             <  1:696890/139‑1 (MQ=255)                                                               cGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGata                                             <  2:247140/139‑1 (MQ=255)                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAAt                                                        >  1:735037/1‑126 (MQ=255)                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAAt                                                        <  2:735037/126‑1 (MQ=255)                                                                            gATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCa                                <  1:643124/139‑1 (MQ=255)                                                                                    tGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCAccc                        >  1:625420/1‑139 (MQ=255)                                                                                                gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                  <  2:106328/101‑1 (MQ=255)                                                                                                gCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                  >  1:106328/1‑101 (MQ=255)                                                                                                          cTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACttt  <  1:663826/139‑1 (MQ=255)                                                                                                                                 |                                                                                                                     CTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTT  >  NC_000913/1139278‑1139522

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