Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,059,654	A→C	E133D (GAA→GAC)	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,654	0	A	C	100.0%	184.3 / NA	47	E133D (GAA→GAC)	torA	trimethylamine N‑oxide reductase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (25/22);  total (25/22)

CGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTG  >  NC_000913/1059519‑1059787                                                                                                                                        |                                                                                                                                      cggattgccggatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGccc                                                                                                                                   <  2:415991/140‑1 (MQ=255) cggattgccggatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGccc                                                                                                                                   <  2:246384/140‑1 (MQ=255)      tgccggatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGAc                                                                                                                              >  2:827780/1‑140 (MQ=255)      tgccggatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGAc                                                                                                                              <  1:827780/140‑1 (MQ=255)         cggatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATg                                                                                                                           <  2:615784/140‑1 (MQ=255)           gatCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGtt                                                                                                                         <  2:602211/140‑1 (MQ=255)             acaggTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCt                                                                                                                       >  2:946158/5‑140 (MQ=255)                cGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATg                                                                                                                    <  1:221107/140‑1 (MQ=255)                  tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATgaa                                                                                                                  <  2:554208/140‑1 (MQ=255)                  tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATgaa                                                                                                                  >  1:612457/1‑140 (MQ=255)                  tACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCATAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATgaa                                                                                                                  <  1:304279/140‑1 (MQ=255)                      caACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATg                                                                                                                           >  1:387642/1‑127 (MQ=255)                      caACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATg                                                                                                                           <  2:387642/127‑1 (MQ=255)                         cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTgg                                                                                                           >  1:652506/1‑140 (MQ=255)                         cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTgg                                                                                                           >  2:527618/1‑140 (MQ=255)                         cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTgg                                                                                                           >  2:809895/1‑140 (MQ=255)                         cgcgGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTgg                                                                                                           >  2:179231/1‑140 (MQ=255)                           cggcgCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGaa                                                                                                         <  1:885748/140‑1 (MQ=255)                              cgcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAAcgc                                                                                                      >  1:268705/1‑140 (MQ=255)                              cgcgTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAAcgc                                                                                                      >  2:927237/1‑140 (MQ=255)                                          aTCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACt                                                                                          >  2:302929/1‑140 (MQ=255)                                                  gTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGcc                                                                                  >  1:54428/1‑140 (MQ=255)                                                     cgcgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAg                                                                               >  1:754184/1‑140 (MQ=255)                                                       cgTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTg                                                                             >  1:729987/1‑140 (MQ=255)                                                               ggCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTg                                                                     <  1:946158/140‑1 (MQ=255)                                                                gCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGa                                                                    >  2:872841/1‑140 (MQ=255)                                                                 cTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGAc                                                                   <  1:179231/140‑1 (MQ=255)                                                                      tAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCa                                                              >  2:833762/1‑140 (MQ=255)                                                                                     cGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGa                                               <  2:754184/140‑1 (MQ=255)                                                                                        tACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgg                                            >  1:917639/1‑140 (MQ=255)                                                                                          ccTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                          <  2:268705/140‑1 (MQ=255)                                                                                          ccTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACgggg                                          <  2:729987/140‑1 (MQ=255)                                                                                            tCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGAt                                        <  2:629633/140‑1 (MQ=255)                                                                                                  cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATg                                       <  1:638796/135‑1 (MQ=255)                                                                                                  cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATg                                       >  2:638796/1‑135 (MQ=255)                                                                                                  cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCa                                  >  2:952163/1‑140 (MQ=255)                                                                                                  cgcgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCa                                  >  2:432253/1‑140 (MQ=255)                                                                                                    cgGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATa                                >  2:866940/1‑140 (MQ=255)                                                                                                                 tttGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGAt                   >  1:792128/1‑140 (MQ=255)                                                                                                                  ttGTGCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATg                  <  2:448310/140‑1 (MQ=255)                                                                                                                   tgtgCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGc                 >  1:668210/1‑140 (MQ=255)                                                                                                                     tgCGCGTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTg               >  1:981629/1‑140 (MQ=255)                                                                                                                        gcgTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCg            <  1:927237/140‑1 (MQ=255)                                                                                                                        gcgTGAGCTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCg            >  2:905597/1‑140 (MQ=255)                                                                                                                               cTGGGATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTa     <  1:809895/140‑1 (MQ=255)                                                                                                                                  ggATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg  <  1:432253/140‑1 (MQ=255)                                                                                                                                  ggATGACGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTg  <  1:952163/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                      CGGATTGCCGGATCACGTACACAACGCGGCGCGTATTCGTTATCCGATGGTACGCGTGGACTGGCTGCGTAAGCGCCATCTCAGCGATACCTCCCAGCGCGGTGATAACCGTTTTGTGCGCGTGAGCTGGGATGAAGCCCTCGACATGTTCTATGAAGAACTGGAACGCGTGCAGAAAACTCACGGGCCGAGTGCCTTGCTGACCGCCAGTGGTTGGCAATCGACGGGGATGTTCCATAACGCTTCGGGGATGCTGGCGAAAGCTATTG  >  NC_000913/1059519‑1059787

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