Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,965,120	G→A	Q358* (CAG→TAG)	flhB ←	flagellar biosynthesis protein FlhB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,965,120	0	G	A	100.0%	109.7 / NA	33	Q358* (CAG→TAG)	flhB	flagellar biosynthesis protein FlhB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/17);  total (16/17)

TTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTC  >  NC_000913/1964997‑1965255                                                                                                                            |                                                                                                                                        ttGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGc                                                                                                                         <  2:115768/140‑1 (MQ=255)  tgtgTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGcc                                                                                                                        >  2:737535/1‑140 (MQ=255)      tCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCg                                                                                                                    >  1:18952/1‑140 (MQ=255)           ttCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCa                                                                                                               <  1:1000863/140‑1 (MQ=255)           ttCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCa                                                                                                               <  1:696765/140‑1 (MQ=255)                      gggCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGAcc                                                                                                    <  2:620470/140‑1 (MQ=255)                      gggCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGAcc                                                                                                    <  2:575699/140‑1 (MQ=255)                          aGGCGCAGCATCGCGGCTAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAgg                                                                                                >  2:28944/1‑140 (MQ=255)                            gcgcAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGcc                                                                                              <  2:39132/140‑1 (MQ=255)                             cgcAGCATCGCGGCCAGGTTACTCATGGGTCTGTTGCTCCTTAATAAACTCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCa                                                                                             <  2:369550/140‑1 (MQ=255)                                            aGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACcg                                                                              <  2:617465/140‑1 (MQ=255)                                                 aCTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCgcggcg                                                                         <  1:737535/140‑1 (MQ=255)                                                       gggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGt                                                                   >  1:32785/1‑140 (MQ=255)                                                        ggTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGtt                                                                  <  2:819511/140‑1 (MQ=255)                                                           tggTTTCTCGTTAATAAGATCCGGGGCTTCCTGCCGCGGAAGATGAGTCGGGTGTACAGGCCGCTATCCACCAGCCAGACGCCAGCGTTTGAGTTGCCAGACCCAGGCTAGCGCTCCCGCCCGCGCGGCGTACAGTTGAc                                                               <  2:395056/139‑1 (MQ=255)                                                            ggTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGAcc                                                              >  1:704838/1‑140 (MQ=255)                                                                    gTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGAttt                                                      <  2:32785/140‑1 (MQ=255)                                                                              tCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACGGTTGACCCGGGATTTGTTGACCAAt                                            <  2:859799/140‑1 (MQ=255)                                                                              tCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAAt                                            >  1:682546/1‑140 (MQ=255)                                                                                  gggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTcc                                        >  2:519891/1‑140 (MQ=255)                                                                                   ggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCg                                       <  1:102026/140‑1 (MQ=255)                                                                                   ggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCg                                       >  1:829198/1‑140 (MQ=255)                                                                                   ggCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCg                                       <  1:212079/140‑1 (MQ=255)                                                                                        ccGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGt                                  >  1:193225/1‑140 (MQ=255)                                                                                                          tAGGTTGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGcc                <  1:28944/140‑1 (MQ=255)                                                                                                              ttGTACAGGGCGCTACCCACCAACCAGACCCCAGCGTCTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcg            <  1:384731/140‑1 (MQ=255)                                                                                                               tGTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcgg           >  1:130176/1‑140 (MQ=255)                                                                                                                gTACAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggc          >  1:719613/1‑140 (MQ=255)                                                                                                                  aCAGGGCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAgcggcgg        <  2:682546/140‑1 (MQ=255)                                                                                                                      ggCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATTTCCGCATTTCGATACATCGCTGGCGCCAGCGGCGGCGCt    >  2:384731/1‑140 (MQ=255)                                                                                                                      ggCGCTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCt    >  1:381657/1‑140 (MQ=255)                                                                                                                       gcgcTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCtt   >  2:911322/1‑140 (MQ=255)                                                                                                                        cgcTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCTCCTGTCGCTACGGCACTCGCGCCAGCGGCGGCGCTTc  >  1:491715/1‑140 (MQ=255)                                                                                                                        cgcTATCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTACGCCTGTTGATACAGCGCTCGCGCCAGCGGCGGCGCTTc  >  2:367645/1‑140 (MQ=255)                                                                                                                            |                                                                                                                                        TTGTGTCGATTTCAGGTTTGCGGGCAGGCGCAGCATCGCGGCCAGATTACTCATGGGTCGGTTTCTCGTTAATAAAATCCAGGGCTTCCGGCACCGGAAGATGAGTAGGTTGTACAGGGCGCTGTCCACCAGCCAGACGCCAGCGTTTCAGTTGCCAGACCCAGGCCAGCACTTCCGCCACCGCGGCGTACAGTTGACCCGGGATTTGTTGACCAATCTCCGCATGTCGATACAGCGCTCGCGCCAGCGGCGGCGCTTC  >  NC_000913/1964997‑1965255

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑