Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	864,955	T→C	T59A (ACG→GCG)	moeB ←	molybdopterin‑synthase adenylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	864,955	0	T	C	100.0%	90.7 / NA	26	T59A (ACG→GCG)	moeB	molybdopterin‑synthase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (13/13);  total (13/13)

TAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGAT  >  NC_000913/864820‑865084                                                                                                                                        |                                                                                                                                  tAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAgg                                                                                                                               <  1:11821/140‑1 (MQ=255)        gATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGa                                                                                                                        >  1:310328/1‑140 (MQ=255)        gATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGa                                                                                                                        >  2:218296/1‑140 (MQ=255)             gtgGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCg                                                                                                                   >  2:19383/1‑140 (MQ=255)                 gTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGgcgc                                                                                                               <  2:310328/140‑1 (MQ=255)                  ttGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCt                                                                                                              >  2:298482/1‑140 (MQ=255)                     aTCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACCCCGGCGCTTGc                                                                                                           <  2:79572/140‑1 (MQ=255)                      tCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGGTTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGcc                                                                                                          >  2:353803/1‑140 (MQ=255)                      tCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGcc                                                                                                          >  2:392678/1‑140 (MQ=255)                       ccGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCa                                                                                                         >  2:233691/1‑140 (MQ=255)                          cgTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGAt                                                                                                      >  2:81883/1‑140 (MQ=255)                          cgTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGAt                                                                                                      >  1:112303/1‑140 (MQ=255)                                 gCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGa                                                                                               <  2:135888/140‑1 (MQ=255)                                 gCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGa                                                                                               <  1:358451/140‑1 (MQ=255)                                                 ccACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAgg                                                                               <  2:82307/140‑1 (MQ=255)                                                  cACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGgc                                                                              <  1:19383/140‑1 (MQ=255)                                                      ttCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGgccgcc                                                                          <  1:81883/140‑1 (MQ=255)                                                          gTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAgg                                                                      <  2:220014/140‑1 (MQ=255)                                                            tGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc                                                                    <  1:218296/140‑1 (MQ=255)                                                                                      tgtCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGTCCACTATCAGCACGCGAGAATCTTTCa                                          >  1:223360/1‑140 (MQ=255)                                                                                      tgtCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGc                                                          >  1:77576/1‑124 (MQ=255)                                                                                      tgtCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGc                                                          <  2:77576/124‑1 (MQ=255)                                                                                        tCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAgc                                        >  1:188951/1‑140 (MQ=255)                                                                                                      ttCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGTTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc                          <  2:112303/140‑1 (MQ=255)                                                                                                          agagCGAAACCGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCaaaa                      >  1:199990/1‑140 (MQ=255)                                                                                                                   ccGTGTCGAAGTCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCcg             <  2:199990/140‑1 (MQ=255)                                                                                                                              tCGAGCAGCGCCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGAt  <  1:353803/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  TAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAGCCGCGCAGAATGAT  >  NC_000913/864820‑865084

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