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%	118.0 / NA	37	L118L (CTC→CTT)	flhA	putative flagellar export pore protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/15);  total (22/15)

GCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAAC  >  NC_000913/1964570‑1964817                                                                                                                                |                                                                                                                         gCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGAcc                                                                                                               >  1:126434/1‑139 (MQ=255)       cATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGc                                                                                                         >  1:342102/1‑139 (MQ=255)        aTCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCt                                                                                                        >  2:299125/1‑139 (MQ=255)                                          cGCCCGGCACCTTTGGTCATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGccgccgccg                                                                               >  2:441168/1‑130 (MQ=255)                                          cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGccgccgccg                                                                               <  1:441168/130‑1 (MQ=255)                                          cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTAt                                                                      >  1:382707/1‑139 (MQ=255)                                          cGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTAt                                                                      >  2:24503/1‑139 (MQ=255)                                              cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                  <  2:15569/139‑1 (MQ=255)                                              cGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGccc                                                                  <  1:98188/139‑1 (MQ=255)                                                gCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCtt                                                                <  1:299125/139‑1 (MQ=255)                                                       tGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTaa                                                         >  1:199490/1‑139 (MQ=255)                                                                 cATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                             <  1:75022/109‑1 (MQ=255)                                                                 cATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCcgccgcg                                                                             >  2:75022/1‑109 (MQ=255)                                                                       gTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGaa                                         <  2:693191/139‑1 (MQ=255)                                                                         tgatgaTCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGc                                       <  2:199490/139‑1 (MQ=255)                                                                          gatgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGcc                                      >  2:221372/1‑139 (MQ=255)                                                                           atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGaa                                         >  2:579557/1‑135 (MQ=255)                                                                           atgatCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGaa                                         <  1:579557/135‑1 (MQ=255)                                                                                    agaATGACAAACACCATGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTg                            <  2:77499/139‑1 (MQ=255)                                                                                        tGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGccgccgccg                                                                               >  1:556141/1‑84 (MQ=255)                                                                                        tGACAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGccgccgccg                                                                               <  2:556141/84‑1 (MQ=255)                                                                                          aCAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGAc                      >  1:51107/1‑139 (MQ=255)                                                                                           cAAACACCACGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACg                     <  1:138857/139‑1 (MQ=255)                                                                                               caccacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCcgccgccgc                                                                              <  2:193992/78‑1 (MQ=255)                                                                                               caccacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCcgccgccgc                                                                              >  1:193992/1‑78 (MQ=255)                                                                                                 ccacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAg               >  1:592262/1‑139 (MQ=255)                                                                                                  cacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGc              >  1:769120/1‑139 (MQ=255)                                                                                                  cacGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGc              >  1:250013/1‑139 (MQ=255)                                                                                                   acGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCg             <  1:221372/139‑1 (MQ=255)                                                                                                    cGATACCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAACGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGt            >  1:19312/1‑139 (MQ=255)                                                                                                        aCCGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTa        <  2:342102/139‑1 (MQ=255)                                                                                                          cGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTaaa      >  2:385299/1‑139 (MQ=255)                                                                                                          cGATAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGACCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTaaa      >  2:552143/1‑139 (MQ=255)                                                                                                            aTAGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAaaca    >  1:652912/1‑139 (MQ=255)                                                                                                              aGCGAAATTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAaacaac  >  2:265517/1‑139 (MQ=255)                                                                                                                   aaTTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCa                                     <  2:83463/99‑1 (MQ=255)                                                                                                                   aaTTGCCACCAACAAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCa                                     >  1:83463/1‑99 (MQ=255)                                                                                                                                |                                                                                                                         GCATACCATCGAGAACAAAGCGCGCACCCACTTCTGCGATACGCCCGGCACCTTTGGTAATGACCATAAAGTTGATGATCACGAGAATGACAAACACCACGATACCGATAGCGAAATTGCCACCAACGAGGAAGTGACCGAACGCTTCGACCACCTTCCCTGCCGCCGCCGCGCCGGTATGCCCTTCCATTAAAATGATACGGGTTGAAGCCACGTTAAGTGCCAGACGCAACAGCGTGGTAAACAAC  >  NC_000913/1964570‑1964817

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