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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	864,935	0	C	T	100.0%	187.6 / NA	53	T65T (ACG→ACA)	moeB	molybdopterin‑synthase adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (29/24);  total (29/24)

TCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAG  >  NC_000913/864798‑865070                                                                                                                                          |                                                                                                                                        tCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTgt                                                                                                                                       >  2:691894/1‑140 (MQ=255)   cagcagTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTgtgt                                                                                                                                     <  2:297509/140‑1 (MQ=255)   cagcagTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTgtgt                                                                                                                                     <  1:796680/140‑1 (MQ=255)         tGCATTGACTGGCGTAATCACGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGt                                                                                                                               >  1:1100340/1‑140 (MQ=255)           cATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCg                                                                                                                             >  1:488845/1‑140 (MQ=255)             tagACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGag                                                                                                                           <  2:1292628/138‑1 (MQ=255)                   ggCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGag                                                                                                                           <  2:1253654/134‑1 (MQ=255)                   ggCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGag                                                                                                                           >  1:1253654/1‑134 (MQ=255)                   ggCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGt                                                                                                                     <  2:1376602/140‑1 (MQ=255)                   ggCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGt                                                                                                                     >  1:729749/1‑140 (MQ=255)                    gCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTc                                                                                                                    <  2:1250905/140‑1 (MQ=255)                    gCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATACGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTc                                                                                                                    <  2:580237/140‑1 (MQ=255)                     cGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCa                                                                                                                   >  2:60067/1‑140 (MQ=255)                         aTCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTc                                                                                                                    <  2:1316385/135‑1 (MQ=255)                         aTCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTc                                                                                                                    >  1:1316385/1‑135 (MQ=255)                          tCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTa                                                                                                              >  1:1366132/1‑140 (MQ=255)                            gcgATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTAcc                                                                                                            <  2:293152/140‑1 (MQ=255)                                     gggTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGgc                                                                                                   <  2:1366132/140‑1 (MQ=255)                                          gATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTg                                                                                              <  2:224581/140‑1 (MQ=255)                                          gATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTg                                                                                              >  2:290355/1‑140 (MQ=255)                                           aTCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGc                                                                                             <  2:307491/140‑1 (MQ=255)                                           aTCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGc                                                                                             >  1:150353/1‑140 (MQ=255)                                             ccGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCa                                                                                           <  2:150353/140‑1 (MQ=255)                                                cgTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGAt                                                                                        >  2:18114/1‑140 (MQ=255)                                                cgTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGAt                                                                                        >  2:773511/1‑140 (MQ=255)                                                 gTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATa                                                                                       >  1:466107/1‑140 (MQ=255)                                                      ggCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGTAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCg                                                                                  >  2:887916/1‑140 (MQ=255)                                                      ggCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCg                                                                                  >  1:270449/1‑140 (MQ=255)                                                           cACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCt                                                                             >  1:699613/1‑140 (MQ=255)                                                                       ccACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAgg                                                                 <  2:77979/140‑1 (MQ=255)                                                                       ccACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAgg                                                                 <  2:167962/140‑1 (MQ=255)                                                                                gTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAgg                                                        >  2:883058/1‑140 (MQ=255)                                                                                gTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAgg                                                        <  2:137660/140‑1 (MQ=255)                                                                                  tGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc                                                      <  1:290355/140‑1 (MQ=255)                                                                                  tGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc                                                      <  2:271883/140‑1 (MQ=255)                                                                                  tGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGcc                                                      <  2:1139057/140‑1 (MQ=255)                                                                                      cGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACt                                                  <  2:699613/140‑1 (MQ=255)                                                                                        aCCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTAt                                                <  2:1100340/140‑1 (MQ=255)                                                                                        aCCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTAt                                                >  2:1336685/1‑140 (MQ=255)                                                                                                            tgtCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCa                            >  1:956812/1‑140 (MQ=255)                                                                                                            tgtCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCa                            >  1:399623/1‑140 (MQ=255)                                                                                                             gtCTGGCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAg                           <  2:466107/140‑1 (MQ=255)                                                                                                                 ggCGTTGCAGATTCGAGAGCGAGACTGTGGCGGAGGAGAGCAGCGTCAGGTTACCACCACCGACGCTTGCCACATACAGCGAGGCTGCACACCCGACGCCCCCCCGGCCCCCTAACAACACGCGAGAATCTTTCAGCGcc                       >  2:1321680/1‑140 (MQ=255)                                                                                                                 ggCGTTGCAGATTCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGcc                       >  1:1238552/1‑140 (MQ=255)                                                                                                                      tGCAGATTCGAGAACGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTg                  >  2:1291970/1‑140 (MQ=255)                                                                                                                            ttCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc            >  2:1083050/1‑140 (MQ=255)                                                                                                                            ttCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTc            <  2:1025592/140‑1 (MQ=255)                                                                                                                             tCGAGAGCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCa           >  1:257515/1‑140 (MQ=255)                                                                                                                                agagCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCaaaa        >  1:720981/1‑140 (MQ=255)                                                                                                                                agagCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCaaaa        >  1:1380996/1‑140 (MQ=255)                                                                                                                                  agCGAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATc      >  2:922968/1‑140 (MQ=255)                                                                                                                                     gAAACTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCaaa   >  1:569217/1‑140 (MQ=255)                                                                                                                                      aaaCTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAg  <  2:488845/140‑1 (MQ=255)                                                                                                                                      aaaCTGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAg  >  1:330452/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                        TCCAGCAGTGCATTGACTGGCGTAATCGCGATATGTGGGTTGATCCGCGTCAGGGCGTCACGGGCGGATTCCACCTTCGGTTGCCCGACCGTGGCATCACTGTGCAGTGTCTGGCGTTGCAGATTCGAGAGCGAAACCGTGTCGAAGTCGAGCAGCGTCAGGTTACCGACACCGGCGCTTGCCAGATACTGCGAGGCTGCACAGCCGAGGCCGCCCAGGCCCACTATCAGCACGCGAGAATCTTTCAGCGCCTCCTGGCCGTCAAAATCAAAG  >  NC_000913/864798‑865070

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