Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,517,508	G→T	V40V (GTG→GTT)	ortT →	orphan toxin OrtT

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,517,508	0	G	T	100.0%	70.2 / NA	21	V40V (GTG→GTT)	ortT	orphan toxin OrtT
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (11/10);  total (11/10)

TTTTCTCAAACTGGTGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCT  >  NC_000913/1517301‑1517736                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     ttttCTCAAACTGGTGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGtttt                                                                                                                                                                                           >  2:91574/1‑251 (MQ=255) tgttCTCAAACTGGTGAAAAAACTTGATGCACGGCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGTACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GAGAAGAAACGCATCCGTTTCTTAAGCGCCTTTCGGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGttt                                                                                                                                                                                            <  2:148263/248‑1 (MQ=255)  tttatcAAACTGGGGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGAGTTTTATTCATCGGGCTGGACGCCCGCACATTTCTCTCTGTCAACCCATGCTTGTCTTTTATGCGGTTATCTCAGTAATAGCAACTCCTAACACCGCGTTTCTTTCT‑ACAGTAAATAAACGCATCTGTTTCTTAAGCACCGTTCTGGTTGGCGTCATATGGCCGATGAGTTTCTCGGTCGCGCTGTTCTTTTc                                                                                                                                                                                          >  2:7323/5‑251 (MQ=255)          aCTGGTGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTg                                                                                                                                                                                                        >  1:203147/1‑229 (MQ=255)          aCTGGTGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTg                                                                                                                                                                                                        <  2:203147/229‑1 (MQ=255)                             gCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCCCCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCtt                                                                                                                                                               >  2:54294/1‑251 (MQ=255)                             gCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCtt                                                                                                                                                               >  2:83062/1‑251 (MQ=255)                                                            gTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTgg                                                                                                                                 >  1:71313/1‑250 (MQ=255)                                                             ttttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                                                                                                                                           >  1:28039/1‑223 (MQ=255)                                                             ttttATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACg                                                                                                                                                           <  2:28039/223‑1 (MQ=255)                                                                                 gCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAg                                                                                                           <  1:63685/251‑1 (MQ=255)                                                                                           gTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAgg                                                                                                  <  1:84129/250‑1 (MQ=255)                                                                                                      aaCACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTc                                                                                      >  2:3039/1‑251 (MQ=255)                                                                                                                    ttttATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGaa                                                                         <  1:83062/250‑1 (MQ=255)                                                                                                                           cGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATcgc                                                                 >  2:6028/1‑251 (MQ=255)                                                                                                                            ggTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATcgct                                                                <  2:72531/251‑2 (MQ=255)                                                                                                                                  ggCAGCAATAGCATTTCTTATCACCTGGATTCTTTCTTCC‑GATAAGACACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCAGACGAGTTTCCCGATGGCGCTGTTGTTTTCACTGTTTTAAATCTTAATCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTCAGGAGAAAGTATGTCGCATATGGATGAAGTCATCGCGCGCGTg                                                          >  2:25358/1‑251 (MQ=255)                                                                                                                                           cGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCcgcg                                                 <  2:27021/251‑1 (MQ=255)                                                                                                                                                                   ttctCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCa                          <  1:6028/250‑1 (MQ=255)                                                                                                                                                                              aaGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTAACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAat               >  1:28815/1‑251 (MQ=255)                                                                                                                                                                                            tttCTTAAGCGCCTTTCTGGTTGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTACTCGCCCt  <  2:205786/250‑1 (MQ=255)                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     TTTTCTCAAACTGGTGAAAAAACTTGATGCACGTCAAAAAATGACGCATATTTGCGCGCGTTTTATTCATCTGGCTGGACGCCCGTACATGTCTCTCTATCAACACATGCTTGTTTTTTATGCGGTTATGGCAGCAATCGCATTTCTTATCACCTGGTTTCTTTCTCAC‑GATAAGAAACGCATCCGTTTCTTAAGCGCCTTTCTGGTGGGGGCAACATGGCCGATGAGTTTCCCGGTGGCGCTGTTGTTTTCACTGTTTTAAATCTTACTCACTGGCTTGACGGCACCACGCAGACTTATATCATTTGGATGAATCGATAAATTTCACAAGTGGCTAAGGAGAAAGTATGTCGCATCTGGATGAAGTCATCGCGCGCGTGGATGCCGCGATTGAAGAGAGCGTCATTGCCCATATGAACGAATTATTAATCGCCCT  >  NC_000913/1517301‑1517736

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