Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,146,949	T→C	I198V (ATA→GTA)	alkA ←	DNA‑3‑methyladenine glycosylase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,146,949	0	T	C	100.0%	64.3 / NA	20	I198V (ATA→GTA)	alkA	DNA‑3‑methyladenine glycosylase 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (10/10);  total (10/10)

TGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCA  >  NC_000913/2146813‑2147081                                                                                                                                         |                                                                                                                                     tGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCCGCGTTTTAATCGCCGGCTCCAAATCGCCCGGTACTGt                                                                                                                                   >  1:282277/1‑140 (MQ=255) tGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGt                                                                                                                                   <  1:517469/140‑1 (MQ=255)     tAATCAGATAATCATCCAGCAGTAAAACATCCTTCACCTTCCAGCCACTCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCAtt                                                                                                                               <  2:256939/140‑1 (MQ=255)             tAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGt                                                                                                                       >  2:677592/1‑140 (MQ=255)             tAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCACCGTTTTCATCGCGTGCCCCACATCGACCGGAACTGTCATAGGTaactg                                                                                                                       >  2:423655/1‑137 (MQ=255)                        cAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAgcg                                                                                                            <  2:597697/140‑1 (MQ=255)                                 aTCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTg                                                                                                   <  2:615807/140‑1 (MQ=255)                                     ttCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAg                                                                                               >  1:42633/1‑140 (MQ=255)                                     ttCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAg                                                                                               >  2:698707/1‑140 (MQ=255)                                                   aCGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCct                                                                                 <  2:146204/140‑1 (MQ=255)                                                               ataatTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCa                                                                     <  2:778199/140‑1 (MQ=255)                                                                        cGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATAc                                                            >  1:108000/1‑140 (MQ=255)                                                                              gcgcCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACg                                                      <  2:42633/140‑1 (MQ=255)                                                                                                  tttGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTct                                  >  2:76134/1‑140 (MQ=255)                                                                                                           ttttCATCGCCTGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGcc                         >  2:63943/1‑140 (MQ=255)                                                                                                                      tGCTCCACATCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAgg              <  1:76134/140‑1 (MQ=255)                                                                                                                            acaTCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCgg        <  1:63943/140‑1 (MQ=255)                                                                                                                              aTCGCCCGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGa      <  2:108000/140‑1 (MQ=255)                                                                                                                                  cccGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa  >  2:125175/1‑140 (MQ=255)                                                                                                                                  cccGGTACTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCa  >  1:485280/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                     TGTTTAATCAGATAATCATCCGGCAGGAAAACATCTTTCGCCTGCCAGCCACGCAAAGCAAAATAATTCGCCGTCCAGCGCCCGATACCCGGAAAAGTTTGCAGCGTTTTCATCGCCTGCTCCACATCGCCCGGTATTGTCATTGGTAAGGTGCCCTCCAGCGCCGCATTTGCCAGATGAATCAGCGCCTCTGCCCGTTTCAACGGCATACCTAACGCTTTTAATGCCTGCGGGTCTGCTGCTGCCAGCCGCTGAGGCGTCGGGAAGCA  >  NC_000913/2146813‑2147081

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