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%	97.8 / NA	33	H11Q (CAC→CAA)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (19/14);  total (19/14)

AACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGT  >  NC_000913/1956837‑1957112                                                                                                                                           |                                                                                                                                          aaCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGtt                                                                                                                                       >  1:1002131/1‑143 (MQ=255)  aCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTATCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCCCCACCAGCGCTCCAGCGGCCATACCCCTTTGttt                                                                                                                                      >  1:498069/1‑143 (MQ=255)     aTCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTtaat                                                                                                                                   <  1:796543/143‑1 (MQ=255)       cTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGCAGCGGTGCTGCCGATGTCACCACCACCGCCCCACCGGCCGCCCCGCTTTCTTtaataa                                                                                                                                 >  2:945601/1‑143 (MQ=255)           ccGTCCCTCACGTCTACGTTCGTCGCTCCCCTACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTc                                                                                                                             <  2:1011363/143‑1 (MQ=255)                     cTTCTACGTTCATTGCTCGCCCACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGtt                                                                                                                   <  2:239950/143‑1 (MQ=255)                         tACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATg                                                                                                               <  1:309029/143‑1 (MQ=255)                                 tCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCCCTTTGTTTAATAACTTCACGTCTTGTTAAGGTCATAACt                                                                                                       >  1:306681/1‑143 (MQ=255)                                    cTCCCCATCGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtcc                                                                                                    <  1:22993/143‑1 (MQ=255)                                                ccGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAAc                                                                                        >  2:827642/1‑143 (MQ=255)                                                    ggTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAgg                                                                                    <  1:842333/143‑1 (MQ=255)                                                        aGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCa                                                                                >  2:556920/1‑143 (MQ=255)                                                         ggATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCAc                                                                               <  1:154199/143‑1 (MQ=255)                                                                  cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCg                                                                      >  1:505672/1‑143 (MQ=255)                                                                  cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCg                                                                      >  1:66202/1‑143 (MQ=255)                                                                     gcgtctcgtttGCCCGCGCGGGTCGCGGTCCTTCCGGGGTCGCCTCCGGCGCGCCTTCGGCAATACCGCTTTGTCTGATATATGCGTGTCTTGTTAATGTCATCGGTTTCTCCCTGATCAACGTGGATCACTGTTTCTCGGTa                                                                   <  2:601084/140‑1 (MQ=255)                                                                           cttctGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGCGGGTCACTGTTTCCCGCGAATATCt                                                             >  2:537468/1‑143 (MQ=255)                                                                            ttctGCGCTTTCGGCGCGCGGTGCTGCCGCGGTTACCGCCAGCGCTCCAGCGGCAATCCCGCCTTGCTCTCTCGCGTCACTTCTTGTTTGTGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCtt                                                            <  2:108442/143‑1 (MQ=255)                                                                                    aTGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTg                                                    <  1:763877/143‑1 (MQ=255)                                                                                          ggAGCGCTGCTGCCGAAGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTgg                                              <  1:60541/143‑1 (MQ=255)                                                                                                     gcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAAGATCTTTGGCGTTGTACTGGAAATACCGCGt                                   >  1:131972/1‑143 (MQ=255)                                                                                                      ccGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACGGGAAATAACGCGtt                                  >  2:427635/1‑143 (MQ=255)                                                                                                           gTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAAt                             <  1:547069/143‑1 (MQ=255)                                                                                                            tCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACTGCGTTAAAATg                            >  2:506009/1‑143 (MQ=255)                                                                                                            tCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGCAACTACTGCGTTAAAATg                            >  2:434229/1‑143 (MQ=255)                                                                                                              accaccAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTCGAAATACCGCGGTAAAAGGTc                          >  2:755144/1‑143 (MQ=255)                                                                                                                caccaGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTCGCGTTGTACTGGAAATACCGCGTTAAAATGACCa                        >  2:510855/1‑143 (MQ=255)                                                                                                                   caGCGCTCCAGCGGCAATACCGCTTTGTTTAATACATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGtt                     <  1:27437/143‑1 (MQ=255)                                                                                                                      cgcTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGAGCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGt                  >  1:246428/1‑143 (MQ=255)                                                                                                                       gcTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCGACGAGGATCTCTGTTTCTCGGTAATATCTTTGTCGTTGTACTGGAAATACCGGGTTTCAATGTCCAGTTCGtt                 >  1:26415/1‑143 (MQ=255)                                                                                                                                   aTACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGCCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATCCt     >  1:541939/1‑143 (MQ=255)                                                                                                                                   aTACCGCTTAGTTTAATCAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGCCTCACTATTTCTCGGTACTCCCTTTGACGTTGTACTGGAAATACCGCGTCAAACTCGCCAGTCCGTTCTCCCtactact     >  2:993823/1‑137 (MQ=255)                                                                                                                                      ccGCTTTGTTTAACAAGTTCACGGCGTGTCAATGTCATAACTTCCTCCCTGATCCACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGt  <  1:405004/143‑1 (MQ=255)                                                                                                                                           |                                                                                                                                          AACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGATCAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGT  >  NC_000913/1956837‑1957112

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