Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,956,975	G→T	H11Q (CAC→CAA)	torZ ←	trimethylamine N‑oxide reductase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,956,975	0	G	T	100.0%	64.0 / NA	24	H11Q (CAC→CAA)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (18/6);  total (18/6)

ACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGC  >  NC_000913/1956838‑1957108                                                                                                                                          |                                                                                                                                      aCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCCCCAGCGCTCCAGCGGCAATGCCGCTTTGttt                                                                                                                                  >  1:226738/1‑143 (MQ=255)            gTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCCTTGTTAAATAAATTCAc                                                                                                                       >  1:332831/1‑143 (MQ=255)                 tCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCtt                                                                                                                  <  2:93760/143‑1 (MQ=255)                       cTACGTTCATCGCTCCCCATCGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAAt                                                                                                            <  2:75945/143‑1 (MQ=255)                          cGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAACGTc                                                                                                         >  1:365130/1‑143 (MQ=255)                               aTCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGCTGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGGCCTGTTAATGTCATAAc                                                                                                    >  2:428434/1‑143 (MQ=255)                                  tcTCCCCAGCGTCCCGCGGTGGGGATTTTACCGCCCTTCTCTTCGGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCAACAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtc                                                                                                 <  2:580845/142‑1 (MQ=255)                                       ccAATGTCCCGCGGTGAGGATTTGACCGCCCTTCTCTTCTGCCCATGCGGGGGGCGGTGCTGCCGATGTCACCACCAGCGCCCCAGCGCCCATACCCCTTTGTTTAATAACTTAACGTCTAGTTAATGTCATACCtcaccgcc                                                                                            >  2:254970/1‑136 (MQ=255)                                            gTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATc                                                                                       <  2:47342/143‑1 (MQ=255)                                              cccGGGGTGGGGCATCTACCTCCGTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCCCCAGCGCTCCAGCGTCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCaa                                                                                     <  2:3157/143‑1 (MQ=255)                                               ccGCGGTGAGGATTTTACCGCCCTTCTCTTGTGCCCATGCGGGGAGCGTTGCTGCCGATGTCAGCACCAGCGCTCCAGCGGCGATACCGCGTCGGTGAATCAATTCACGTCTGGCTAATGTCATAACTTCCTCCCTGGTTAAc                                                                                    >  2:613563/1‑143 (MQ=255)                                                           ttttACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGt                                                                        >  1:212268/1‑143 (MQ=255)                                                                gcGCCCTTCACTTCCGCCCCTGCGGGGAGCGGTGCTGCCGATGTCCCCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTATTACATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTtctc                                                                   <  1:227042/142‑1 (MQ=255)                                                                    ccTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCCCCACCAGCGCTCCAGCGGCAATACTGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGCTCACTGTTGCTCGGtc                                                               >  2:220537/1‑142 (MQ=255)                                                                                      cGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGCGGATCACTGTTGCTCGGTACTATCTTTGGGGTCGTAc                                             >  2:228541/1‑143 (MQ=255)                                                                                         ggAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTgg                                          >  1:495173/1‑143 (MQ=255)                                                                                                           tCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTGCTGGAAATACGCCGTTAAAATg                        >  1:363812/1‑143 (MQ=255)                                                                                                            caccacCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGGTGAGGTCATAACTTCCTCCCTGATCAACGCGGATCACTGTTTCTCGGTAAGATCTTTGGCGTTGTGCTGGCAATGCCGCGCTaaacggc                       >  1:272475/1‑139 (MQ=255)                                                                                                                  caGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTCAACTGTCCCGtt                 >  1:210563/1‑143 (MQ=255)                                                                                                                  caGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTCAGAATAGCCAGtt                 >  1:137830/1‑143 (MQ=255)                                                                                                                        tCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACCTCCTCCCCGATCGACGAGTAGCACGGTTTCTCGGTAATATCTTTGGCGTTGCCCTGGAAATACCGCGTTATAATGTGCATGTCGtttt           >  2:74790/1‑143 (MQ=255)                                                                                                                          cAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTCAAAATGTCCAGTTCGTTTTCg         >  2:342697/1‑143 (MQ=255)                                                                                                                               gCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGCTGCACTGGAAATACCGCGTCACAAAGTCCAGGTCTTTTTCGCTccg    >  2:430228/1‑141 (MQ=255)                                                                                                                                 aaTACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGCTTTCTCTCATGc  >  1:13230/1‑143 (MQ=255)                                                                                                                                          |                                                                                                                                      ACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGC  >  NC_000913/1956838‑1957108

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