Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	826,574	C→T	18.1%	D168N (GAT→AAT)	yvyE ←	hypothetical membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	826,574	0	C	T	18.1%	51.3 / 3.4	22	D168N (GAT→AAT)	yvyE	hypothetical membrane protein
Reads supporting (aligned to +/- strand):  ref base C (12/6);  new base T (2/2);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGTATATCGAAAGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGTTCGACACGCCCGATATCATAAATGA  >  CP000730/826442‑826708                                                                                                                                     |                                                                                                                                       ttGGTATATCGAAAGGTAACTGCTTTGGAGTTTCTTCAGTCAGATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAg                                                                                                                                >  1:227454/1‑141 (MQ=255) ttGGTATATCGAAAGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGAGTAAg                                                                                                                                >  1:189451/1‑141 (MQ=255)  tGGTATATCGAAAGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGt                                                                                                                               <  2:121962/141‑1 (MQ=255)            aaaGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCtt                                                                                                                     <  2:150999/141‑1 (MQ=255)             aaGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTa                                                                                                                    >  1:207818/1‑141 (MQ=255)              aggcaaCTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCAGCAGTCATTTAATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATTAGTGTAAGTTTGGTTTCTTaa                                                                                                                   >  1:200939/5‑141 (MQ=255)              aGGCAACTGCTTTGGAGTATCTTCAGTCAAATCGAAATTGTTAGCAGGAATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATTAGTGTAAGTTTGGTTTCTTaa                                                                                                                   >  2:92287/1‑141 (MQ=255)               ggTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAAc                                                                                                                  <  1:243126/141‑1 (MQ=255)               ggTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAAc                                                                                                                  <  1:243151/141‑1 (MQ=255)                    cTGCTTTGGCGTTTCTTCAGTCAAATCGAAAGTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCCCTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTGTCTTAACATaaa                                                                                                             >  2:140515/1‑141 (MQ=255)                                                 aaTTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAAtt                                                                                >  2:110209/1‑141 (MQ=255)                                                 aaTTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAAtt                                                                                >  2:79554/1‑141 (MQ=255)                                                             tCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGAt                                                                    <  2:4342/141‑1 (MQ=255)                                                                               taataaACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACta                                                  >  1:54403/1‑141 (MQ=255)                                                                                        cTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATTAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTa                                         <  2:200939/141‑1 (MQ=255)                                                                                        cTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATTAGTGTAAGTTTGGTTTATTAACATCAAATGAGTTGAAGCTAATACATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTa                                         <  1:92287/141‑1 (MQ=255)                                                                                             tAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCt                                    >  1:168583/1‑141 (MQ=255)                                                                                                       tAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGtt                          >  2:215979/1‑141 (MQ=255)                                                                                                               tCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGTTCGACACGc                  <  1:249249/141‑1 (MQ=255)                                                                                                               tCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGTTCGACACGc                  >  1:187699/1‑141 (MQ=255)                                                                                                                            ttACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATACTTCAGTGTTACTATGCATGGTACTGCTTCTCGTAGTTCGCCACGCCCAATATCATaaa     >  1:208439/1‑141 (MQ=255)                                                                                                                               ccTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGTTCGACACGCCCGATATCATAAATGa  >  2:222056/1‑141 (MQ=255)                                                                                                                                     |                                                                                                                                       TTGGTATATCGAAAGGTAACTGCTTTGGAGTTTCTTCAGTCAAATCGAAATTTCCTGCAGTCATTTGATTTAAAAAGTTAATAAACGCTTCATAGTCACTTTTAACGACATCGATATAGTAGCTTACCTTATCAGTGTAAGTTTGGTTTCTTAACATAAAATGAGTTGAAGCTAATTCATATTCAAATTTACCAGTTTGATCATAATTCAGTGTTACTATACATGGTACTGCTTCTCGTAGTTCGACACGCCCGATATCATAAATGA  >  CP000730/826442‑826708

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