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%	153.8 / NA	45	W59G (TGG→GGG)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (20/25);  total (20/25)

CGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑GGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATC  >  NC_000913/1059292‑1059562                                                                                                                                            |                                                                                                                                     cGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑g                                                                                                                                     <  2:1105813/140‑1 (MQ=255)      tCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑gggggg                                                                                                                                >  1:1283631/1‑140 (MQ=255)                cTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcg                                                                                                                       >  1:1376186/1‑140 (MQ=255)                cTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcg                                                                                                                       <  1:451309/140‑1 (MQ=255)                  cggcggCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCcgcgcg                                                                                                                     <  1:493073/140‑1 (MQ=255)                   ggcggcTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGa                                                                                                                    <  1:968181/140‑1 (MQ=255)                   ggcggcTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGa                                                                                                                    <  2:1176942/140‑1 (MQ=255)                     cggcTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACg                                                                                                                  >  2:777926/1‑140 (MQ=255)                        cTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTg                                                                                                               >  2:1139725/1‑140 (MQ=255)                                           cTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTgg                                                                                            >  2:376022/1‑140 (MQ=255)                                                   gTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCaaaa                                                                                    >  1:1010350/1‑140 (MQ=255)                                                    tCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAAc                                                                                   >  2:783846/1‑140 (MQ=255)                                                    tCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAAc                                                                                   >  2:1024720/1‑140 (MQ=255)                                                    tCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAAc                                                                                   <  2:1010350/140‑1 (MQ=255)                                                          ttAACGCCGCGACGTGCGACTGCGGCGCAAGGGGGGCCTCACGCTGGCATCTGGAAAGAGGGGATTCTTACCGGGTCGGACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                             <  2:222430/140‑1 (MQ=255)                                                          ttAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                             <  2:125539/140‑1 (MQ=255)                                                          ttAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCg                                                                             >  2:1131108/1‑140 (MQ=255)                                                            aaCGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGaa                                                                           <  1:257186/140‑1 (MQ=255)                                                             aCGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAAc                                                                          <  1:1024720/140‑1 (MQ=255)                                                              cgccgcGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACt                                                                         >  2:317376/1‑140 (MQ=255)                                                              cgccgcGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACt                                                                         >  1:300843/1‑140 (MQ=255)                                                                 cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGa                                                                      >  1:935/1‑140 (MQ=255)                                                                 cgcgACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGa                                                                      >  2:381511/1‑140 (MQ=255)                                                                       gTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAata                                                                <  2:935/140‑1 (MQ=255)                                                                        tGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAatat                                                               >  2:810953/1‑140 (MQ=255)                                                                         gCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATc                                                              <  1:457606/140‑1 (MQ=255)                                                                              tGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCg                                                         <  2:360930/140‑1 (MQ=255)                                                                                cggcgCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGaa                                                       >  2:48528/1‑140 (MQ=255)                                                                                  gcgcAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCCATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGaaaa                                                     >  2:200943/1‑140 (MQ=255)                                                                                     cAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATga                                                  >  2:845272/1‑140 (MQ=255)                                                                                           gcgACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccg                                            <  1:48528/140‑1 (MQ=255)                                                                                            cgACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCCATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATgattgccgg                                           <  1:200943/140‑1 (MQ=255)                                                                                                    gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  1:317376/140‑1 (MQ=255)                                                                                                    gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  2:386250/140‑1 (MQ=255)                                                                                                    gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  1:12079/140‑1 (MQ=255)                                                                                                    gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  1:672274/140‑1 (MQ=255)                                                                                                    gCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGattgccggattgccgg                                   <  2:83488/140‑1 (MQ=255)                                                                                                      tGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATtgccggattgccggat                                 <  1:1267339/140‑1 (MQ=255)                                                                                                                aaaGAGGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTAcaca                       >  1:611783/1‑140 (MQ=255)                                                                                                                   gagGGCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAAcg                    <  1:1139725/140‑1 (MQ=255)                                                                                                                     gggCATTCTTACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAAcgcg                  <  1:539437/140‑1 (MQ=255)                                                                                                                              tACCGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTAtt         <  1:777926/140‑1 (MQ=255)                                                                                                                                ccGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg       >  2:679210/1‑140 (MQ=255)                                                                                                                                ccGGGTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCg       <  1:845272/140‑1 (MQ=255)                                                                                                                                    gTCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTAt   <  2:463513/140‑1 (MQ=255)                                                                                                                                     tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  >  2:181836/1‑140 (MQ=255)                                                                                                                                     tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  >  1:639731/1‑140 (MQ=255)                                                                                                                                     tCGCACTGGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATc  <  1:810953/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                     CGTTTTCTGGCACAACTCGGCGGCTTAACCGTCGCCGGGATGCTGGGGCCGTCATTGTTAACGCCGCGACGTGCGACTGCGGCGCAAGCGGCGACTGACGCTGTCATCTCGAAAGAGGGCATTCTTACCGGGTCGCACT‑GGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATC  >  NC_000913/1059292‑1059562

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: ‑