Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,366,668	C→T	A36T (GCC→ACC)	dcuA ←	C4‑dicarboxylate transporter DcuA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,366,668	0	C	T	93.8%	47.9 / NA	16	A36T (GCC→ACC)	dcuA	C4‑dicarboxylate transporter DcuA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (7/8);  minor base G (1/0);  total (8/8)
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

ATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAT  >  NC_000913/4366430‑4366907                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                aTAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATCGCGGTAAGAACCAGCAc                                                                                                                                                                                                                                     >  1:206042/1‑251 (MQ=255)                                                    gggaaaagggTCCCGATCGGTGCGAGGAGTTTGATTTATTGGGGTTTAGGGCGCAGCAGTTTTTTTGTCGAATAACCCAAATAGCCCAGCCCGCCAGCAACCTGAATGGAAGAAATGGCGGCAAAAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaa                                                                                                                                                                                  <  1:17012/243‑1 (MQ=255)                                                       aaaaggctcccggtgggggcggggATCGTTATGTATTTCGGTTCCCGGCGCCGCATCTTTTCGGTCGGATAAACAAGATAGCCCAGCCCCCCAGCAACCTGAAGGGCAGAAATGGCGTCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaagaa                                                                                                                                                                               <  1:101058/229‑1 (MQ=255)                                                                 cGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTAACGGGTTTAACGCTAATAGCGGGAAGAACCAGCACCCACAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGccgg                                                                                                                                                                      >  2:166318/1‑248 (MQ=255)                                                                                                    cGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAACTAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGc                                                                                                                                  >  2:187928/1‑251 (MQ=255)                                                                                                                     tgtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAAc                                                                                                                 >  2:302870/1‑251 (MQ=255)                                                                                                                      gtctgataaacaagataccccaccccGCCAGCACCCTGCAGGGCAAAAATAGGGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTa                                                                                                                   <  1:9485/224‑1 (MQ=255)                                                                                                                               cccaaaTAGCCCAGACCGCCAGCACCCGGCAGGGCAAAAATGGCGGCGAAACCCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGatat                                                                                                       <  1:54597/246‑1 (MQ=255)                                                                                                                                  aGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAACCCAATACCTATTCCCCCCAATCTGGCGCCCAAGGAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATAtt                                                                                                      >  2:269248/1‑249 (MQ=255)                                                                                                                                        tCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGtt                                                                                              >  1:123282/1‑251 (MQ=255)                                                                                                                                            gACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAttt                                                                                          >  1:116740/1‑251 (MQ=255)                                                                                                                                                                                                        gatggagagggATGTGACAGGGTGTAACGCCAATAGCGGTAAGAACCAGCGTCCCCAATCCGCCTGGAAAACCAATACCTATTCGCCCCAAGCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATGAGCCTTCCTTGTTTTTGAACAAGTTGATATTAGATTGTTAGTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGtact                              <  2:62372/243‑1 (MQ=255)                                                                                                                                                                                                                     gTTACCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGAtt                 <  1:43993/251‑1 (MQ=255)                                                                                                                                                                                                                        aCCGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCAGTACGATTACt              <  1:228165/251‑1 (MQ=255)                                                                                                                                                                                                                          cGGGTTTAACGCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGt            <  1:92267/251‑1 (MQ=255)                                                                                                                                                                                                                                    gCCAATAGCGGTAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTtcat  >  1:7408/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                ATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAGAAAAAAGGCACGTCATCTGACGTGCCTTTTTTATTTGTACTACCCTGTACGATTACTGTTCGCTTTCAT  >  NC_000913/4366430‑4366907

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