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

AAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTC  >  NC_000913/1956836‑1957114                                                                                                                                            |                                                                                                                                             aaaCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTttgt                                                                                                                                           >  2:789403/1‑143 (MQ=255)   aCAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGttt                                                                                                                                         <  2:793438/143‑1 (MQ=255)    cAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTGCCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTta                                                                                                                                        <  2:244603/143‑1 (MQ=255)     aaTCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTtaa                                                                                                                                       <  2:581715/143‑1 (MQ=255)        cctGCCGTCCTTCACTTCTACGTTCATCGCTCCCCCACGTCCCGCGGTGAGGATTTTACCGCCCTCCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTtaataa                                                                                                                                    <  2:461337/141‑1 (MQ=255)         tgtccGTCCTTCTCTTCTACGTTCATCGCTCCCACAGGACTCTGTAGAAGGATTTTACCGCCCTTCTCTCCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATaaa                                                                                                                                   <  2:503848/140‑1 (MQ=255)          cgCCGTCCTTCAGTTCTACGTTCATCGCTCCCCTACGTCCCGCGGTGACGATTTTACCGCCCTTCTCTTCGGCACATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGTGCTGCAGCGGCAATACCGCTTTGTTTAATAAAt                                                                                                                                  <  2:246446/142‑1 (MQ=255)                      cTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCTGTGCTGCCGCTGTCACCACCAGCGCTCCAGCGGCAATACCGCTTCGTTTAATAAATTCACGTCTTGtt                                                                                                                      >  2:396877/1‑143 (MQ=255)                                   cGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCCTGTTAATGTCATCACtt                                                                                                         >  2:48673/1‑143 (MQ=255)                                        cccAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtcctcc                                                                                                    >  1:388423/1‑143 (MQ=255)                                        cccAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTtcctcc                                                                                                    >  1:1098970/1‑143 (MQ=255)                                        cccAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATATAATAACTtcctcc                                                                                                    >  2:800291/1‑143 (MQ=255)                                          cAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTGTAATTAATTCCCGTCTTGTTAATGTCATAACTTCCTCCCt                                                                                                  >  1:95461/1‑143 (MQ=255)                                            ccGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGa                                                                                                <  1:413419/142‑1 (MQ=255)                                              ggcccGCGTCGAGGATGTGACCGCCCTTCTCTTCTGCCCTTGCGCGGAGCTGTGCTGCCTATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑c                                                                                             <  2:130118/141‑1 (MQ=255)                                                    cgGTGAGGATTTTACCGCCCTGCTCCTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAg                                                                                       <  2:691433/143‑1 (MQ=255)                                                          ggAGTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCAc                                                                                 <  1:1124575/143‑1 (MQ=255)                                                              tttACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATACCTGCCTCCCTGATACAAC‑AGGATCACTgct                                                                             >  2:393915/1‑141 (MQ=255)                                                              tttACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGtt                                                                             >  1:699274/1‑143 (MQ=255)                                                                 aCCGCCCCTCTCGTCTGCCCATGCGGGGAGCGGTGCCGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTtct                                                                          <  1:426273/143‑1 (MQ=255)                                                                   cGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTCAT‑CAACGAGGATCACTGATtctct                                                                        >  2:856976/1‑142 (MQ=255)                                                                    gCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATCACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCgg                                                                       >  1:1048270/1‑143 (MQ=255)                                                                                       gCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTGTCTCGGTAATATCTGTGGCGTTGTa                                                    >  2:413301/1‑143 (MQ=255)                                                                                                    ctgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGCTACTGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTAAGTCCTATCTTTGGCGTCGGACAGGAAATcccac                                       >  2:679490/1‑141 (MQ=255)                                                                                                     tgcCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGGTTAATAAATTCACGTCTTGTTAAGGTCATAACTTCCTCCCTGAT‑CAACGCGGGTCACTGTTTCTCGGTCCTATATCTCGCCTTGTCCTCGCAATACcgcg                                      >  2:879791/1‑143 (MQ=255)                                                                                                         gATGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTATAATAGATACACGTCTTGTTAATGTAATAACGTCCTCCCTGAT‑CAACGAGGATCACTGTGTCTCGGTAATATATTTGGCGATGTACTGGAAATAccaaatcag                                  >  2:120247/1‑136 (MQ=255)                                                                                                          aTGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaaa                                 <  1:83874/143‑1 (MQ=255)                                                                                                          aTGTCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTaaa                                 >  2:709081/1‑143 (MQ=255)                                                                                                             tCACCACCAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTACTATCTTTGGCGTTGTACTGGAAATACCGCGTTAACATg                              >  2:83302/1‑143 (MQ=255)                                                                                                               accaccAGCGCTCCATCGGCAATACCGCTTTGTTTAATAAATTCTCGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTc                            <  1:591288/143‑1 (MQ=255)                                                                                                               accaccAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATACCTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAGTACCGCGTTCAAATggc                            >  2:821368/1‑141 (MQ=255)                                                                                                               accaccAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGCAAATAAAGCGTTACAATGTc                            >  2:809253/1‑143 (MQ=255)                                                                                                               accaccAGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATAGCTTTGGCGGTGTACTGGAAATCGCGCGTTaaaggtgc                            >  2:485366/1‑138 (MQ=255)                                                                                                                ccaccaGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTcc                           <  1:444030/143‑1 (MQ=255)                                                                                                                     aGCGCTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGGAATACCGCGTTAAAATGTCCAGTTc                      >  1:338518/1‑143 (MQ=255)                                                                                                                         cTCCAGCGGCAATACCGCTTTGTTTAATAAATTCGCGTCGTGTTAATGTCCTAACTTCCTCCCTGAC‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGttt                  <  2:855919/143‑1 (MQ=255)                                                                                                                         cTCCAGCGGCAATACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGttt                  <  1:710389/143‑1 (MQ=255)                                                                                                                                   aaTACCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTACTATCTTTGGCGTGGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCAtac        >  2:695229/1‑141 (MQ=255)                                                                                                                                      aCCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTATTATCTTTGGCGTTGTTCTGGAAATACCGCGTTAAAATGTACAGTTCGTTTTAGCTAAtaatgg     >  1:6955/1‑138 (MQ=255)                                                                                                                                      aCCGCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTgg     <  1:1059339/143‑1 (MQ=255)                                                                                                                                        cGCTTTGTTTAATAAATGCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGtt   <  1:885764/143‑1 (MQ=255)                                                                                                                                         gCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTc  >  2:266779/1‑143 (MQ=255)                                                                                                                                         gCTTTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGATAAAATACAAAGATAGTTTTAGATCATGATGGTTc  >  1:1045783/1‑143 (MQ=255)                                                                                                                                            |                                                                                                                                             AAACAATCTTGCCGTCCTTCACTTCTACGTTCATCGCTCCCCAACGTCCCGCGGTGAGGATTTTACCGCCCTTCTCTTCTGCCCATGCGGGGAGCGGTGCTGCCGATGTCACCACCAGCGCTCCAGCGGCAATACCGCTGTGTTTAATAAATTCACGTCTTGTTAATGTCATAACTTCCTCCCTGAT‑CAACGAGGATCACTGTTTCTCGGTAATATCTTTGGCGTTGTACTGGAAATACCGCGTTAAAATGTCCAGTTCGTTTTCGCTCATGCTGGTTC  >  NC_000913/1956836‑1957114

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