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

CTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGA  >  NC_000913/1139278‑1139530                                                                                                                                 |                                                                                                                             cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   >  2:322183/1‑139 (MQ=255) cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   <  1:1268124/139‑1 (MQ=255) cTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   >  1:511970/1‑139 (MQ=255)      gcTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTg                                                                                                              <  1:1358581/139‑1 (MQ=255)           cATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACg                                                                                                         >  2:1020404/1‑139 (MQ=255)                                   aaCACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCGATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTACCACGGTAACTGATGCCTCCGCGGTACTg                                                                                 >  2:994669/1‑139 (MQ=255)                                    acacAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                               <  1:1235032/108‑1 (MQ=255)                                    acacAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGt                                                                                                               >  2:1235032/1‑108 (MQ=255)                                          aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAg                                                                          >  1:1095402/1‑139 (MQ=255)                                          aGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAg                                                                          >  1:1168761/1‑139 (MQ=255)                                           gCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCCGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGa                                                                         <  2:511970/139‑1 (MQ=255)                                               gATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTAc                                                                     >  2:257014/1‑139 (MQ=255)                                                aTTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACa                                                                    <  2:1095402/139‑1 (MQ=255)                                                        cTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCg                                                            <  2:1339165/139‑1 (MQ=255)                                                          aaCGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGtt                                                          <  1:977763/139‑1 (MQ=255)                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcgcg                                                                                       >  1:426328/1‑103 (MQ=255)                                                                 aTGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcgcg                                                                                       <  2:426328/103‑1 (MQ=255)                                                                     cGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                             <  1:812869/109‑1 (MQ=255)                                                                     cGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGa                                                                             >  2:812869/1‑109 (MQ=255)                                                                                 cTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCa                                   >  2:1023694/1‑139 (MQ=255)                                                                                            gTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGcc                        <  1:1020404/139‑1 (MQ=255)                                                                                                   cgcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATa                 >  2:1048407/1‑139 (MQ=255)                                                                                                   cgcgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATa                 >  1:1152372/1‑139 (MQ=255)                                                                                                     cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   <  1:525831/39‑1 (MQ=255)                                                                                                     cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   <  2:1112631/39‑1 (MQ=255)                                                                                                     cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   >  2:525831/1‑39 (MQ=255)                                                                                                     cgGTTCTGCAAAACACGAAAAACAAAGG‑GACGTTGCGAt                                                                                                                   >  1:1112631/1‑39 (MQ=255)                                                                                                            gCAAAACACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTa        <  1:257014/139‑1 (MQ=255)                                                                                                                aaCACGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGt    <  1:1048407/139‑1 (MQ=255)                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcg                                                                                         <  1:1369066/52‑1 (MQ=255)                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCcg                                                                                         >  2:1369066/1‑52 (MQ=255)                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATAc                >  2:267103/1‑125 (MQ=255)                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATAc                <  1:267103/125‑1 (MQ=255)                                                                                                                  cacGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGa  >  1:603091/1‑139 (MQ=255)                                                                                                                    cGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCaa                   >  2:480426/1‑120 (MQ=255)                                                                                                                    cGAAAAACAAAGG‑GACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCaa                   <  1:480426/120‑1 (MQ=255)                                                                                                                                 |                                                                                                                             CTGGCGCTGGCATTTGCCGAGGCTTTCAACACCCAACACAAAGCCGGATTTGATGCTAACGGCGATGCCGGTGAAGATTTCTTTGCTATCGGTAAGCCCGCGGTTCTGCAAAACACGAAAAACAAAGGTGACGTTGCGATCGGTGCCACGGTAACTGATGCCTCCGCGGTACTGGCGACAGATTACAAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGA  >  NC_000913/1139278‑1139530

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