Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,139,597	T→C	V407A (GTT→GCT)	flgK →	flagellar hook‑filament junction protein 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,139,597	0	T	C	100.0%	134.6 / NA	39	V407A (GTT→GCT)	flgK	flagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/21);  total (18/21)

AAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGTTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCT  >  NC_000913/1139464‑1139732                                                                                                                                      |                                                                                                                                        aaaaTCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGa                                                                                                                                   <  1:754899/140‑1 (MQ=255)        cGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTc                                                                                                                            >  2:690162/1‑140 (MQ=255)        cGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTc                                                                                                                            <  1:767247/140‑1 (MQ=255)              ataataATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGTATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTg                                                                                                                      <  1:92492/140‑1 (MQ=255)              ataataATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTg                                                                                                                      <  2:315804/140‑1 (MQ=255)               taataaTCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGa                                                                                                                     >  1:821317/1‑140 (MQ=255)                     tCAGTGGCAGGTCGCCCGCCTTGCCAGCAATCCCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAg                                                                                                               <  2:752688/140‑1 (MQ=255)                         tGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAg                                                                                                           >  2:810917/1‑140 (MQ=255)                         tGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAg                                                                                                           >  1:584747/1‑140 (MQ=255)                         tGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAg                                                                                                           >  1:732444/1‑140 (MQ=255)                            cAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGa                                                                                                                     <  1:250792/127‑1 (MQ=255)                            cAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGa                                                                                                                     >  2:250792/1‑127 (MQ=255)                              ggTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGACATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACg                                                                                                      <  2:313675/140‑1 (MQ=255)                                tCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGcc                                                                                                    <  2:584747/140‑1 (MQ=255)                                   cccGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCa                                                                                                   <  2:501403/138‑1 (MQ=255)                                   cccGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCa                                                                                                   >  1:501403/1‑138 (MQ=255)                                   cccGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATc                                                                                                 >  2:492916/1‑140 (MQ=255)                                        cTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCaa                                                                                            >  2:386968/1‑140 (MQ=255)                                           gCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACAt                                                                                         <  2:821317/140‑1 (MQ=255)                                           gCCAGCAATACCACTTTTACGGTGACGCCGGACGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACAt                                                                                         <  2:732444/140‑1 (MQ=255)                                             cAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATgg                                                                                       <  1:127955/140‑1 (MQ=255)                                             cAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATgg                                                                                       <  1:690162/140‑1 (MQ=255)                                             cAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATgg                                                                                       <  1:27102/140‑1 (MQ=255)                                              aGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATgg                                                                                       <  1:595347/139‑1 (MQ=255)                                              aGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATgg                                                                                       >  2:595347/1‑139 (MQ=255)                                                                  gACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAg                                                                  <  1:810917/140‑1 (MQ=255)                                                                                  ggTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGgc                                                  <  1:386968/140‑1 (MQ=255)                                                                                            cATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCgaagaa                                        >  2:578915/1‑140 (MQ=255)                                                                                                      tCTGGAGTTGACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTg                              >  1:637932/1‑140 (MQ=255)                                                                                                               gACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATa                     >  2:206352/1‑140 (MQ=255)                                                                                                               gACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATa                     >  2:526764/1‑140 (MQ=255)                                                                                                               gACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATa                     >  1:462271/1‑140 (MQ=255)                                                                                                               gACGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGAGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATa                     >  2:800782/1‑140 (MQ=255)                                                                                                                aCGTTTACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATaa                    <  2:462271/140‑1 (MQ=255)                                                                                                                    ttACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGc                >  2:139448/1‑140 (MQ=255)                                                                                                                     tACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCa               <  2:646370/140‑1 (MQ=255)                                                                                                                     tACAGGAACGCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCa               <  1:492916/140‑1 (MQ=255)                                                                                                                              gCCTGCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAgg      <  2:637932/140‑1 (MQ=255)                                                                                                                                  gCCGCTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCct  >  1:721967/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                        AAAATCTCGTTCGATAATAATCAGTGGCAGGTCACCCGCCTTGCCAGCAATACCACTTTTACGGTGACGCCGGATGCCAACGGTAAAGTGGCATTTGATGGTCTGGAGTTGACGTTTACAGGAACGCCTGCCGTTAACGACAGCTTCACGCTGAAACCAGTAAGTGACGCCATCGTCAACATGGATGTATTAATCACCGACGAAGCGAAAATAGCGATGGCGAGCGAAGAAGATGCGGGTGATAGCGATAACCGCAACGGTCAGGCCCT  >  NC_000913/1139464‑1139732

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