Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	W3110S.gb	1,135,828	T→C	71.4%	D150D (GAT→GAC)	flgF →	flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	W3110S.gb	1,135,828	0	T	C	71.4%	17.3 / 6.0	14	D150D (GAT→GAC)	flgF	flagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base C (2/8);  total (4/10)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

GCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAG  >  W3110S.gb/1135760‑1135889                                                                     |                                                              gCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGAccc                                                             >  1:405169/1‑71 (MQ=255) gCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGAccc                                                             >  1:725202/1‑71 (MQ=255)       ccGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCaaa                                                       <  1:1731250/71‑1 (MQ=255)                  ggAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTgcgc                                            <  1:910616/71‑1 (MQ=255)                      aTCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGt                                        <  1:971759/71‑1 (MQ=255)                           tATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATAc                                                    <  1:2078562/54‑1 (MQ=255)                             ttGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGt                                                 >  1:2354893/1‑55 (MQ=255)                             ttGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGt                                 >  1:1209827/1‑71 (MQ=255)                              tgctgtCGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTGCGCCAGt                                        <  1:2590864/63‑1 (MQ=255)                                   ccGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTa                             <  1:1265357/69‑1 (MQ=255)                                   ccGATGGCACAATCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTa                             <  1:1643676/69‑1 (MQ=255)                                      aTGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAg                                         <  1:307152/54‑1 (MQ=255)                                              aTCTCGGCGCTCAATCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGc                <  1:1678740/71‑1 (MQ=255)                                                 tcGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCAc             >  1:2578704/1‑71 (MQ=255)                                                          aaTCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCa              >  1:376107/1‑61 (MQ=255)                                                            tCCGGGCGACCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAg  <  1:86441/71‑1 (MQ=255)                                                                     |                                                              GCTGTGCCGGAAGGGGCGGAAATCACTATTGCTGCCGATGGCACAATCTCGGCGCTCAATCCGGGCGATCCGGCAAATACGGTTGCGCCAGTAGGGCGTCTTAAACTGGTGAAAGCCACGGGCAGCGAAG  >  W3110S.gb/1135760‑1135889

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