Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,136	C→T	E328K (GAG→AAG)	lplA ←	lipoate‑‑protein ligase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,623,136	0	C	T	94.1%	42.6 / NA	17	E328K (GAG→AAG)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (10/6);  minor base G (0/1);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAA  >  NC_000913/4623079‑4623382                                                          |                                                                                                                                                                                                                                                       cctggatattttaatcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGAGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAGACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAATTGCCGTTCCAcgccc                                                       >  2:335490/24‑250 (MQ=255)    ggatattttaatcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGccc                                                    >  2:223649/21‑251 (MQ=255)                cagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCCGTCGCCCGGCGAGGGCTTCCCGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGCAATGCAGTTCCACGCCGCCCCAGCTAAAGCGt                                        >  1:424977/9‑251 (MQ=255)                cagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGt                                        >  1:19389/9‑251 (MQ=255)                   cgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTc                                      >  1:538451/6‑250 (MQ=255)                   cgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCa                                     <  2:191239/246‑1 (MQ=255)                   cgagttACCTTACCGCCCCCGCTATCCACGTCGAAAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGAGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGGGATATGGCCTTTTTCAACGTCGAAATACAGTTCCACGCCGCCCCAGCTAAAGCGTTc                                      >  1:573596/6‑250 (MQ=255)                    gagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCAt                                    >  1:582504/5‑251 (MQ=255)                     agttcccTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATc                                   <  2:30660/246‑1 (MQ=255)                        tACCTTCCCGCCCCCGCTATCCACGTCCACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCTAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCca                                 <  2:582504/250‑1 (MQ=255)                        tACCTTACCGCCCCCGCTATCCACTTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGCCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATcc                                  >  1:355084/1‑249 (MQ=255)                         cccTTACCGCCCCCGATATCCACGTCGACAACTGCCGTAGCGCGTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCGTCCAGCGGCGCTGGGTTGAGGCTGTCGGTAAACACCTGGTCGCGGGTGATATGGTCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCTTTTATCcagc                               <  2:84018/250‑1 (MQ=255)                                 gCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGCAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAACTGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCCGCAGCTGCgc                       >  1:27708/1‑250 (MQ=255)                                 gCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCgg                       >  1:136349/1‑250 (MQ=255)                                    cccGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGaa                    <  2:254906/251‑1 (MQ=255)                                            ccacgtcgacaaCTTCCGTAGCTCTTGTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCTCTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCTGTGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGAGATGTCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGa             <  1:556682/240‑1 (MQ=255)                                                      aaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTGCGCACTCGTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGTGTTGAGGCTGTCGGTAAACACCTTGGCGCGGGTGATATGGCCTGTGTCAACGTCTAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGaa  <  1:389350/251‑1 (MQ=255)                                                          |                                                                                                                                                                                                                                                       AGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAA  >  NC_000913/4623079‑4623382

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