Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,060,354	A→G	T367A (ACG→GCG)	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,060,354	0	A	G	100.0%	113.7 / NA	32	T367A (ACG→GCG)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/16);  total (16/16)

TGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGC  >  NC_000913/1060217‑1060479                                                                                                                                          |                                                                                                                              tgGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCg                                                                                                                             >  1:764838/1‑140 (MQ=255) tgGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCg                                                                                                                             >  1:451135/1‑140 (MQ=255)      tttGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCg                                                                                                                             >  1:53853/1‑135 (MQ=255)             aGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTgct                                                                                                                 <  1:122455/140‑1 (MQ=255)               ttCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTgctgg                                                                                                               >  2:4626/1‑140 (MQ=255)               ttCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTgctgg                                                                                                               >  1:834741/1‑140 (MQ=255)                         aTCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTgctggc                                                                                                              <  1:596658/131‑1 (MQ=255)                         aTCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTgctggc                                                                                                              >  2:596658/1‑131 (MQ=255)                         aTCTGCTGGGTGAGAAACACGATCACCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGACGCCGAAACCATCCGTGCCCTGGCACGGCCGCGGGCGGCGAACAGAGAGCAGATTTTTGCTGGGTGgggcgtg                                                                                                     >  2:651577/1‑140 (MQ=255)                          tCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGgtgcgtgc                                                                                                    <  1:215294/140‑1 (MQ=255)                             gctgGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGc                                                                                                 <  1:510322/140‑1 (MQ=255)                               tgGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGt                                                                                               <  2:51284/140‑1 (MQ=255)                                  gTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATg                                                                                            >  1:552028/1‑140 (MQ=255)                                    gagaAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgca                                                                                          <  2:451135/140‑1 (MQ=255)                                    gagaAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgca                                                                                          <  2:764838/140‑1 (MQ=255)                                      gaAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATgcagc                                                                                        >  1:841961/1‑140 (MQ=255)                                        aaGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCAc                                                                                      <  2:552028/140‑1 (MQ=255)                                            cGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCTGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTg                                                                                  <  2:834741/140‑1 (MQ=255)                                            cGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTg                                                                                  <  1:4626/140‑1 (MQ=255)                                            cGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTg                                                                                  >  1:570724/1‑140 (MQ=255)                                               tCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAAc                                                                               >  2:549327/1‑140 (MQ=255)                                                                 aTGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATtg                                                             <  2:841961/140‑1 (MQ=255)                                                                     gCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGt                                                         >  1:925667/1‑140 (MQ=255)                                                                             aCTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTggcgg                                                 <  2:925667/140‑1 (MQ=255)                                                                              cTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTggcggc                                                >  2:643483/1‑140 (MQ=255)                                                                                gACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGTGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCAGGGTGAACGGTGGGCGTGGATGATTGTGGTTCTGGGGGCGa                                              >  2:838917/1‑140 (MQ=255)                                                                                                 aaaCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGc                             <  1:835170/140‑1 (MQ=255)                                                                                                   aCCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCt                           >  2:183068/1‑140 (MQ=255)                                                                                                      aTTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGcc                        <  1:549327/140‑1 (MQ=255)                                                                                                       ttCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCa                       >  1:133296/1‑140 (MQ=255)                                                                                                           tGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGc                             >  2:246023/1‑130 (MQ=255)                                                                                                           tGGGCTGGCGCGGCAGATGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGc                             <  1:246023/130‑1 (MQ=255)                                                                                                                            tGGCGGCGAACAGAGCGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTtggc  <  1:838917/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                              TGGGTTTTGAGCAGTTCCTGCCGTATCTGCTGGGTGAGAAAGACGGTCAGCCGAAAGATGCCGCATGGGCTGAAAAACTGACCGGCATTGATGCCGAAACCATTCGTGGGCTGGCGCGGCAGATGGCGGCGAACAGAACGCAAATTATTGCTGGCTGGTGCGTGCAGCGTATGCAGCACGGTGAACAGTGGGCGTGGATGATTGTGGTTCTGGCGGCGATGCTGGGGCAAATTGGCCTGCCAGGTGGTGGTTTTGGTTTTGGC  >  NC_000913/1060217‑1060479

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