Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	819,747	G→A	*151* (TAG→TAA)	moaE →	molybdopterin synthase, large subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	819,747	0	G	A	100.0%	194.2 / NA	59	*151* (TAG→TAA)	moaE	molybdopterin synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/30);  total (29/30)

CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819614‑819878                                                                                                                                      |                                                                                                                                    cGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAAttttt                                                                                                                                >  1:1242117/1‑139 (MQ=255)   tttGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTg                                                                                                                              >  2:511162/1‑139 (MQ=255)       aaGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  2:844673/139‑1 (MQ=255)        aGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGc                                                                                                                         <  2:866553/139‑1 (MQ=255)                            aTTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGt                                                                                                     >  1:1074671/1‑139 (MQ=255)                              tATCTCAAAACCCGCGCACCTTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATGAGGCACCGCTGCGTcc                                                                                                   >  2:997569/1‑139 (MQ=255)                              tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGa                                                                                                                     <  1:85240/121‑1 (MQ=255)                              tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGa                                                                                                                     >  2:85240/1‑121 (MQ=255)                              tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                   >  1:283996/1‑139 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  1:102180/114‑1 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          >  2:102180/1‑114 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        <  1:1113466/116‑1 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        >  2:1113466/1‑116 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                 >  2:613472/1‑139 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                 >  2:375672/1‑139 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                 >  2:328834/1‑139 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCAGGAAGGCGAACGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                 >  1:906281/1‑139 (MQ=255)                                tctcAAAACCCGCGCACCGTTCTGGAAGCACGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGATCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCgg                                                                                                 >  2:983070/1‑139 (MQ=255)                                   cAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCaa                                                                                              >  2:1306018/1‑139 (MQ=255)                                      aaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAttt                                                                                           >  1:522514/1‑139 (MQ=255)                                      aaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAttt                                                                                           <  2:1242117/139‑1 (MQ=255)                                        cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  2:188624/106‑1 (MQ=255)                                        cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          >  1:188624/1‑106 (MQ=255)                                        cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         <  2:283996/139‑1 (MQ=255)                                        cccGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         <  2:1031372/139‑1 (MQ=255)                                            cgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         >  2:17842/1‑135 (MQ=255)                                            cgcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         <  1:17842/135‑1 (MQ=255)                                             gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  1:375672/139‑1 (MQ=255)                                             gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  2:522514/139‑1 (MQ=255)                                             gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  1:122294/139‑1 (MQ=255)                                              cACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTc                                                                                   >  2:961873/1‑139 (MQ=255)                                                   ttCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTAc                                                                              >  2:379394/1‑139 (MQ=255)                                                              cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    >  2:1140018/1‑122 (MQ=255)                                                              cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  1:1140018/122‑1 (MQ=255)                                                              cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                   >  1:640527/1‑139 (MQ=255)                                                                aaGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCt                                                                 >  1:365176/1‑139 (MQ=255)                                                                 aGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                                                <  1:983070/139‑1 (MQ=255)                                                                 aGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTg                                                                <  1:379394/139‑1 (MQ=255)                                                                      cGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGCTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCg                                                           <  2:906281/139‑1 (MQ=255)                                                                          ggAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGc                                                       <  1:511162/139‑1 (MQ=255)                                                                              ggCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGt                                                             <  1:1336261/129‑1 (MQ=255)                                                                              ggCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGt                                                             >  2:1336261/1‑129 (MQ=255)                                                                                        gggTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCCGCGTCGATGCACAGCCTCATCAtt                                         <  1:997569/139‑1 (MQ=255)                                                                                              aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         <  1:760587/85‑1 (MQ=255)                                                                                              aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCa                                                                                         >  2:760587/1‑85 (MQ=255)                                                                                               aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                                  <  1:613472/139‑1 (MQ=255)                                                                                                   cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGt                                                                        >  1:525948/1‑97 (MQ=255)                                                                                                   cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGt                                                                        <  2:525948/97‑1 (MQ=255)                                                                                                   cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGcac                                                     >  2:836599/1‑116 (MQ=255)                                                                                                   cGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGcac                                                     <  1:836599/116‑1 (MQ=255)                                                                                                      gagaGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  1:98488/82‑1 (MQ=255)                                                                                                      gagaGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    >  2:98488/1‑82 (MQ=255)                                                                                                        gagCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAg                         <  2:1074671/139‑1 (MQ=255)                                                                                                           cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAtt                                         >  2:242042/1‑120 (MQ=255)                                                                                                           cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAtt                                         <  1:242042/120‑1 (MQ=255)                                                                                                                 gcagGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCtt                <  2:365176/139‑1 (MQ=255)                                                                                                                     gcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg                                   <  1:706434/116‑1 (MQ=255)                                                                                                                     gcggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAg                                   >  2:706434/1‑116 (MQ=255)                                                                                                                               gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  1:328834/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                    CGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819614‑819878

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