Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,059,436	(G)6→7	coding (181/2547 nt)	torA →	trimethylamine N‑oxide reductase 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,059,430	1	.	G	100.0%	138.6 / NA	39	W59G (TGG→GGG)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (22/17);  total (22/17)

GCGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑GGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATC  >  NC_000913/1059291‑1059562                                                                                                                                             |                                                                                                                                     gCGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACt                                                                                                                                       <  1:898853/140‑1 (MQ=255)  cGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑g                                                                                                                                     <  1:437917/140‑1 (MQ=255)  cGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑g                                                                                                                                     <  2:94271/140‑1 (MQ=255)  cGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑g                                                                                                                                     <  1:779781/140‑1 (MQ=255)    ttttCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑ggg                                                                                                                                   <  2:877420/140‑1 (MQ=255)               aaCTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCG‑‑GCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCGCTGGGGGGGCTATCcgcg                                                                                                                       >  1:493998/1‑140 (MQ=255)                 cTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcg                                                                                                                       >  2:621715/1‑140 (MQ=255)                 cTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcg                                                                                                                       >  2:308719/1‑140 (MQ=255)                 cTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcg                                                                                                                       >  1:912846/1‑140 (MQ=255)                   cggcggCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcgcg                                                                                                                     <  2:881144/140‑1 (MQ=255)                   cggcggCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcgcg                                                                                                                     <  1:689277/140‑1 (MQ=255)                     gcggcTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGAc                                                                                                                   <  2:998734/140‑1 (MQ=255)                     gcggcTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGAAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGAc                                                                                                                   >  2:962851/1‑140 (MQ=255)                                          tGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTtgt                                                                                              <  1:621715/140‑1 (MQ=255)                                                     tCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAAc                                                                                   >  2:289/1‑140 (MQ=255)                                                      cATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCCAAAcc                                                                                  >  1:518289/1‑140 (MQ=255)                                                        ttGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGt                                                                                <  1:850620/140‑1 (MQ=255)                                                           ttAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                             <  2:1077236/140‑1 (MQ=255)                                                             aaCGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                           <  2:641655/140‑1 (MQ=255)                                                             aaCGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                           <  2:912846/140‑1 (MQ=255)                                                             aaCGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                           <  1:871420/140‑1 (MQ=255)                                                               cgccgcGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACt                                                                         >  2:806061/1‑140 (MQ=255)                                                               cgccgcGACGTGCGACTGCGGCACAAGCCGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACt                                                                         >  2:902497/1‑140 (MQ=255)                                                                  cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGa                                                                      >  2:650393/1‑140 (MQ=255)                                                                    cgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCCAAACCGTTCGaa                                                                           >  1:447348/1‑133 (MQ=255)                                                                    cgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                           <  2:447348/133‑1 (MQ=255)                                                                    cgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATa                                                                    >  2:1038715/1‑140 (MQ=255)                                                                                gcggcgCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGa                                                        <  2:77115/140‑1 (MQ=255)                                                                                   gcgcAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGaaaa                                                     <  1:1038715/140‑1 (MQ=255)                                                                                                     gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  2:518289/140‑1 (MQ=255)                                                                                                     gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  1:650393/140‑1 (MQ=255)                                                                                                         tCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATct                               >  2:871440/1‑139 (MQ=255)                                                                                                               cGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTaca                         >  2:164138/1‑140 (MQ=255)                                                                                                                 aaaGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTAcaca                       >  2:217988/1‑140 (MQ=255)                                                                                                                       ggCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACgcgg                 >  2:549964/1‑140 (MQ=255)                                                                                                                        gCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACgcggc                >  2:219192/1‑140 (MQ=255)                                                                                                                         cATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACgcggcg               <  1:308719/140‑1 (MQ=255)                                                                                                                           ttCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGcgcg             <  1:549964/140‑1 (MQ=255)                                                                                                                              ttACCGGATCGCACTGGGGGGGCAATACGTCCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTAt          <  2:764513/140‑1 (MQ=255)                                                                                                                                  cGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGt      >  2:610719/1‑140 (MQ=255)                                                                                                                                     gTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGgcgc              <  2:461089/129‑1 (MQ=255)                                                                                                                                     gTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGgcgc              >  1:461089/1‑129 (MQ=255)                                                                                                                                     gTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACCTACACAACGCGGCCCGTATTCGTTAt   >  1:764513/1‑140 (MQ=255)                                                                                                                                      tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  >  2:1085164/1‑140 (MQ=255)                                                                                                                                      tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  >  1:529806/1‑140 (MQ=255)                                                                                                                                      tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  >  1:120873/1‑140 (MQ=255)                                                                                                                                             |                                                                                                                                     GCGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑GGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATC  >  NC_000913/1059291‑1059562

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