Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,005,021	(A)7→6	coding (1150/1407 nt)	fliD →	flagellar filament capping protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,005,015	0	A	.	100.0%	120.6 / NA	29	coding (1144/1407 nt)	fliD	flagellar filament capping protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (12/17);  total (12/17)

AACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACC  >  NC_000913/2004938‑2005095                                                                              |                                                                                 aactggaTGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGTaaa                                                                                >  1:356738/1‑80 (MQ=255)        tGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGc                                                                       <  1:356696/81‑1 (MQ=255)             aCAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCAc                                                                  <  1:356697/81‑1 (MQ=255)               aaaCTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGa                                                                >  2:356731/1‑81 (MQ=255)                 aCTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGacc                                                              <  2:356698/81‑1 (MQ=255)                  cTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGacca                                                             >  2:356743/1‑81 (MQ=255)                   tCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGaccac                                                            >  2:356741/1‑81 (MQ=255)                       cGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATc                                                        >  1:356734/1‑81 (MQ=255)                          tGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCgg                                                      >  1:356748/1‑80 (MQ=255)                           gCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCgg                                                      >  2:356742/1‑79 (MQ=255)                           gCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCgg                                                      >  1:356737/1‑79 (MQ=255)                            cACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGgcag                                                   >  1:356736/1‑81 (MQ=255)                            cACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGgc                                                     <  2:356699/79‑1 (MQ=255)                            cACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGgc                                                     <  2:356700/79‑1 (MQ=255)                             aCTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGgcagc                                                  <  2:356701/81‑1 (MQ=255)                              cTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGgcagca                                                 <  2:356702/81‑1 (MQ=255)                                        aTGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAg                                       <  1:356703/81‑1 (MQ=255)                                         tGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGt                                      <  2:356704/81‑1 (MQ=255)                                            cAGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTgg                                   <  2:356705/81‑1 (MQ=255)                                             aGCGGCGTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGc                                  >  1:356740/1‑81 (MQ=255)                                                 gcgTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTc                              <  2:356706/81‑1 (MQ=255)                                                  cgTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCg                             <  1:356708/81‑1 (MQ=255)                                                  cgTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCg                             <  1:356707/81‑1 (MQ=255)                                                   gTAGGTGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGa                            <  1:356709/81‑1 (MQ=255)                                                        tGCATTGATTGTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACg                       <  1:356710/81‑1 (MQ=255)                                                              gattgTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCat                  >  1:356745/1‑80 (MQ=255)                                                               attgTTGGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCatt                 <  2:356711/80‑1 (MQ=255)                                                                  gttgGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCattatt              <  2:356712/80‑1 (MQ=255)                                                                    tgGCGATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCAttatta             >  2:356756/1‑79 (MQ=255)                                                                        gATGGT‑AAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCg       >  2:356750/1‑81 (MQ=255)                                                                              aaaaaaaCCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTAcc  >  2:356746/1‑81 (MQ=255)                                                                              |                                                                                 AACTGGATGCCGACAAACTCACCGCTGCACTGAAAAAAGATGCCAGCGGCGTAGGTGCATTGATTGTTGGCGATGGTAAAAAAACCGGCATCACGACCACCATCGGCAGCAACCTGACCAGTTGGCTTTCGACAACGGGCATTATTAAAGCCGCTACC  >  NC_000913/2004938‑2005095

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