Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	781,828	Δ1 bp	intergenic (+176/‑257)	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,826	0	A	.	100.0%	129.5 / NA	32	intergenic (+174/‑259)	lysQ/nadA	tRNA‑Lys/quinolinate synthase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (14/18);  total (14/18)

CGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCC  >  NC_000913/781707‑781945                                                                                                                        |                                                                                                                        cGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCa                                                                                                    >  2:804382/1‑140 (MQ=17)  gAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAg                                                                                                   <  2:953786/140‑1 (MQ=18)        gcagcaGGTTTGAGTTGAGCACAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATcc                                                                                             >  2:849998/1‑140 (MQ=255)         cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCt                                                                                            <  2:202875/140‑1 (MQ=255)         cagcagGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCt                                                                                            <  2:748477/140‑1 (MQ=255)                        gagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTa                                                                             <  1:804382/140‑1 (MQ=255)                        gagCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTa                                                                             <  2:351802/140‑1 (MQ=255)                         agCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTaa                                                                            >  1:75144/1‑140 (MQ=255)                         agCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTaa                                                                            >  2:438517/1‑140 (MQ=255)                                       acGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCg                                                              <  1:928066/140‑1 (MQ=255)                                         ggAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTc                                                            <  2:356692/140‑1 (MQ=255)                                            gCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCaga                                                         >  2:703153/1‑140 (MQ=255)                                                cTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCtt                                                     >  1:261069/1‑140 (MQ=255)                                                cTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCtt                                                     >  2:1018389/1‑140 (MQ=255)                                                 tCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTa                                                    <  1:288570/140‑1 (MQ=255)                                                 tCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTa                                                    <  2:75144/140‑1 (MQ=255)                                                          cGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaag                                           >  2:695714/1‑140 (MQ=255)                                                             cccGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagcaa                                        <  2:261069/140‑1 (MQ=255)                                                                        agcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagc                                          >  1:236154/1‑127 (MQ=255)                                                                        agcgaagcgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTgcaagc                                          <  2:236154/127‑1 (MQ=255)                                                                               cgaGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAAtt                      <  1:1018389/140‑1 (MQ=255)                                                                                             cGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc        <  1:695714/140‑1 (MQ=255)                                                                                             cGACCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc        <  1:378961/140‑1 (MQ=255)                                                                                               aCCCACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGc      <  1:845831/140‑1 (MQ=255)                                                                                                cccACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTATTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa     >  1:714524/1‑140 (MQ=255)                                                                                                cccACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa     >  2:41995/1‑140 (MQ=255)                                                                                                cccACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa     <  1:849998/140‑1 (MQ=255)                                                                                                cccACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa     >  1:523170/1‑140 (MQ=255)                                                                                                cccACCAATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCa     <  1:438517/140‑1 (MQ=255)                                                                                                   accaATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc        <  2:341377/134‑1 (MQ=255)                                                                                                   accaATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGAc        >  1:341377/1‑134 (MQ=255)                                                                                                   accaATGTAAAAAAGCGCCCT‑AAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATcc  >  2:571338/1‑140 (MQ=255)                                                                                                                        |                                                                                                                        CGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATTCAGGCATCCTCAATTTCACTTTGTAAACCTGATGACATCGTCAGAGCTTACTGTGCAAGCAACTCTATGTCGGTGGAATTAGGCGTAAAATGACGCATCC  >  NC_000913/781707‑781945

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