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

ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTG  >  NC_000913/4622966‑4623223                                                                                                                                        |                                                                                                                           aTCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAc                                                                                                                         >  1:340451/1‑139 (MQ=255)  tCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATAGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTAcc                                                                                                                        >  2:974102/1‑139 (MQ=255)      gCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTAc                                                                                                                    >  2:8246/1‑139 (MQ=255)         ggCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGc                                                                                                                 <  2:46788/139‑1 (MQ=255)                     cTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGccatcca                                                                                                     <  1:1473347/139‑1 (MQ=255)                     cTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGccatcca                                                                                                     <  2:1607798/139‑1 (MQ=255)                           aaCGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGa                                                                                               >  2:1455273/1‑139 (MQ=255)                             cGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTGTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGACAGTTGCCCGCATGGGCGGGTAATTACCTTGCAGCCCCCGCCATCCATGCCGATa                                                                                             >  1:675626/1‑139 (MQ=255)                                                   cAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          <  2:287315/136‑1 (MQ=255)                                                   cAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTg                                                                       >  1:1175000/1‑139 (MQ=255)                                                   cAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGAGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          >  1:287315/1‑136 (MQ=255)                                                    aGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTa                                                                                                                     <  1:1505988/92‑1 (MQ=255)                                                    aGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTa                                                                                                                     >  2:1505988/1‑92 (MQ=255)                                                           cGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                               <  2:340451/139‑1 (MQ=255)                                                             tttCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAg                                                             <  2:1447778/139‑1 (MQ=255)                                                                        tGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAg                                                             >  1:15487/1‑128 (MQ=255)                                                                        tGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAg                                                             <  2:15487/128‑1 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGccatccat                                                                                                    <  1:308497/87‑1 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGccatccat                                                                                                    >  2:308497/1‑87 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCt                                                                        >  1:1189533/1‑115 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCt                                                                        <  2:1189533/115‑1 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAgcg                                                >  1:1198654/1‑139 (MQ=255)                                                                          cgcAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTC‑TTTTCCTGTTCCGGGAAGTCAACCAACAgcgc                                               >  1:654390/1‑139 (MQ=255)                                                                               gAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          <  2:211306/108‑1 (MQ=255)                                                                               gAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          >  1:211306/1‑108 (MQ=255)                                                                               gAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCg                                           <  2:90670/139‑1 (MQ=255)                                                                               gAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTC‑TTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGc                                          <  2:654390/139‑1 (MQ=255)                                                                                        gAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGt                                  <  1:1391249/139‑1 (MQ=255)                                                                                                  ggCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCatat                        <  2:1283473/139‑1 (MQ=255)                                                                                                    cAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCt                      <  1:1455273/139‑1 (MQ=255)                                                                                                      aaTCGAAGAGAAAGTTGCCCGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTgc                    <  2:1198654/139‑1 (MQ=255)                                                                                                                       cccGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          >  2:1462943/1‑68 (MQ=255)                                                                                                                       cccGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTc                                                                          <  1:1462943/68‑1 (MQ=255)                                                                                                                        ccGCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTg  >  2:533344/1‑139 (MQ=255)                                                                                                                          gCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                               >  1:473094/1‑76 (MQ=255)                                                                                                                          gCATGGGCGGGTAATTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGa                                                               <  2:473094/76‑1 (MQ=255)                                                                                                                                        |                                                                                                                           ATCCCGCTGGCGCAGCAGTTCTGGCAAACGGTCGCCAGAGATCCGCGTATCAGCGAGGCGTTTCGCCACATTGCGCAAGAAATGCCGGAAAAAATCCGGCAAATCGAAGAGAAAGTTGCCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTG  >  NC_000913/4622966‑4623223

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