Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,955,440	T→C	Q523R (CAG→CGG)	torZ ←	trimethylamine N‑oxide reductase system III, catalytic subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,955,440	0	T	C	100.0%	127.0 / NA	37	Q523R (CAG→CGG)	torZ	trimethylamine N‑oxide reductase system III, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/15);  total (22/15)

AAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATG  >  NC_000913/1955292‑1955575                                                                                                                                                     |                                                                                                                                        aaaTTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCgg                                                                                                                                       >  2:124614/1‑151 (MQ=255)           cACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAAt                                                                                                                             <  1:242018/151‑1 (MQ=255)                cATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCg                                                                                                                        >  1:354233/1‑151 (MQ=255)                   ttCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTc                                                                                                                     >  1:462042/1‑151 (MQ=255)                   ttCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTc                                                                                                                     >  2:223898/1‑151 (MQ=255)                    tCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTc                                                                                                                     >  1:205223/1‑150 (MQ=255)                     cATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCa                                                                                                                    <  2:593060/150‑1 (MQ=255)                        cTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGt                                                                                                                >  1:2818/1‑151 (MQ=255)                          ttACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTc                                                                                                               >  2:235251/1‑150 (MQ=255)                                       aGATCTCTTTGCCGCCAGGTTGGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTc                                                                                                 <  2:292919/151‑1 (MQ=255)                                               tttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTGAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGAtg                                                                                         >  1:337803/1‑151 (MQ=255)                                               tttCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGAt                                                                                          <  1:377258/150‑1 (MQ=255)                                                  ccGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATgtgg                                                                                      >  2:318618/1‑151 (MQ=255)                                                    gCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATgtggtg                                                                                    <  1:468579/151‑1 (MQ=255)                                                             gAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAAt                                                                           <  2:66222/151‑1 (MQ=255)                                                             gAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAAt                                                                           >  2:382825/1‑151 (MQ=255)                                                                     ccGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGatat                                                                    <  2:265871/150‑1 (MQ=255)                                                                     ccGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGatat                                                                    <  2:2818/150‑1 (MQ=255)                                                                      cGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCt                                                                  <  2:573444/151‑1 (MQ=255)                                                                                          cAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGt                                              >  2:456123/1‑151 (MQ=255)                                                                                                   tACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGc                                      <  2:317677/150‑1 (MQ=255)                                                                                                     cgcgCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCAtt                                   >  1:547429/1‑151 (MQ=255)                                                                                                      gcgcTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTc                                  >  1:71572/1‑151 (MQ=255)                                                                                                       cgcTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCa                                 <  2:133682/151‑1 (MQ=255)                                                                                                          ttCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGa                               >  2:125926/1‑150 (MQ=255)                                                                                                          ttCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGCGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGa                               >  2:275543/1‑150 (MQ=255)                                                                                                            cAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGaaa                             <  2:1475/150‑1 (MQ=255)                                                                                                            cAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAac                            >  1:202808/1‑151 (MQ=255)                                                                                                             aaaTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAacg                           >  2:392279/1‑151 (MQ=255)                                                                                                              aaTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAacga                          <  2:337796/151‑1 (MQ=255)                                                                                                              aaTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAacga                          <  2:314960/151‑1 (MQ=255)                                                                                                               aTTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAacgac                         >  1:356117/1‑151 (MQ=255)                                                                                                                   cGGAGCAACAGCCTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGat                      <  2:340151/150‑1 (MQ=255)                                                                                                                              ccTGCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGttt          >  1:355367/1‑151 (MQ=255)                                                                                                                                 gCTTCATCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGATGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTg       >  1:266429/1‑151 (MQ=255)                                                                                                                                      aTCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATg  >  1:238969/1‑151 (MQ=255)                                                                                                                                      aTCGGCACAATATGCCGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATg  >  2:166941/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                        AAATTTCAGCCACGCCATTTCATCTTTACCTTCGGTATAGATCTCTTTTCCGCCAGGTTTGAGTAATTCCGCAAGATCGGCAAACACGTCAAAATCGTTACGCGCTTCAAATTGCGGAGCAACAGCCTGCTTCATCGGCACAATATGCTGGTTGCTGTAATCACCGGTCATCGTCAAGTCATTGCGCTCAAACGATGTGGTGATCGGTAATACGATATCTGCGTGTTTGGCTGCTGCGGTCCAGTAGCATTCAGAAACGACGATCATCTCCGGTTTCTGCCATG  >  NC_000913/1955292‑1955575

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