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%	222.2 / NA	61	G287G (GGT→GGC)	moeA	molybdopterin molybdenumtransferase, molybdopterin biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/31);  total (30/31)

CTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCC  >  NC_000913/865372‑865636                                                                                                                                     |                                                                                                                                     cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaac                                                                                                                                >  1:922862/1‑139 (MQ=255) cTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCAAGCAACGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaac                                                                                                                                >  1:1678526/1‑139 (MQ=255)   ggCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgc                                                                                                                              >  2:1429697/1‑139 (MQ=255)    gcgcgcGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcg                                                                                                                             <  1:952723/139‑1 (MQ=255)        gcgGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaacgcgaacg                                                                                                                         <  1:1429697/139‑1 (MQ=255)              aGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTAc                                                                                                                   >  2:694513/1‑139 (MQ=255)                      ggCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTa                                                                                                           >  2:944058/1‑139 (MQ=255)                        cgctgtTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGGTGAGTTTGCCGAACGCGAACGGGCTACCTGGTTTAAt                                                                                                         <  1:1189811/136‑1 (MQ=255)                          gtgtTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCg                                                                                                       >  2:764489/1‑139 (MQ=255)                               aCCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGc                                                                                                  >  1:1313662/1‑139 (MQ=255)                                  gCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTc                                                                                               >  1:1543386/1‑139 (MQ=255)                                        ttttGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAg                                                                                         <  1:944058/139‑1 (MQ=255)                                            gCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACC‑GGTTTAATCGCCAGCTTCCAGAAGGCGAt                                                                                    <  1:694513/139‑1 (MQ=255)                                             ccAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTa                                                                                                                    >  2:1272217/1‑107 (MQ=255)                                             ccAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGACGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTa                                                                                                                    <  1:1272217/107‑1 (MQ=255)                                              cAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc                                                                                   <  1:735948/139‑1 (MQ=255)                                                 cAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                        >  2:413125/1‑131 (MQ=255)                                                 cAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                        <  1:413125/131‑1 (MQ=255)                                                 cAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAAGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCAcc                                                                                >  1:1187758/1‑139 (MQ=255)                                                   aaGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTcccc                                                                              <  2:240409/139‑1 (MQ=255)                                                          tACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGctct                                                                       >  1:547834/1‑139 (MQ=255)                                                                ttGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGc                                                                                                      <  2:223678/102‑1 (MQ=255)                                                                ttGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGc                                                                                                      >  1:223678/1‑102 (MQ=255)                                                                ttGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGa                                                                 >  2:1658478/1‑139 (MQ=255)                                                                 tGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGaa                                                                <  1:1459345/139‑1 (MQ=255)                                                                  gATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAAt                                                               >  1:308420/1‑139 (MQ=255)                                                                    tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                        >  2:13707/1‑112 (MQ=255)                                                                    tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAgg                                                                                        <  1:13707/112‑1 (MQ=255)                                                                    tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGa                                                                                     >  2:1234435/1‑115 (MQ=255)                                                                    tAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGa                                                                                     <  1:1234435/115‑1 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc                                                                                   <  1:698442/111‑1 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGAtc                                                                                   >  2:698442/1‑111 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTgg                                                       >  1:854003/1‑139 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTgg                                                       >  2:51894/1‑139 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTgg                                                       >  1:642734/1‑139 (MQ=255)                                                                          gTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTgg                                                       >  2:1184821/1‑139 (MQ=255)                                                                             aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaac                                                                                                                                >  2:592535/1‑63 (MQ=255)                                                                             aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCcgaac                                                                                                                                <  1:592535/63‑1 (MQ=255)                                                                             aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAg                                                                           >  2:336817/1‑116 (MQ=255)                                                                             aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAg                                                                           <  1:336817/116‑1 (MQ=255)                                                                             aGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGgtgt                                                    >  1:814775/1‑139 (MQ=255)                                                                               cGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTaa                                                  <  1:1523049/139‑1 (MQ=255)                                                                                  cGCTGAAACCGGGATGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATcc                                               >  1:1206668/1‑139 (MQ=255)                                                                                    cTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGc                                             >  2:1183486/1‑139 (MQ=255)                                                                                    cTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGc                                             >  1:71975/1‑139 (MQ=255)                                                                                        aaCCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGACGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCa                                         <  2:1678526/139‑1 (MQ=255)                                                                                          ccGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCAcc                                       <  2:1313662/139‑1 (MQ=255)                                                                                                   ccGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAac                              <  2:854003/139‑1 (MQ=255)                                                                                                     ggCAGGCCGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAacac                            <  1:1390144/139‑1 (MQ=255)                                                                                                           ccGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCgg                      <  1:1183486/139‑1 (MQ=255)                                                                                                            cGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTaa                                                  >  2:1528609/1‑110 (MQ=255)                                                                                                            cGCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTaa                                                  <  1:1528609/110‑1 (MQ=255)                                                                                                             gCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGaa                    <  2:1118536/139‑1 (MQ=255)                                                                                                             gCAGAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGaa                    <  2:642734/139‑1 (MQ=255)                                                                                                                gAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTg                 <  1:764489/139‑1 (MQ=255)                                                                                                                gAACCAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTg                 <  2:1187758/139‑1 (MQ=255)                                                                                                                   ccAGCTATTGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATc              <  2:349995/139‑1 (MQ=255)                                                                                                                          ttGCTGAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACAt       <  2:1543386/139‑1 (MQ=255)                                                                                                                               gAGTTTGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGcc  <  2:308420/139‑1 (MQ=255)                                                                                                                                  tttGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGc                                                                                       >  2:1387308/1‑51 (MQ=255)                                                                                                                                  tttGCCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCt                                           <  1:101637/95‑1 (MQ=255)                                                                                                                                     |                                                                                                                                     CTGGCGCGCGGGCAGGCCGCTGGCGGTGTTACCGCTTAGTTTTGCCAGCAAAGGCTGTACCAGTTGATAGAAGGTCAGCGTCGCTGAAACCGGGTTGCCCGGCAGGCCGCAGAACCAGCTATTGCTGAGTTTACCGAACGCGAACGGTTTACCTGGTTTAATCGCCAGCTTCCAGAAGGCGATCTCCCCCAGCTCTTCAAGAATCGTTTTGGTGTAATCCGCCTCACCCACTGAAACACCGCCGGAACTGATCACCACATCCGCC  >  NC_000913/865372‑865636

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