Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,891,299	C→A	intergenic (‑32/+26)	tsaB ← / ← yoaA	N(6)‑L‑threonylcarbamoyladenine synthase, TsaB subunit/putative 5' to 3' DNA helicase implicated in DNA repair

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,891,299	0	C	A	100.0%	104.5 / NA	30	intergenic (‑32/+26)	tsaB/yoaA	N(6)‑L‑threonylcarbamoyladenine synthase, TsaB subunit/putative 5' to 3' DNA helicase implicated in DNA repair
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (13/17);  total (13/17)

AGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGC  >  NC_000913/1891163‑1891427                                                                                                                                         |                                                                                                                                 aGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGca                                                                                                                               <  2:68180/140‑1 (MQ=255) aGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGca                                                                                                                               <  1:140786/140‑1 (MQ=255) aGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGca                                                                                                                               <  1:319003/140‑1 (MQ=255) aGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGca                                                                                                                               <  2:110883/140‑1 (MQ=255) aGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGca                                                                                                                               <  1:386238/140‑1 (MQ=255)   tATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGcatc                                                                                                                             <  1:163197/140‑1 (MQ=255)           cGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATcataccata                                                                                                                     >  1:447196/1‑140 (MQ=255)             aGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACt                                                                                                                   >  2:131870/1‑140 (MQ=255)                gCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCg                                                                                                                >  1:344090/1‑140 (MQ=255)                                gCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTgga                                                                                                >  1:135808/1‑140 (MQ=255)                                   ttGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTggagga                                                                                             <  2:135808/140‑1 (MQ=255)                                      aCAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATgg                                                                                          <  2:67325/140‑1 (MQ=255)                                         gTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTAt                                                                                       <  2:27807/140‑1 (MQ=255)                                          tACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATc                                                                                      >  1:348524/1‑140 (MQ=255)                                           aCCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCg                                                                                     <  2:415098/140‑1 (MQ=255)                                               tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAg                                                                                 >  1:46078/1‑140 (MQ=255)                                               tcgtTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAg                                                                                 <  1:104449/140‑1 (MQ=255)                                                         ccgCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAAcc                                                                       >  1:138827/3‑140 (MQ=255)                                                                 gagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGc                                                               <  1:128024/140‑1 (MQ=255)                                                                 gagCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGc                                                               <  1:131870/140‑1 (MQ=255)                                                                      ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGt                                                          >  1:44726/1‑140 (MQ=255)                                                                      ggCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGt                                                          <  2:447196/140‑1 (MQ=255)                                                                               cgcgGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGc                                                 >  1:429805/1‑140 (MQ=255)                                                                                 cgGTATCGATAGCCAGAATTCACATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGt                                               >  2:92186/1‑140 (MQ=255)                                                                                                gAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGa                                <  2:121128/140‑1 (MQ=255)                                                                                                                tCGTGCTTAGATCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCACACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTcg                <  2:56548/140‑1 (MQ=255)                                                                                                                          aTCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAgg      >  2:253695/1‑140 (MQ=255)                                                                                                                          aTCAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCACACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAgg      >  1:460120/1‑140 (MQ=255)                                                                                                                            cAATAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGAc    <  2:348524/140‑1 (MQ=255)                                                                                                                              aTAAAAAGGCGAGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGc  >  2:347325/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                 AGTATGTTCACGAGGGCAAAGCTCAAAATGAGCGTTGACAGTACCGTCGTTCCACAGGGCGACAGAGCAGGCCTCTGTCGCGGTATCGATAGCCAGAATTCGCATGGGTCTTCGTGCTTAGATCAATAAAAAGGCGCGCATCATACCATACTCCGTAACAAATTACCTGGAGGATGGTATCGCAAGGAAACGAACCGCACGGGCAATGTCACGGGTGCGTGGCGCGGGCGGCAGACTGGCGAGAAACGTCGCGCCGTAAGGACGC  >  NC_000913/1891163‑1891427

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