Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,139,092	A→G	S239G (AGT→GGT)	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,092	0	A	G	96.4%	95.3 / ‑6.1	28	S239G (AGT→GGT)	flgK	flagellar hook‑filament junction protein 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (13/14);  minor base C (1/0);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACC  >  NC_000913/1138968‑1139208                                                                                                                             |                                                                                                                     tCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAACGGTTACTTACTCGTACAGGGACGTACGgcg                                                                                                              >  1:233295/1‑133 (MQ=255)    gctgGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCgg                                                                                                           >  2:66042/1‑133 (MQ=255)             aCGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTggcg                                                                                                  <  2:475853/133‑1 (MQ=255)              cgcgATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTggcgg                                                                                                 <  1:454343/133‑1 (MQ=255)               gcgATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTggcggc                                                                                                >  2:555028/1‑133 (MQ=255)                   tCAACTGGTGAGCGATTTAAACCACATTCTTGTTGTAGAAGCCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGtt                                                                                            <  2:233295/133‑1 (MQ=255)                    cAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGttc                                                                                           >  1:332643/1‑133 (MQ=255)                                     aaaCCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCt                                                                          <  1:555028/133‑1 (MQ=255)                                     aaaCCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCt                                                                          <  2:425316/133‑1 (MQ=255)                                     aaaCCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCt                                                                          <  2:354159/133‑1 (MQ=255)                                                ttgGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGAACCTTCTCGTACGAc                                                               >  1:603790/1‑133 (MQ=255)                                                   gtgtAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGt                                                            <  1:565412/133‑1 (MQ=255)                                                              aGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTg                                                 >  1:535197/1‑133 (MQ=255)                                                              aGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTg                                                 >  1:547390/1‑133 (MQ=255)                                                              aGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTg                                                 >  2:607119/1‑133 (MQ=255)                                                                        aTGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCATCTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggc                                       <  2:627609/133‑1 (MQ=255)                                                                        aTGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggc                                       <  2:332643/133‑1 (MQ=255)                                                                             ggcACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggcaggca                                  <  1:66042/133‑1 (MQ=255)                                                                             ggcACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACggcaggca                                  <  1:607119/133‑1 (MQ=255)                                                                                         aTCACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATcc                      >  1:565750/1‑133 (MQ=255)                                                                                           cACGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCg                    <  2:535197/133‑1 (MQ=255)                                                                                            acgttGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCgg                   >  2:595225/5‑133 (MQ=255)                                                                                            aCGATGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTAATGGGACGGCAAGCGATATTGAGATCCCgg                   >  1:180490/1‑133 (MQ=255)                                                                                               aTGGCCAATGGTTACTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGgaga                <  2:565750/133‑1 (MQ=255)                                                                                                    cAATGGTTACTCACTCGTTCAGGGTGGTCCGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTa           <  2:180490/133‑1 (MQ=255)                                                                                                            aCTCACTGGTTCATGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACTGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATAc   >  2:396795/1‑133 (MQ=255)                                                                                                            aCTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGGAATATTGAGATCCCGGGGAAATTACTGAATAc   >  2:168935/1‑133 (MQ=255)                                                                                                             cTCACTGGTTCAGGGAGGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAAGATTGAGAGCCCAGAGAAATTAGTGAATAcc  >  1:181339/1‑133 (MQ=255)                                                                                                                             |                                                                                                                     TCTGCTGGATCAACGCGATCAACTGGTGAGCGAATTAAACCAGATTGTTGGTGTAGAAGTCAGCGTTCAGGATGGCGGCACTTATAACATCACGATGGCCAATGGTTACTCACTGGTTCAGGGAAGTACGGCGCGGCAACTGGCGGCAGTTCCTTCCAGCGCTGACCCTTCTCGTACGACTGTCGCTTATGTTGATGGGACGGCAGGCAATATTGAGATCCCGGAGAAATTACTGAATACC  >  NC_000913/1138968‑1139208

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