Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	865,504	A→G	G287G (GGT→GGC)	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,504	0	A	G	100.0%	173.3 / NA	49	G287G (GGT→GGC)	moeA	molybdopterin molybdenumtransferase, molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/26);  total (23/26)

GCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGG  >  NC_000913/865369‑865639                                                                                                                                        |                                                                                                                                        gcgcTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcg                                                                                                                                      <  1:672752/139‑1 (MQ=255)    cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaac                                                                                                                                   >  2:69697/1‑139 (MQ=255)      ggCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgc                                                                                                                                 <  1:69697/139‑1 (MQ=255)       gcgcgcGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                <  1:479393/139‑1 (MQ=255)             gggCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGt                                                                                                                          >  1:1372267/1‑139 (MQ=255)                cAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTa                                                                                                                       <  2:568499/139‑1 (MQ=255)                                        tAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAg                                                                                               <  2:384825/139‑1 (MQ=255)                                             ttGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                <  1:431146/101‑1 (MQ=255)                                             ttGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                >  2:431146/1‑101 (MQ=255)                                               gCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      >  2:372975/1‑109 (MQ=255)                                               gCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      <  1:372975/109‑1 (MQ=255)                                               gCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGa                                                                                        <  2:1372267/139‑1 (MQ=255)                                                ccAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAt                                                                                       <  1:1587706/139‑1 (MQ=255)                                                 cAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc                                                                                      <  1:1475464/139‑1 (MQ=255)                                                     aaaGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTcc                                                                                   >  2:1206104/1‑138 (MQ=255)                                                     aaaGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTcc                                                                                   <  1:1206104/138‑1 (MQ=255)                                                              aCCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTa                                                                                                                       <  1:515902/93‑1 (MQ=255)                                                              aCCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTa                                                                                                                       >  2:515902/1‑93 (MQ=255)                                                                  gTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                           <  2:432801/117‑1 (MQ=255)                                                                  gTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                           >  1:432801/1‑117 (MQ=255)                                                                   ttGATAGAAGGTCAGCGTCGTTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                <  1:1322185/79‑1 (MQ=255)                                                                   ttGATAGAAGGTCAGCGTCGTTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                >  2:1322185/1‑79 (MQ=255)                                                                   ttGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcg                                                                                                                                      <  1:632454/73‑1 (MQ=255)                                                                   ttGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcg                                                                                                                                      >  2:632454/1‑73 (MQ=255)                                                                       tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGTcgaacgcg                                                                                                                                <  1:1075342/75‑1 (MQ=255)                                                                       tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                                >  2:1075342/1‑75 (MQ=255)                                                                       tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCg                                                                                                          >  2:1597534/1‑97 (MQ=255)                                                                       tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCg                                                                                                          <  1:1597534/97‑1 (MQ=255)                                                                          aaGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGttt                                                             <  2:386098/139‑1 (MQ=255)                                                                            ggTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGc                                                                                                         >  1:1552957/1‑93 (MQ=255)                                                                            ggTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGc                                                                                                         <  2:1552957/93‑1 (MQ=255)                                                                            ggTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGaa                                                                   <  1:49070/131‑1 (MQ=255)                                                                            ggTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGaa                                                                   >  2:49070/1‑131 (MQ=255)                                                                             gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      <  2:514470/79‑1 (MQ=255)                                                                             gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      >  1:514470/1‑79 (MQ=255)                                                                             gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGa                                                                                              <  1:665224/103‑1 (MQ=255)                                                                             gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTgg                                                          >  2:1613252/1‑139 (MQ=255)                                                                             gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGGTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGa                                                                                              >  2:665224/1‑103 (MQ=255)                                                                                aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGTTCTTCAAGAATCGTTTTGgtgt                                                       >  2:1586901/1‑139 (MQ=255)                                                                                       cTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGc                                                >  1:989549/1‑139 (MQ=255)                                                                                        tGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGATTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGcc                                               <  2:1187581/139‑1 (MQ=255)                                                                                           aaCCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCa                                            <  2:989549/139‑1 (MQ=255)                                                                                                 gTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      >  1:902019/1‑59 (MQ=255)                                                                                                 gTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                      <  2:902019/59‑1 (MQ=255)                                                                                                                          cTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATcacca             >  2:72877/1‑139 (MQ=255)                                                                                                                              tGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCt                                              >  2:393942/1‑102 (MQ=255)                                                                                                                              tGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCt                                              <  1:393942/102‑1 (MQ=255)                                                                                                                                 tGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGc      >  1:156668/1‑139 (MQ=255)                                                                                                                                     tttGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTgg  >  1:814305/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                        GCGCTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCCTGG  >  NC_000913/865369‑865639

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