Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,964,697	G→A	L118L (CTC→CTT)	flhA ←	putative flagellar export pore protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,964,697	0	G	A	100.0%	101.1 / NA	32	L118L (CTC→CTT)	flhA	putative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/13);  total (19/13)

ATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACA  >  NC_000913/1964572‑1964815                                                                                                                              |                                                                                                                       aTACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGa                                                                                                           >  1:201285/1‑139 (MQ=255)                 aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGaa                                                                                                                   >  1:1072607/1‑115 (MQ=255)                 aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGaa                                                                                                                   <  2:1072607/115‑1 (MQ=255)                 aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCg                                                                                                            >  1:1353489/1‑122 (MQ=255)                 aGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCg                                                                                                            <  2:1353489/122‑1 (MQ=255)                                  gCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcgc                                                                          >  2:128606/1‑139 (MQ=255)                                      tACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGt                                                                      <  1:340115/139‑1 (MQ=255)                                        cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTAt                                                                    >  2:227534/1‑139 (MQ=255)                                           ccGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGcc                                                                 >  2:122318/1‑139 (MQ=255)                                             ggCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCt                                                               <  2:201285/139‑1 (MQ=255)                                                  cTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCAt                                                          <  1:128606/139‑1 (MQ=255)                                                   tttGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCAtt                                                         <  1:122318/139‑1 (MQ=255)                                                            gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                >  2:972649/1‑139 (MQ=255)                                                            gACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATa                                                >  2:1266862/1‑139 (MQ=255)                                                                             tCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGtt                               <  1:446255/139‑1 (MQ=255)                                                                                  agaATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTg                          <  1:972649/139‑1 (MQ=255)                                                                                      tGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAg                      >  1:386662/1‑139 (MQ=255)                                                                                       gACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCg                                                                                                            >  2:725659/1‑52 (MQ=255)                                                                                       gACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCg                                                                                                            <  1:725659/52‑1 (MQ=255)                                                                                              accacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACa              >  1:1344205/1‑139 (MQ=255)                                                                                              accacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACa              >  1:1022551/1‑139 (MQ=255)                                                                                                cacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                <  2:1314697/87‑1 (MQ=255)                                                                                                cacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                >  1:1314697/1‑87 (MQ=255)                                                                                                  cGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCTCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGt          >  2:838781/1‑139 (MQ=255)                                                                                                  cGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                        >  1:333464/1‑125 (MQ=255)                                                                                                  cGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGcc                        <  2:333464/125‑1 (MQ=255)                                                                                                     tACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACg                                                                                                        <  2:457888/42‑1 (MQ=255)                                                                                                     tACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACg                                                                                                        >  1:457888/1‑42 (MQ=255)                                                                                                       ccGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTaa     >  1:194331/1‑139 (MQ=255)                                                                                                          aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAaaca  >  2:899854/1‑139 (MQ=255)                                                                                                             gCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTa      <  2:91644/132‑1 (MQ=255)                                                                                                             gCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTa      >  1:91644/1‑132 (MQ=255)                                                                                                                              |                                                                                                                       ATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACA  >  NC_000913/1964572‑1964815

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