Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,623,188	C→T	Q310Q (CAG→CAA)	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,188	0	C	T	100.0%	60.6 / NA	22	Q310Q (CAG→CAA)	lplA	lipoate‑‑protein ligase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (10/12);  total (10/12)

CCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTT  >  NC_000913/4623084‑4623427                                                                                                         |                                                                                                                                                                                                                                                atattttaatcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGCAAACACCTGGGCGCGGGGGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGccccc                                                                                               >  1:302782/19‑250 (MQ=255)        aatcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAg                                                                                        <  1:42890/240‑1 (MQ=255)         atcagcgagttACCTTCCCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGc                                                                                       <  1:486875/241‑1 (MQ=255)         atcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAg                                                                                              >  1:94597/11‑244 (MQ=255)         atcagcgagttACCTTACCGCCCCCGCTATCCACGTCGACAACTGCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAg                                                                                              <  2:94597/234‑1 (MQ=255)                   tACCTTACCGCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCca                                                                              <  1:275069/250‑1 (MQ=255)                        ctccgCCCCCGCTATCCACGGCGACAACTTCCGTAGCTCGTTTTCCTGGTCCGGGAAGTCAACCAACAGCTCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCcagcag                                                                          <  2:359812/245‑1 (MQ=255)                           cgCCCCCGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGACATGCAGTTCACCGCCGCCCCCGCTAACGCGGTCATCCAGAAGCTGCg                                                                     >  2:322012/2‑251 (MQ=255)                             cccccGCTATCCACGTCGACAACTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCggg                                                                   >  2:173744/1‑249 (MQ=255)                                  gctatccacgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCACCCAACAGCGCTTCGCAATCTTGTTGCAGCATATCAGCGCGGTACAGGCATCCTTGCCGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACCCCTGGGCGAGGGAGCTGTGGCCTTTTTCAAAGTCGAAATGCAGTTACACGCCGCCCCAGCTAAAGAGTTCATCCAGCAGCTGCGAGAATGc                                                              >  2:192745/16‑251 (MQ=255)                                     atccacgtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCATATGCGAGAATGCCg                                                            <  1:397664/238‑1 (MQ=255)                                           gtcgacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTg                                                     <  2:126163/245‑1 (MQ=255)                                             cgacaaCTTCCGTAGCTCTTTTGCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTGTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCAGCCAGCAGATGCGAGAATGCCGGAGCCTGAc                                                   <  1:192745/247‑1 (MQ=255)                                              gacaaCTTCCGTAGCTCTTTTTCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGATAATGCCGGAGCCTGAcc                                                  <  1:366736/248‑1 (MQ=255)                                                                ttttCCTGGTCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAg                                  <  1:190959/249‑1 (MQ=255)                                                                  ttccttgtCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCTCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCTGTAAACACCGGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCATATGCGAGAATGCCGGAGCCTGACCTAAGGTCCATACCCAGCTa                               <  2:65649/243‑1 (MQ=255)                                                                        gtCCGGGAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGCGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTgac                        >  1:78756/2‑249 (MQ=255)                                                                         tCCGGAAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGATACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGCGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGACATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCCACATCTGCAAGCATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTActctcacg                       >  2:260821/1‑247 (MQ=255)                                                                            gggAAGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCCGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAAAGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGACGGAGCCGGGCCGACGTTCCATTCCCCGCTACTCTGACGGGc                    >  2:298107/1‑251 (MQ=255)                                                                               aaGTCAACCAACAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCaaa                 <  1:119855/251‑1 (MQ=255)                                                                                      ccaacAGCGCTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGACGGGCAAAGGTTTCg          >  1:147835/1‑251 (MQ=255)                                                                                              gcTTCGCACTCTTGTTGCAGCATATCCGCGCGGTACAGGCAGCCTTGCAGTCGCCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGCGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGCTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGCCGGGCACAGTTTTCGGCGAAGtt  >  2:152508/1‑251 (MQ=255)                                                                                                         |                                                                                                                                                                                                                                                CCCGCATGGGCGGGTAACTACCTTACAGCCCCCGCCATCCATGCCGATAACTCCCGTAGCTCTTTTTCCTGTTCCGGGAAGTCAACCAACAGCGCTTCGCACTCCTGTTGCAGCATATCTGCGCGGTACAGGCAGCCTTGCAGTCGTCCGGCGAGGGCTTCCAGCGGCGCGGGGTTGAGGCTGTCGGTAAACACCTGGGCGCGGGTGATATGGCCTTTTTCAACGTCGAAATGCAGTTCCACGCCGCCCCAGGTAAAGCGTTCATCCAGCAGATGCGAGAATGCCGGAGCCTGACCGAAGTTCCATTCCCAGCTACTCTGGCGGGCAAAGGTTTCGGCGAAGTT  >  NC_000913/4623084‑4623427

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