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	100.0%	56.7 / NA	19	*339* (TAG→TAA)	lplA	lipoate‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (11/8);  total (11/8)

ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCA  >  NC_000913/4622966‑4623225                                                                                                                                        |                                                                                                                             aTCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAc                                                                                                                           >  2:218353/1‑139 (MQ=255)     cGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTa                                                                                                                       >  1:58275/1‑139 (MQ=255)      gCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTAc                                                                                                                      <  2:10451/139‑1 (MQ=255)              agcagTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCg                                                                                                              >  2:510834/1‑139 (MQ=255)                           aaCGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGa                                                                                                 >  2:607081/1‑139 (MQ=255)                                                tATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCg                                                                                                              <  2:664523/105‑1 (MQ=255)                                                tATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCg                                                                                                              >  1:664523/1‑105 (MQ=255)                                                        aGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCg                                                                    >  2:643123/1‑139 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                            >  1:164566/1‑128 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                            <  2:164566/128‑1 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                                 <  2:785058/139‑1 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                                 >  1:335962/1‑139 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                            >  2:584475/1‑113 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                            <  1:584475/113‑1 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAgcg                                                  <  2:383728/139‑1 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAgcg                                                  >  2:399676/1‑139 (MQ=255)                                                                                                       aTCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTgcg                     <  2:146476/139‑1 (MQ=255)                                                                                                             gagaAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTAc               >  1:365000/1‑139 (MQ=255)                                                                                                                          gCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCa  <  2:365000/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                             ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCA  >  NC_000913/4622966‑4623225

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