Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	781,812	T→C	intergenic (+160/‑273)	lysQ → / → nadA	tRNA‑Lys/quinolinate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	781,812	0	T	C	100.0%	36.7 / NA	13	intergenic (+160/‑273)	lysQ/nadA	tRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (7/6);  total (7/6)

TTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATG  >  NC_000913/781723‑781937                                                                                          |                                                                                                                              ttgagttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCt                                                                                    >  2:43352/1‑133 (MQ=255)      ttgagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAACGCCCTAAGGGCGCTTATTTGCTATACAGGCATCCTCAAtt                                                                               >  2:149367/1‑133 (MQ=14)         agCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGTTATTCAGGCATCCTCAATTTCa                                                                            >  1:350046/1‑133 (MQ=255)         agCGCAGCCAAACAACGGGGCCTCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCa                                                                            <  1:174019/133‑1 (MQ=255)                               gCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCGCGACCCACCAATGCAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCg                                                      <  2:421530/133‑1 (MQ=255)                               gCTCGCGGAGACGGCCCGAAGGCCGAGCGACGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCg                                                      <  1:194705/133‑1 (MQ=255)                                cTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCTTCAATTTCACTTTGTACACCTGATGACATCGt                                                     >  2:199198/1‑133 (MQ=255)                                 tCGCGGCGACGGCCCGAAGGGCGAGCGAGGCGGGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTc                                                    <  2:205056/133‑1 (MQ=255)                                                  aGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag                                   <  1:98183/133‑1 (MQ=255)                                                      cgagcgaAGCGAGTCATCCTGCACGACCCACCAATGCAAAAAAGCTCACTACAGGCGCTTTTTTGCTATTCAGGGATCCTCACTTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAAc                               >  2:381992/1‑133 (MQ=25)                                                                      tCCTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAt               <  1:432426/133‑1 (MQ=255)                                                                       ccTGCACGACCCACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAtt              >  1:304775/1‑133 (MQ=255)                                                                            aCGACCCACCCATGCAAAAAAGCGCCCTAAAGTCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCg         <  2:222479/133‑1 (MQ=255)                                                                                cccACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTaaa     >  2:513106/1‑133 (MQ=255)                                                                                cccACCAATGCAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCA‑‑‑CTATGTCGGTGGAATTAGGCGTAAAATg  >  2:3973/1‑133 (MQ=38)                                                                                          |                                                                                                                              TTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATG  >  NC_000913/781723‑781937

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