Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	CP006881	2,185,791	C→T	M1M (GTG→ATG) †	U712_10950 ←	Hypothetical Protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	2,185,791	0	C	T	90.0%	46.1 / ‑1.7	20	M1M (GTG→ATG) †	U712_10950	Hypothetical Protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (8/10);  total (9/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TTGTATAAAGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGTTGGATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATAC  >  CP006881/2185692‑2185926                                                                                                    |                                                                                                                                        gtaagtaAAGCACATCCCACCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGc                                                                                    <  1:134920/133‑1 (MQ=255)         aGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGaagaag                                                                           >  1:135544/1‑139 (MQ=255)         aGCACATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGaagaag                                                                           >  4:207313/1‑139 (MQ=255)           cacATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGt                                                                         <  5:73831/137‑1 (MQ=255)           cacATCCCCCCTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGt                                                                         <  4:408029/137‑1 (MQ=255)                ccccccTTTTCG‑TTTTTCTATCCGTTCTATAAGACAAGCTATGTTTTTTCGCCTCGAATTATGA‑‑‑‑‑‑‑‑‑‑‑‑‑TCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGAt                                                                    >  8:80215/1‑139 (MQ=255)                                                      gCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGt                                                                         >  5:35321/1‑111 (MQ=255)                                                      gCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGt                                                                         <  6:35321/111‑1 (MQ=255)                                                                                          atgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        >  3:287422/5‑139 (MQ=255)                                                                                          atgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        >  3:281509/5‑139 (MQ=255)                                                                                           tgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        <  3:315626/135‑1 (MQ=255)                                                                                           tgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        <  4:157767/135‑1 (MQ=255)                                                                                           tgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        <  6:228438/135‑1 (MQ=255)                                                                                           tgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        <  7:210510/135‑1 (MQ=255)                                                                                           tgatCTGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAg        <  8:250184/135‑1 (MQ=255)                                                                                                tGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGctcct                       >  4:291376/1‑116 (MQ=255)                                                                                                tGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGctcct                       <  3:291376/118‑3 (MQ=255)                                                                                                tGCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTGGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGTTATACAGGAAGCAATAc  >  3:82397/1‑139 (MQ=255)                                                                                                 gCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGtt                  >  3:291729/1‑122 (MQ=255)                                                                                                 gCATAAAAAATCCTGAAGTAAAA‑CTTCAGGATTTTTTCGAAAGCTTATTTAAATGCTCGAAGAAGGTTGGATAAGAAACGTGAATGGCATCCGTGTGCTCGATTTCAATCGGCTCCTCCGtt                  <  4:291729/122‑1 (MQ=255)                                                                                                    |                                                                                                                                        TTGTATAAAGCGCATCCCCCCTTTTCGTTTTTTCTTTTCTCTCTATAAGACAAGCTATGTTTTTTTGCCTCGGATTATGATTCTCACAGGCATTCTGCACAAAAAATCCTGAAGGAAAATCTTCAGGATTTGTTTGAAAGCTTATTTAAATGCTCGAAAAAGGTTGGATAAGAAACGTGAATGGCATCTGTCTGCTCAATTTCAATTGGCTCTTCAGTTATACACGAAGCGATAC  >  CP006881/2185692‑2185926

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