Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,954,604	A→G	F802L (TTT→CTT)	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,954,604	0	A	G	100.0%	103.0 / NA	31	F802L (TTT→CTT)	torZ	trimethylamine N‑oxide reductase 2
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/13);  total (18/13)

TTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTT  >  NC_000913/1954476‑1954729                                                                                                                                 |                                                                                                                              ttCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCg                                                                                                                             >  1:308649/1‑131 (MQ=255)           tGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCg                                                                                                                 >  2:604035/1‑133 (MQ=255)                  aCTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGTCTGATCAAGCGCCGTTAACTTTGGTGCAtt                                                                                                          >  1:283289/1‑133 (MQ=255)                     aTGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGGACAGCTGGATGAACAAGCGCCGTTATCTTCGGTGCATTGcc                                                                                                       >  2:230198/1‑133 (MQ=39)                            ccTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTatat                                                                                                >  2:475194/1‑133 (MQ=255)                              tGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATAttt                                                                                              >  1:46825/1‑133 (MQ=255)                                aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATAttttt                                                                                            <  1:142616/133‑1 (MQ=255)                                aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATAttttt                                                                                            >  1:368488/1‑133 (MQ=255)                                aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATAttttt                                                                                            >  1:18967/1‑133 (MQ=255)                                aTAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTAGGTGCATTGCCCGAATAttttt                                                                                            >  1:214536/1‑133 (MQ=255)                                       tAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGAtata                                                                                     <  2:332310/133‑1 (MQ=255)                                       tAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGAtata                                                                                     <  2:317965/133‑1 (MQ=255)                                          cgcATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATacac                                                                                  >  2:569260/1‑133 (MQ=255)                                              tCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAAc                                                                              <  1:24271/133‑1 (MQ=255)                                                                gCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTAc                                                                    <  1:133088/125‑1 (MQ=255)                                                                gCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTAc                                                                    >  2:133088/1‑125 (MQ=255)                                                                 cATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAgg                                                           <  2:130018/133‑1 (MQ=255)                                                                   tAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCa                                                         <  1:475194/133‑1 (MQ=255)                                                                           gcgcCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGcc                                                   >  1:474488/1‑131 (MQ=255)                                                                           gcgcCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGcc                                                   <  2:474488/131‑1 (MQ=255)                                                                              cccGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTg                                              >  1:603332/1‑133 (MQ=255)                                                                              cccGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTg                                              <  1:569260/133‑1 (MQ=255)                                                                                   ttCCGGGCGGTATTCTCGTCTTACTCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGc                                         <  2:283289/133‑1 (MQ=255)                                                                                               ttCACGTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATccg                             <  2:603332/133‑1 (MQ=255)                                                                                                    gTCTTACGCCTGAACAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGtt                        >  2:20503/1‑133 (MQ=255)                                                                                                                aaCAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCa            >  2:455829/1‑133 (MQ=255)                                                                                                                aaCAGCTGGCTGATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCa            >  1:36741/1‑133 (MQ=255)                                                                                                                    gctggctgATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGc        <  1:20503/133‑1 (MQ=255)                                                                                                                     ctggctgATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTg         <  2:462017/131‑1 (MQ=255)                                                                                                                     ctggctgATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTg         >  1:462017/1‑131 (MQ=255)                                                                                                                          tgATCAAGCGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGtttt  >  1:44355/1‑133 (MQ=255)                                                                                                                                 |                                                                                                                              TTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGCAGCTTTACGTTTGCATAACCTCAGCGCCCGTTTCCGGGCGCTATTCACGTCTTACGCCTGAACAGCTGGCTGATCAAACGCCGTTAACTTCGGTGCATTGCCCGTATATTTTTCGATATACACCAACGCTGAGTTACCGGCACAGGCATTTGCCAGCTGCGAGCTGGGGATATCCGCCGTTAACACGTTCGCACTGCCGTTTT  >  NC_000913/1954476‑1954729

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑