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%	163.1 / NA	51	*151* (TAG→TAA)	moaE	molybdopterin synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (24/27);  total (24/27)

TGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819612‑819878                                                                                                                                        |                                                                                                                                    tGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAAttt                                                                                                                                  >  2:69241/1‑139 (MQ=255)     tttGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTg                                                                                                                              >  2:1069309/1‑139 (MQ=255)         aaGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  2:214958/139‑1 (MQ=255)         aaGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  2:11997/139‑1 (MQ=255)         aaGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAg                                                                                                                          <  1:1111702/139‑1 (MQ=255)           gCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        <  2:630775/139‑1 (MQ=255)                gCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATa                                                                                                                   <  1:43519/139‑1 (MQ=255)                              aTTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGt                                                                                                     >  1:549042/1‑139 (MQ=255)                                tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                   >  1:111162/1‑139 (MQ=255)                                tATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTcc                                                                                                   >  1:905769/1‑139 (MQ=255)                                       aaaCCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAAtt                                                                                            >  2:495906/1‑139 (MQ=255)                                            cgcgCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCAcc                                                                                       >  2:1469665/1‑139 (MQ=255)                                               gcACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  2:476377/139‑1 (MQ=255)                                                      tCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACa                                                                             <  1:1069309/139‑1 (MQ=255)                                                          gAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACg                                                                         >  2:1529359/1‑139 (MQ=255)                                                               gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCTGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTAc                                                                    >  2:240509/1‑139 (MQ=255)                                                               gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATa                                                                                                                   <  1:1059865/92‑1 (MQ=255)                                                               gcgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATa                                                                                                                   >  2:1059865/1‑92 (MQ=255)                                                                cgAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACg                                                                   >  1:260727/1‑139 (MQ=255)                                                                      cACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTCCACCGCTCAGTACATACGTTTACGCTGCGt                                                             <  1:495906/139‑1 (MQ=255)                                                                             gAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGca                                                      >  1:978984/1‑139 (MQ=255)                                                                              aaGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACa                                                                             >  1:385187/1‑115 (MQ=255)                                                                              aaGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACa                                                                             <  2:385187/115‑1 (MQ=255)                                                                               aGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGcaca                                                    >  1:954069/1‑139 (MQ=255)                                                                                     ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca                                              <  2:330546/139‑1 (MQ=255)                                                                                     ccGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTAGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCtca                                              <  2:905769/139‑1 (MQ=255)                                                                                       gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTa                                                                                                                           >  2:256340/1‑60 (MQ=255)                                                                                       gCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTa                                                                                                                           <  1:256340/60‑1 (MQ=255)                                                                                          gggTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCAtt                                         <  1:1529359/139‑1 (MQ=255)                                                                                               gAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCa                                    <  1:1469665/139‑1 (MQ=255)                                                                                                aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGa                                                                                                                     <  1:824445/57‑1 (MQ=255)                                                                                                aaGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGa                                                                                                                     >  2:824445/1‑57 (MQ=255)                                                                                                 aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                                  <  2:111162/139‑1 (MQ=255)                                                                                                 aGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGt                                  <  2:260727/139‑1 (MQ=255)                                                                                                   cTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTAt                                >  1:15638/1‑139 (MQ=255)                                                                                                           agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAgg                                                                                                                >  2:318233/1‑51 (MQ=255)                                                                                                           agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAgg                                                                                                                <  1:318233/51‑1 (MQ=255)                                                                                                           agCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGa                        >  2:239458/1‑139 (MQ=255)                                                                                                            gCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGat                       <  1:1229593/139‑1 (MQ=255)                                                                                                            gCGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGat                       <  2:954069/139‑1 (MQ=255)                                                                                                             cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    >  1:1254608/1‑77 (MQ=255)                                                                                                             cGATCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCt                                                                                    <  2:1254608/77‑1 (MQ=255)                                                                                                               aTCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        >  1:1237932/1‑39 (MQ=37)                                                                                                               aTCAGCAGGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        <  2:1237932/39‑1 (MQ=255)                                                                                                                  agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        <  2:837620/36‑1 (MQ=255)                                                                                                                  agcagGCGGCAAAACGCTGGTAATTTTTTGTTAGcc                                                                                                                        >  1:837620/1‑36 (MQ=37)                                                                                                                         ggcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCt          >  2:1013958/1‑139 (MQ=255)                                                                                                                          gcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat              >  1:1518796/1‑134 (MQ=255)                                                                                                                          gcAAAACGCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTtat              <  2:1518796/134‑1 (MQ=255)                                                                                                                                 gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  2:15638/139‑1 (MQ=255)                                                                                                                                 gCTGGTAATTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGgaga  <  1:1013958/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                    TGCGTTTGAAGCCGGGCAGTTTATTATGGATTATCTCAAAACCCGCGCACCGTTCTGGAAGCGCGAAGCCACGCCGGAAGGCGACCGCTGGGTTGAAGCTCGGGAGAGCGATCAGCAGGCGGCAAAACGCTGGTAGTTTTTTGTTAGCCGGATAAGGCACCGCTGCGTCCGGCAATTTCACCGCTCAGTACATACGTTTACGCTGCGTCGATGCACAGCCTCATCATTTTGCAGTATCCTTAAGATATTCCTTATATCTTCAGGAGA  >  NC_000913/819612‑819878

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