Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,982,601	A→G	intergenic (+16/‑5)	hisR → / → leuT	tRNA‑His/tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,982,601	0	A	G	100.0%	77.3 / NA	22	intergenic (+16/‑5)	hisR/leuT	tRNA‑His/tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/8);  total (14/8)

ATACCGCGTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCG  >  NC_000913/3982485‑3982736                                                                                                                     |                                                                                                                                        aTACCGCGTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTgg                                                                                                                  >  1:401678/1‑140 (MQ=255)   aCCGCGTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTa                                                                                                                >  2:154189/1‑140 (MQ=255)    ccGCGTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTAGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAg                                                                                                               <  2:418452/140‑1 (MQ=255)      gcgTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGAc                                                                                                             >  2:177731/1‑140 (MQ=255)                  ttGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCa                                                                                                 >  1:46298/1‑140 (MQ=255)                          gtggCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAg                                                                                         >  1:486632/1‑140 (MQ=255)                             gCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAgtgt                                                                                      >  2:396726/1‑140 (MQ=255)                                  aGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTa                                                                                 >  1:249564/1‑140 (MQ=255)                                       aGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGAc                                                                            >  1:169455/1‑140 (MQ=255)                                                 gCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTc                                                                  >  1:177508/1‑140 (MQ=255)                                                 gCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTc                                                                  <  1:407141/140‑1 (MQ=255)                                                  cccTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCa                                                                 <  2:177508/140‑1 (MQ=255)                                                        aTTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTccc                                                           <  2:486632/140‑1 (MQ=255)                                                                 ccAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGc                                                  >  2:60124/1‑140 (MQ=255)                                                                      tGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTTGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGcaccac                                             >  1:225363/1‑140 (MQ=255)                                                                             ggTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTa                                      <  1:60124/140‑1 (MQ=255)                                                                             ggTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTa                                      <  1:257213/140‑1 (MQ=255)                                                                                                   ccATTATTAGAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCaa                >  1:298487/1‑140 (MQ=255)                                                                                                            gAAGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTa       <  2:320892/140‑1 (MQ=255)                                                                                                              aGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa          <  2:335478/135‑1 (MQ=255)                                                                                                              aGTTGTGGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCa          >  1:335478/1‑135 (MQ=255)                                                                                                                 tgtgGCAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCg  >  2:436248/1‑140 (MQ=255)                                                                                                                     |                                                                                                                                        ATACCGCGTAACAAGATTTGTAGTGGTGGCTATAGCTCAGTTGGTAGAGCCCTGGATTGTGATTCCAGTTGTCGTGGGTTCGAATCCCATTAGCCACCCCATTATTAGAAGTTGTGACAATGCGAAGGTGGCGGAATTGGTAGACGCGCTAGCTTCAGGTGTTAGTGTCCTTACGGACGTGGGGGTTCAAGTCCCCCCCCTCGCACCACGACTTTAAAGAATTGAACTAAAAATTCAAAAAGCAGTATTTCG  >  NC_000913/3982485‑3982736

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