Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	865,423	A→G	P314P (CCT→CCC)	moeA ←	molybdopterin molybdenumtransferase, molybdopterin biosynthesis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	865,423	0	A	G	100.0%	147.1 / NA	39	P314P (CCT→CCC)	moeA	molybdopterin molybdenumtransferase, molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/21);  total (18/21)

CCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATC  >  NC_000913/865289‑865555                                                                                                                                       |                                                                                                                                     ccATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCt                                                                                                                                  >  2:326251/1‑139 (MQ=255)     cGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTAc                                                                                                                              <  1:737109/139‑1 (MQ=255)     cGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTAc                                                                                                                              <  2:389342/139‑1 (MQ=255)     cGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTAc                                                                                                                              <  1:791384/139‑1 (MQ=255)       gcgttgcgttGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCa                                                                                                                            >  1:532535/1‑139 (MQ=255)                        gccgcGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGc                                                                                                           <  1:675294/139‑1 (MQ=255)                          cgcgCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAggg                                                                                                                                    >  1:123914/1‑112 (MQ=255)                                         gACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGtt                                                                                          >  1:230680/1‑139 (MQ=255)                                             tCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGccc                                                                                      >  2:835978/1‑139 (MQ=255)                                              ccTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCg                                                                                     >  2:648333/1‑139 (MQ=255)                                                 ggCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCa                                                                                  <  1:58338/139‑1 (MQ=255)                                                  gCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAg                                                                                 <  2:532535/139‑1 (MQ=255)                                                       tttttCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTTTACCAGTTGATAGAAGGTCAGCGTCGc                                                                                                      <  1:703073/113‑1 (MQ=255)                                                       tttttCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCACTTGATAGAAGGTCAGCGTCGc                                                                                                      >  2:703073/1‑113 (MQ=255)                                                           tCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGaa                                                                        <  2:230680/139‑1 (MQ=255)                                                             aGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCg                                                                             >  2:857086/1‑132 (MQ=255)                                                             aGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCg                                                                             <  1:857086/132‑1 (MQ=255)                                                                cGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGc                                                                   <  1:648333/139‑1 (MQ=255)                                                                       cTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCt                                                            <  2:846213/139‑1 (MQ=255)                                                                         gtgCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGa                                                          >  1:218714/1‑139 (MQ=255)                                                                                cgcgCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTAc                                                                                                                              <  1:845438/64‑1 (MQ=255)                                                                                cgcgCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTAc                                                                                                                              >  2:845438/1‑64 (MQ=255)                                                                                cgcgCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTAcc                                                   >  2:740533/1‑139 (MQ=255)                                                                                  cgcTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCg                                                                                                          <  1:829433/82‑1 (MQ=255)                                                                                  cgcTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCg                                                                                                          >  2:829433/1‑82 (MQ=255)                                                                                    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGc                                                                   <  1:1124437/119‑1 (MQ=255)                                                                                    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGc                                                                   >  2:1124437/1‑119 (MQ=255)                                                                                    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaac                                               >  2:365923/1‑139 (MQ=255)                                                                                    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaac                                               >  2:739963/1‑139 (MQ=255)                                                                                    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaac                                               >  1:1029304/1‑139 (MQ=255)                                                                                       gcgcgcGGGCAGGCTGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaacgcg                                            <  2:1056148/139‑1 (MQ=255)                                                                                       gcgcgcGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaacgcg                                            <  1:740533/139‑1 (MQ=255)                                                                                       gcgcgcGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACcgaacgcg                                            <  1:42636/139‑1 (MQ=255)                                                                                                 aGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTAc                                  <  2:102745/139‑1 (MQ=255)                                                                                                      gCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGt                             <  2:74100/139‑1 (MQ=255)                                                                                                                         aGTTTTGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGa          >  1:601333/1‑139 (MQ=255)                                                                                                                              tGCCAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCg     >  1:281024/1‑139 (MQ=255)                                                                                                                                ccAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAt   >  2:583457/1‑139 (MQ=255)                                                                                                                                 cAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc  <  2:601333/139‑1 (MQ=255)                                                                                                                                 cAGCAAGGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc  <  2:1029304/139‑1 (MQ=255)                                                                                                                                       |                                                                                                                                     CCATCGGCGTTGCGTTGCAGCACGCCGCGCTGGAAATCAAGACGTCCTGGCGTTTTTTTCAGGCGGGACGCTGTGCGTACGCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATC  >  NC_000913/865289‑865555

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