Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,005,319	G→A	G6D (GGC→GAC)	fliS →	flagellar biosynthesis protein FliS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,005,319	0	G	A	100.0%	170.1 / NA	50	G6D (GGC→GAC)	fliS	flagellar biosynthesis protein FliS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/21);  total (29/21)

AACAATTTACCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGGCACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGACAAC  >  NC_000913/2005184‑2005449                                                                                                                                        |                                                                                                                                   aaCAATTTACCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACAcc                                                                                                                                <  1:809154/140‑1 (MQ=255)    aaTTTACCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAg                                                                                                                             >  2:570852/1‑140 (MQ=255)      tttACCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGc                                                                                                                           >  1:499624/1‑140 (MQ=255)         aCCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTa                                                                                                                        >  1:905927/1‑140 (MQ=255)           ccAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATg                                                                                                                      >  2:450235/1‑140 (MQ=255)              aCTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGcac                                                                                                                   >  2:157547/1‑140 (MQ=255)                 ggACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACaaa                                                                                                                <  2:499624/140‑1 (MQ=255)                       tttAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCg                                                                                                          >  2:893215/1‑140 (MQ=255)                       tttAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCg                                                                                                          >  1:153023/1‑140 (MQ=255)                        ttAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGt                                                                                                         >  2:661858/1‑140 (MQ=255)                           aTGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAAGATTCATCGGGAGACAGGTCATGTGCGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGa                                                                                                      <  2:508335/140‑1 (MQ=255)                                 tCGTTAAACAGCCCCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAgcgc                                                                                                >  1:653989/1‑140 (MQ=255)                                  cGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGcc                                                                                               >  1:607928/1‑140 (MQ=255)                                   gTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCg                                                                                              >  2:424721/1‑140 (MQ=255)                                             aCCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAgcgc                                                                                    >  1:492157/1‑140 (MQ=255)                                             aCCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAgcgc                                                                                    >  2:792641/1‑140 (MQ=255)                                                      tACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCcagca                                                                           >  1:606636/1‑140 (MQ=255)                                                       aCTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCcagcag                                                                          >  1:45012/1‑140 (MQ=255)                                                         ttAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCcagcagca                                                                        <  2:905927/140‑1 (MQ=255)                                                         ttAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCcagcagca                                                                        <  1:157547/140‑1 (MQ=255)                                                               cagcagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGt                                                                  <  1:450235/140‑1 (MQ=255)                                                                 gcagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCa                                                                <  1:836417/140‑1 (MQ=255)                                                                 gcagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCa                                                                <  1:35397/140‑1 (MQ=255)                                                                 gcagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCa                                                                <  2:45012/140‑1 (MQ=255)                                                                  cagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCcagccagc                                                                            >  1:574340/1‑127 (MQ=255)                                                                  cagTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCcagccagc                                                                            <  2:574340/127‑1 (MQ=255)                                                                    gTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCa                                                             <  1:424721/140‑1 (MQ=255)                                                                     ttCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCAt                                                            <  2:606636/140‑1 (MQ=255)                                                                     ttCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCAt                                                            >  2:471994/1‑140 (MQ=255)                                                                      tCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGATCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATg                                                           <  1:231136/140‑1 (MQ=255)                                                                           aacaacAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTAtt                                                      <  1:661858/140‑1 (MQ=255)                                                                             caacaGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTg                                                    <  2:753190/140‑1 (MQ=255)                                                                                        ctaaGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTg                                         >  1:78692/3‑140 (MQ=255)                                                                                        ccAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTg                                         >  2:271789/1‑140 (MQ=255)                                                                                              aaGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCa                                   <  2:153023/140‑1 (MQ=255)                                                                                              aaGAAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTCTTTGATGGAGTGCTGAGCGCa                                   <  2:653989/140‑1 (MQ=255)                                                                                               aGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCAc                                  <  2:78692/140‑1 (MQ=255)                                                                                               aGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCAACATGCTATTTGATGGAGTGCTGAGCGCAc                                  >  1:876882/1‑140 (MQ=255)                                                                                                           cGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAgagcga                      >  1:511225/1‑140 (MQ=255)                                                                                                           cGGGAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAgagcga                      >  1:176079/1‑140 (MQ=255)                                                                                                            gggAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAgagcgag                     <  2:511225/140‑1 (MQ=255)                                                                                                            gggAGACAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAgagcgag                     <  1:893215/140‑1 (MQ=255)                                                                                                              gagaCAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGcc                   >  1:372639/1‑140 (MQ=255)                                                                                                              gagaCAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGcc                   >  2:526420/1‑140 (MQ=255)                                                                                                              gagaCAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGcc                   >  2:421652/1‑140 (MQ=255)                                                                                                                 aCAGGTCATGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGt                >  2:779709/1‑140 (MQ=255)                                                                                                                        aTGTACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCa         <  2:607928/140‑1 (MQ=255)                                                                                                                           tACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGa      <  2:855821/140‑1 (MQ=255)                                                                                                                           tACGCGGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGa      <  1:779709/140‑1 (MQ=255)                                                                                                                               cgGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGacaac  >  1:293846/1‑140 (MQ=255)                                                                                                                               cgGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGacaac  >  2:547163/1‑140 (MQ=255)                                                                                                                               cgGCAAAAGACACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGacaac  >  2:49560/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                   AACAATTTACCCAACTGGACGTTTTAATGACCTCGTTAAACAGCACCAGCAGCTACTTAACGCAGCAGTTCGAAAACAACAGTAATTCCAAGTAAGCAATATTCATCGGGAGACAGGTCATGTACGCGGCAAAAGGCACCCAGGCCTATGCACAAATTGGCGTCGAAAGCGCCGTAATGAGCGCCAGCCAGCAGCAGCTGGTCACCATGCTATTTGATGGAGTGCTGAGCGCACTGGTTAGAGCGAGCCTGTTTATGCAGGACAAC  >  NC_000913/2005184‑2005449

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