Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,101	C→T	*339* (TAG→TAA)	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,101	0	C	T	93.8%	41.5 / ‑5.7	16	*339* (TAG→TAA)	lplA	lipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (8/7);  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

ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTG  >  NC_000913/4622966‑4623223                                                                                                                                        |                                                                                                                           aTCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAc                                                                                                                         >  1:755155/1‑139 (MQ=255) aTCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAc                                                                                                                         >  2:447035/1‑139 (MQ=255)  tCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAcc                                                                                                                        >  1:616301/1‑139 (MQ=255)                          aaaCGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCg                                                                                                <  1:780575/139‑1 (MQ=255)                                                tATCAGCGAGGCGTTTCGCCAAATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          <  2:616301/139‑1 (MQ=255)                                                         ggCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCgg                                                                 >  2:12273/1‑139 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                               <  1:447035/139‑1 (MQ=255)                                                                        tGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCCGTTCCGGGAAGTCAACCAACAg                                                  >  2:85734/1‑139 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAgcg                                                <  2:755155/139‑1 (MQ=255)                                                                                   tGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACt                                       <  1:403728/139‑1 (MQ=255)                                                                                    gCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTc                                      >  1:610125/1‑139 (MQ=255)                                                                                                       aTCGAAGAGAAAGTTGTCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTgcg                   <  2:1051655/139‑1 (MQ=255)                                                                                                                aaaGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACAAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAgg          <  2:610125/139‑1 (MQ=255)                                                                                                                        ccGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTg  >  1:35288/1‑139 (MQ=255)                                                                                                                                     tAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGt                                                                                     >  1:507367/1‑43 (MQ=255)                                                                                                                                     tAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGt                                                                                     <  2:507367/43‑1 (MQ=255)                                                                                                                                        |                                                                                                                           ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTG  >  NC_000913/4622966‑4623223

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