Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	865,008	C→T	G41D (GGC→GAC)	moeB ←	molybdopterin synthase sulfurylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	865,008	0	C	T	100.0%	91.5 / NA	27	G41D (GGC→GAC)	moeB	molybdopterin synthase sulfurylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/13);  total (14/13)

TCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGG  >  NC_000913/864894‑865135                                                                                                                   |                                                                                                                                tCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGGGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCAcgc                                                                                                         <  1:613599/139‑1 (MQ=255)   aCTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGcc                                                                                                                       <  2:61499/123‑1 (MQ=255)   aCTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGcc                                                                                                                       >  1:61499/1‑123 (MQ=255)      gtgCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGaa                                                                                                    <  2:230495/139‑1 (MQ=255)       tgCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAAt                                                                                                   >  1:223516/1‑139 (MQ=255)         cAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCt                                                                                                 >  1:8750/1‑139 (MQ=255)            tgtCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTc                                                                                              >  1:183597/1‑139 (MQ=255)              tCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAg                                                                                            <  2:183597/139‑1 (MQ=255)               cTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAgc                                                                                           <  1:781236/139‑1 (MQ=255)                                  agCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                                                                             >  1:595012/1‑134 (MQ=255)                                  agCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                                                                             <  2:595012/134‑1 (MQ=255)                                                                            ccGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGccc                                                                                                                            >  1:439261/1‑45 (MQ=255)                                                                            ccGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGccc                                                                                                                            <  2:439261/45‑1 (MQ=255)                                                                                cgcTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCgcagc                                  >  1:1000251/1‑131 (MQ=255)                                                                                cgcTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCgcagc                                  <  2:1000251/131‑1 (MQ=255)                                                                                cgcTTGCCAGATACTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTg                          >  1:402002/1‑139 (MQ=255)                                                                                            aCTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTc              <  1:569576/139‑1 (MQ=255)                                                                                             cTGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATGAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGAGTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTcc             <  1:754473/139‑1 (MQ=255)                                                                                              tGCGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCg            >  2:415700/1‑139 (MQ=255)                                                                                                cGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                                                                             <  2:986790/72‑1 (MQ=255)                                                                                                cGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                                                                             >  1:986790/1‑72 (MQ=255)                                                                                                cGAGGCTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGcc          >  1:351319/1‑139 (MQ=255)                                                                                                     cTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGGCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTAc     >  2:791306/1‑139 (MQ=255)                                                                                                     cTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGcc                                                <  2:285165/96‑1 (MQ=255)                                                                                                     cTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGcc                                                >  1:285165/1‑96 (MQ=255)                                                                                                     cTGCACAGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTAc     >  2:214258/1‑139 (MQ=255)                                                                                                        cacaGCCGAGGTCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAgg  <  2:402002/139‑1 (MQ=255)                                                                                                                   |                                                                                                                                TCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGATTTGCCGGTTGTAGCGCAGCATCTCCTGATCGCTGAGTTCCGCCATTACAGG  >  NC_000913/864894‑865135

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