Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,059,014	G→A	V311M (GTG→ATG)	torC →	trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,059,014	0	G	A	100.0%	209.6 / NA	64	V311M (GTG→ATG)	torC	trimethylamine N‑oxide (TMAO) reductase I, cytochrome c‑type subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/33);  total (31/33)

AACGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGGTGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCACGCTCAACGGCATGAT  >  NC_000913/1058880‑1059144                                                                                                                                       |                                                                                                                                   aaCGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                <  2:117668/139‑1 (MQ=255) aaCGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                <  1:590808/139‑1 (MQ=255)  aCGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                <  1:566852/138‑1 (MQ=255)  aCGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                >  2:566852/1‑138 (MQ=255)   cGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACg                                                                                                                              >  2:1477008/1‑139 (MQ=255)      aTCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGata                                                                                                                           <  2:1419637/139‑1 (MQ=255)              tgCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCg                                                                                                                   >  2:980627/1‑139 (MQ=255)                              gATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCgg                                                                                                   <  1:46013/139‑1 (MQ=255)                               aTGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGa                                                                                                  <  1:1009606/139‑1 (MQ=255)                                 gtgCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGAtt                                                                                                <  2:371730/139‑1 (MQ=255)                                       cagcagGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGAt                                                                                                                 <  1:312316/116‑1 (MQ=255)                                       cagcagGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGAt                                                                                                                 >  2:312316/1‑116 (MQ=255)                                         gcagGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCaaa                                                                                                                      >  1:801939/1‑109 (MQ=255)                                         gcagGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCaaa                                                                                                                      <  2:801939/109‑1 (MQ=255)                                              tAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACgg                                                                                   <  1:980627/139‑1 (MQ=255)                                                    cGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                >  1:582938/1‑88 (MQ=255)                                                    cGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGaa                                                                                                                                <  2:582938/88‑1 (MQ=255)                                                                  accGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACg                                                               <  1:1570977/139‑1 (MQ=255)                                                                      ttGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGTCACGgcgc                                                           >  1:241823/1‑139 (MQ=255)                                                                              acCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAGAGGCGACATGATGAACGATATCaaa                                                                                                                      >  1:423569/1‑72 (MQ=255)                                                                              acCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAGAGGCGACATGATGAACGATATCaaa                                                                                                                      <  2:423569/72‑1 (MQ=255)                                                                                   tACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGAt                                                                                                                 <  2:229396/72‑1 (MQ=255)                                                                                   tACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGAt                                                                                                                 >  1:229396/1‑72 (MQ=255)                                                                                       gAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGcc                                              <  1:1388305/135‑1 (MQ=255)                                                                                       gAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGcc                                              >  2:1388305/1‑135 (MQ=255)                                                                                        aGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCaaa                                                                                                                      <  1:1405176/62‑1 (MQ=255)                                                                                        aGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCaaa                                                                                                                      >  2:1405176/1‑62 (MQ=255)                                                                                             aGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    >  1:112070/1‑139 (MQ=255)                                                                                             aGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    >  2:934590/1‑139 (MQ=255)                                                                                                    tGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACgg                             <  2:112070/139‑1 (MQ=255)                                                                                                      cAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAAc                                                                         <  1:915502/93‑1 (MQ=255)                                                                                                      cAGGCCACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAAc                                                                         >  2:915502/1‑93 (MQ=255)                                                                                                           cACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACgg                                                                                   >  1:189713/1‑78 (MQ=255)                                                                                                           cACTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACgg                                                                                   <  2:189713/78‑1 (MQ=255)                                                                                                            aCTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCg                     >  2:97735/1‑139 (MQ=255)                                                                                                            aCTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGC‑CCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCgg                    >  1:1234336/1‑139 (MQ=255)                                                                                                             cTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    >  1:252949/1‑123 (MQ=255)                                                                                                             cTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    <  2:252949/123‑1 (MQ=255)                                                                                                             cTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    >  2:1539810/1‑123 (MQ=255)                                                                                                             cTGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACg                                    <  1:1539810/123‑1 (MQ=255)                                                                                                              tGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAAc                                                                         >  1:1296551/1‑85 (MQ=255)                                                                                                              tGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAAc                                                                         <  2:1296551/85‑1 (MQ=255)                                                                                                              tGCGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGC‑CCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCa                  <  2:1234336/139‑1 (MQ=255)                                                                                                                cGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCgg                                                                                                   >  1:517203/1‑57 (MQ=255)                                                                                                                cGTGGATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCgg                                                                                                   <  2:517203/57‑1 (MQ=255)                                                                                                                   ggATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCAc                 <  1:1309614/136‑1 (MQ=255)                                                                                                                   ggATGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCAc                 >  2:1309614/1‑136 (MQ=255)                                                                                                                     aTGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGAtt                                                                                                <  1:332289/55‑1 (MQ=255)                                                                                                                     aTGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGAtt                                                                                                >  2:332289/1‑55 (MQ=255)                                                                                                                     aTGAAGAAAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGAACGGCACGCTCa            <  1:934590/139‑1 (MQ=255)                                                                                                                        aagaaAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCt                                                                              >  2:304696/1‑70 (MQ=255)                                                                                                                        aagaaAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCt                                                                              <  1:304696/70‑1 (MQ=255)                                                                                                                          gaaAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGgcgc                                                           <  2:448034/87‑1 (MQ=255)                                                                                                                          gaaAGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGgcgc                                                           >  1:448034/1‑87 (MQ=255)                                                                                                                             aGGCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCACGCTCAACGGCATg    >  2:1464957/1‑139 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGTCTTa                                                                                                        >  1:956457/1‑37 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTa                                                                                                        <  2:956457/37‑1 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACgg                                                                                   >  2:32650/1‑58 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACgg                                                                                   <  1:32650/58‑1 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACgg                             <  2:848367/112‑1 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACgg                             >  1:848367/1‑112 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCAc                 <  2:648903/124‑1 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCAc                 >  1:648903/1‑124 (MQ=255)                                                                                                                               gCGACATGATGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCACGCTCAACGGCATGAt  >  1:693210/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                   AACGCATCTTTGTTGCCTCGATTCGTGGTGATGTGCAGCAGCAGGTAAAAACGCTGGAGAAAACCACCGTTGCCGACACCAATACCGAGTGGAGCAAGTTGCAGGCCACTGCGTGGATGAAGAAAGGCGACATGGTGAACGATATCAAACCGATCTGGGCTTATGCGGATTCGTTGTACAACGGCACCTGTAACCAGTGCCACGGCGCACCGGAAATCGCCCACTTTGACGCTAACGGTTGGATCGGCACGCTCAACGGCATGAT  >  NC_000913/1058880‑1059144

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