Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,061,935	G→A	S46N (AGC→AAC)	torD →	trimethylamine‑N‑oxide reductase‑specific chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,061,935	0	G	A	83.3%	32.3 / ‑4.8	19	S46N (AGC→AAC)	torD	trimethylamine‑N‑oxide reductase‑specific chaperone
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base A (10/5);  minor base C (2/1);  total (12/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCGAT  >  NC_000913/1061689‑1062181                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        ctggttggaaTTGAGAAGTACAACGGAACAGTGGAGCAGATGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGAAGAGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGAACACGCTGACAGCACAACAGATTGCCTGTGTTTACGCATGGCTAGCGACGTTGTTCTCCCGTGAGCTGGACGGTGAACAGCTAACGCAAATCGCCAGTGCGCAGATGGCAGCAGGGTTTTCGTT‑GCTGAAAAACGa                                                                                                                                                                                                                                                      >  2:236775/9‑250 (MQ=255)        aaaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCAGTGGAGATGGTGGCGCAGTGCGCATATGTTCCCGCGTCGCAGGTGAAATCATGCCCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGTCTAGCGCAGTTGTTCTCCCGTGAGCTGGAAGATGAACAACTGACGCCAATCGCCAGTGCGCAGCTGGCTGAATGGTTTACTTT‑GCTGAAAACCGAACcgccgc                                                                                                                                                                                                                                              >  2:205768/1‑251 (MQ=255)         aaTTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCt                                                                                                                                                                                                                                             >  1:221173/1‑251 (MQ=255)           ttGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCAAAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGAAGCAAACCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAAAGAAccgccgacc                                                                                                                                                                                                                                            >  2:220598/1‑247 (MQ=255)                           aaCAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACg                                                                                                                                                                                                                            >  1:15023/1‑250 (MQ=255)                                    gCAGGTGACGGCGTATCACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGACGCCGGTGCAATCATGACCACGCTGACAGCACATCAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGTGCTGTACGATGCACAACTGACGCAAATCTCCCGTGCGAGGATGGCTGAATGGTTTTCGTT‑GTTTAAAACCGAACCGCCGCTCACTGATGCTGTGAACgcgttggtca                                                                                                                                                                                                                   >  2:99986/1‑241 (MQ=255)                                                    accgccccgggggggaggtgggCGCAGTGAAAATATGCTCCCCCGTCGCAGGTGAAAACATGCCCCCGCTGACAGCACAACAGGTGGCCGTTGTTCACGCCGGCCTCGCGCAGTTTTCCTCCCGGGGGGTGGACAATGAACACCGGACGCAAACCGCCGGTGCGCAGATGGCTGAAGGGTTTGTGTG‑GCTAAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTg                                                                                                                                                                                                  <  2:54519/231‑1 (MQ=255)                                                                       tggCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TTACGAAGccc                                                                                                                                                                              >  1:259910/1‑251 (MQ=255)                                                                                      atGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAAAGAGGTGGAAAACCGTATTGCCACGCTGACAGTTCG‑TGAAGATGCCCGTCTGGaagatgccg                                                                                                                                                               >  1:285262/1‑244 (MQ=255)                                                                                                cgctagggtaaaaTCATGCCCACCCTGACAGCAAAACAGATGTCCTTTTTTGACGCCGGGCTGGCGTATTTGTTCCCCGGTGAGGTGGACGTGGAACACCGGACACAATCCGCCAGGGCGCAGATGGCGGAATGGTGTTTGTT‑GCAAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTg                                                                                                                                                     <  2:45257/241‑1 (MQ=255)                                                                                                   gCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGCACCGCCGCTCACTGCGGCGGGGAGCGAGCTGGAAAACCGACTTGCCAAGCGGACAGGACG‑AGACGAGTCCCGGCTGGACAAGGCAGCGGACTTTTTCgg                                                                                                                                                  >  2:109346/1‑251 (MQ=255)                                                                                                                                acaacaGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCACGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGGTGGAAAACCGAGTTGCCACGCTGACAGTACGGTG‑CGATTGCCGTCTGGAACTGGCCGCGGGCTTTTGCGGCGTGGTTCTGGTGACAGGCaagcaagc                                                                                                                       >  1:13013/1‑249 (MQ=255)                                                                                                                                                                   gtttgttCCCCCGTGAGCTGGACGATGAACAACTGCCCAAAACCGCCAGTGCGCATATGGCTGAATGGTTTGCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAg                                                                                  <  1:264397/249‑1 (MQ=255)                                                                                                                                                                                          gATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGGGGACGAGGTGGAAAAACGGATTGCCACGCTGACAGTACG‑TGACGGTGCCCGTCTGGCACTTGCCGCGGACTTTTGCGGCCTGGTTCTGATGACCGCCGAACAAACGGCGGGGCCTTATGCATCGGGCTACAAAAAGGCCGAGCAAGGGAttgacagttttgtatt                                                           >  2:46516/1‑237 (MQ=255)                                                                                                                                                                                             aaccactggACGAAAATCGCCAGTGCCAAGA‑GGGTGAATGGTTTTTGTTTGCTAAAAACCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGa                                                        <  1:1617/243‑1 (MQ=255)                                                                                                                                                                                                        aaaaccccccAGTGCGCAGATGGTTGATTGTTTTTCTTT‑GCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATgg                                             <  1:220598/243‑1 (MQ=255)                                                                                                                                                                                                                                 gAATGGTTTTCGTT‑GCTGAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAAcc                    >  1:6117/1‑251 (MQ=255)                                                                                                                                                                                                                                     gtttttctgt‑GCTAAAAAACGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGc                 <  1:260847/241‑1 (MQ=255)                                                                                                                                                                                                                                                    aaaaACGATCCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCGat  <  1:205768/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        CTGGTGGAAATTGAGAAGTACAACGGAACAGTGGAGCAGGTGACGGCGTTTAACGGCCCCGTGGAGATGGTGGCGCAGTGCGAATATGTTCCCGCGTCGCAGGTGAAATCATGACCACGCTGACAGCACAACAGATTGCCTGTGTTTACGCCTGGCTAGCGCAGTTGTTCTCCCGTGAGCTGGACGATGAACAACTGACGCAAATCGCCAGTGCGCAGATGGCTGAATGGTTTTCGTT‑GCTGAAAAGCGAACCGCCGCTCACTGCGGCGGTGAACGAGCTGGAAAACCGTATTGCCACGCTGACAGTACG‑TGACGATGCCCGTCTGGAACTGGCCGCGGACTTTTGCGGCCTGTTTCTGATGACCGACAAACAAGCGGCGCTGCCGTATGCATCGGCCTACAAACAGGACGAGCAAGAGATTAAACGCTTGTTAGTTGAGGCAGGGATGGAAACCAGCGGCAATTTCAACGAACCGGCAGATCATCTGGCGAT  >  NC_000913/1061689‑1062181

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