Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	42,108	G→A	intergenic (‑177/‑295)	caiT ← / → fixA	L‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	42,108	0	G	A	93.3%	45.0 / ‑5.8	15	intergenic (‑177/‑295)	caiT/fixA	L‑carnitine:gamma‑butyrobetaine antiporter/putative electron transfer flavoprotein FixA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/6);  minor base C (0/1);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACAT  >  NC_000913/41863‑42296                                                                                                                                                                                                                                                      |                                                                                                                                                                                             gCCGACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTGCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCg                                                                                                                                                                                          <  1:48634/250‑1 (MQ=255)  ccGACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGtt                                                                                                                                                                                        >  2:169107/1‑251 (MQ=255)      cgaTTATTAACGGCGGAAAGAAAACCTGCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCGCTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTGATTTTTGTTAACATTTAATATAAGTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTGGTTTCGGGTGAATAGGCGGCGtttttt                                                                                                                                                                                    <  2:146813/251‑1 (MQ=255)                                   ggTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGtt                                                                                                                                                       >  2:50594/1‑251 (MQ=255)                                                                    aTGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATACGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCtgtgcttctccagccgtggtcgcgtct                                                                                                                      >  2:150679/1‑227 (MQ=255)                                                                                         aaCGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATa                                                                                                 <  1:50594/251‑1 (MQ=255)                                                                                                                  tttgtttTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGTTCCATAAAACAATATTGAAAATTTCTTTTTg                                                                        >  1:192336/1‑251 (MQ=255)                                                                                                                       tttATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGTACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTAc                                                                    <  2:243492/250‑1 (MQ=255)                                                                                                                        ttATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACg                                                                   <  1:298641/250‑1 (MQ=255)                                                                                                                                                  tattattaACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTGATGCTCTGTTGGGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTGGCTACGCCGTGTTTTCAATATTGGTGAGGatc                                         <  1:194991/250‑3 (MQ=255)                                                                                                                                                           ccTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAatat                               >  1:172998/1‑251 (MQ=255)                                                                                                                                                           ccTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAata                                >  1:80255/1‑250 (MQ=255)                                                                                                                                                           ccTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTGTGCAATAGTGGTGAGGAACTTAACAata                                <  2:80255/250‑1 (MQ=255)                                                                                                                                                                                        aaTGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGAGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGaaat  >  2:422893/1‑248 (MQ=255)                                                                                                                                                                                                   aaTAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAAGGCGTTTTTTCGTTAATTTTGATTAATAATCa                                                                                                                                                              >  1:426198/1‑84 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                             GCCGACGATTATTAACGGCGGAAAGAAAACCTTCGGTTCTATTCCCGTTTTTCTCTTTTCATTCTTCATGAGTTAATTCCACTGTGAAAACGAATATTTATTTTGCGTTCCCGTTTGTTTTATTTTTGTTAACATTTAATATAATTATTATTAACCTCGTGGACGCGTTAATGGCTAACTCATAATGGGTATTCAATAAGCTGTATTCTGTGATTGGTATCACATTTTTGTTTCGGGTGAATAGAGGGCGTTTTTTCGTTAATTTTGATTAATAATCAGTTTGTTATGCTCTGTTGTGAGTAAAAAATAACATCTGACTTTCAATATTGGTGATCCATAAAACAATATTGAAAATTTCTTTTTGCTACGCCGTGTTTTCAATATTGGTGAGGAACTTAACAATATTGAAAGTTGGATTTATCTGCGTGTGACAT  >  NC_000913/41863‑42296

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