Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	819,624	C→T	A110A (GCC→GCT)	moaE →	molybdopterin synthase catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	819,624	0	C	T	100.0%	81.9 / NA	24	A110A (GCC→GCT)	moaE	molybdopterin synthase catalytic subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (15/9);  total (15/9)

TGTTGATGAAGCGCGTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAG  >  NC_000913/819492‑819758                                                                                                                                     |                                                                                                                                       tgttgATGAAGCGCGTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGt                                                                                                                                 <  1:613747/140‑1 (MQ=255)   ttgATGAAGCGCGTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGttt                                                                                                                               >  2:496669/1‑140 (MQ=255)   ttgATGAAGCGCGTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGttt                                                                                                                               >  1:86186/1‑140 (MQ=255)             gcgTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCACCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGAtt                                                                                                                     <  1:457522/140‑1 (MQ=255)                         gccgctggGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAAccc                                                                                                         <  2:555921/140‑1 (MQ=255)                            gctggGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCcgcg                                                                                                      <  2:686405/140‑1 (MQ=255)                                   cgcgTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCa                                                                                                    >  1:417028/1‑135 (MQ=255)                                   cgcgTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCa                                                                                                    <  2:417028/135‑1 (MQ=255)                                                          cGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCgg                                                                        >  2:661206/1‑140 (MQ=255)                                                                    tGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCg                                                              <  1:496669/140‑1 (MQ=255)                                                                           gCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGtt                                                       >  1:522999/1‑140 (MQ=255)                                                                           gCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGtt                                                       >  2:140546/1‑140 (MQ=255)                                                                            cGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTg                                                      <  1:140546/140‑1 (MQ=255)                                                                                           tGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCg                                       >  2:167184/1‑140 (MQ=255)                                                                                             tCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGAt                                     >  2:648211/1‑140 (MQ=255)                                                                                                           gcgcATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCGGGGTTGATGCTCGGGAGAGCGCTCAGCAGGCGGCaaa                       >  1:287755/1‑140 (MQ=255)                                                                                                              cATCGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACg                    >  1:138051/1‑140 (MQ=255)                                                                                                                 cGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACg                    >  1:541886/1‑137 (MQ=255)                                                                                                                 cGCAGCAGTGCGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACg                    <  2:541886/137‑1 (MQ=255)                                                                                                                           cGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGtttttt       >  1:490816/1‑140 (MQ=255)                                                                                                                           cGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGtttttt       >  2:653594/1‑140 (MQ=255)                                                                                                                           cGTTTGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGtttttt       >  1:417283/1‑140 (MQ=255)                                                                                                                             tttGAAGCTGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGTGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGt     >  1:186551/1‑140 (MQ=255)                                                                                                                                gAAGCTGGGCAGTTTATTATGGTTTATCTCATAACTCGCCCACCGTTCCGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAg  <  2:211131/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                       TGTTGATGAAGCGCGTAACCGCTGGCCGCTGGGGCGCGTCACTGTGATTCACCGCATCGGGGAATTATGGCCGGGCGATGAAATCGTTTTTGTCGGTGTCACCAGTGCGCATCGCAGCAGTGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAG  >  NC_000913/819492‑819758

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: ‑