Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,059,533	A→G	H93R (CAC→CGC)	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,533	0	A	G	100.0%	135.5 / NA	38	H93R (CAC→CGC)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/16);  total (22/16)

ATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTT  >  NC_000913/1059397‑1059668                                                                                                                                         |                                                                                                                                        aTCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGTTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGt                                                                                                                                      >  1:377803/1‑140 (MQ=255) aTCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATcgcgg                                                                                                                                      >  2:865308/1‑139 (MQ=255)         agagGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAAcg                                                                                                                              >  2:146714/1‑140 (MQ=255)           agGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAAcgcg                                                                                                                            >  2:693280/1‑140 (MQ=255)             ggCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGCTAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGCACACAACgcggc                                                                                                                          >  2:408252/1‑140 (MQ=255)             ggCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACgcggc                                                                                                                          >  1:505230/1‑140 (MQ=255)                aTTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGcgcg                                                                                                                       <  1:146714/140‑1 (MQ=255)                            tCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATcc                                                                                                           >  2:674140/1‑140 (MQ=255)                            tCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATcc                                                                                                           >  2:155078/1‑140 (MQ=255)                             cGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCg                                                                                                          >  1:610634/1‑140 (MQ=255)                              gCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATTATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGa                                                                                                         >  2:1002818/1‑140 (MQ=255)                              gCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGa                                                                                                         <  2:774004/140‑1 (MQ=255)                                 cTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATgg                                                                                                      <  1:84022/140‑1 (MQ=255)                                    gggggCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTAc                                                                                                   >  2:68847/1‑140 (MQ=255)                                         cTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTg                                                                                              >  1:229919/1‑140 (MQ=255)                                               gcgcgACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTgg                                                                                        <  1:161749/140‑1 (MQ=255)                                                    aCGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCg                                                                                   <  2:288818/140‑1 (MQ=255)                                                      ggTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTa                                                                                 >  2:742352/1‑140 (MQ=255)                                                               tGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCAtc                                                                        >  2:171868/1‑140 (MQ=255)                                                                ggTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCAtct                                                                       >  1:351007/1‑140 (MQ=255)                                                                 gTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCAtctc                                                                      >  1:523408/1‑140 (MQ=255)                                                                    gCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGc                                                                   >  2:552883/1‑140 (MQ=255)                                                                     cTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCg                                                                  <  1:68847/140‑1 (MQ=255)                                                                       ttGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGAt                                                                >  2:165463/1‑140 (MQ=255)                                                                                   aaCCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAgcgc                                                    <  1:778379/140‑1 (MQ=255)                                                                                      cGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGt                                                 <  2:377803/140‑1 (MQ=255)                                                                                      cGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGt                                                 <  1:674140/140‑1 (MQ=255)                                                                                        ttCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGa                                               <  2:610634/140‑1 (MQ=255)                                                                                          cGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATa                                             <  1:693280/140‑1 (MQ=255)                                                                                                ggATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGtt                                       >  1:133019/1‑140 (MQ=255)                                                                                                        atCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTgcgc                               <  2:523408/140‑1 (MQ=255)                                                                                                         tCCGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGcgcg                              <  1:171868/140‑1 (MQ=255)                                                                                                           cGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTg                            >  2:121943/1‑140 (MQ=255)                                                                                                           cGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTg                            >  1:544006/1‑140 (MQ=255)                                                                                                           cGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTg                            >  2:600175/1‑140 (MQ=255)                                                                                                           cGTCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGAGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTg                            >  2:880114/1‑140 (MQ=255)                                                                                                             tCGAAAATGATTGCCGGATTGCCGGATCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAg                          >  2:737628/1‑140 (MQ=255)                                                                                                                       ttgccggattgccggaTCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGaa                >  1:89728/1‑140 (MQ=255)                                                                                                                           cggattgccggatCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATGCCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGccc            <  2:65953/140‑1 (MQ=255)                                                                                                                                tgccggatCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGAc       <  1:633124/140‑1 (MQ=255)                                                                                                                                  ccggatCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACAt     <  2:545654/140‑1 (MQ=255)                                                                                                                                     gatCGCGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGtt  <  1:188064/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        ATCTCGAAAGAGGGCATTCTTACCGGGTCGCACTGGGGGGCTATCCGCGCGACGGTGAAGGATGGTCGCTTTGTGGCGGCAAAACCGTTCGAACTGGATAAATATCCGTCGAAAATGATTGCCGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTT  >  NC_000913/1059397‑1059668

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