Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,712,919	C→T	Q51* (CAA→TAA)	dtpA →	dipeptide/tripeptide:H(+) symporter DtpA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,712,919	0	C	T	86.2%	71.9 / 6.3	29	Q51* (CAA→TAA)	dtpA	dipeptide/tripeptide:H(+) symporter DtpA
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (17/8);  total (19/10)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.50e-01

GTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGTCGGTAACGGCCTGTTTAAAGCTA  >  NC_000913/1712674‑1713162                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    gTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGACTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTg                                                                                                                                                                                                                                                >  2:192798/1‑251 (MQ=255)                ataACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAGAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAg                                                                                                                                                                                                                                   <  2:81632/249‑1 (MQ=255)                 taACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCgg                                                                                                                                                                                                                                >  1:85445/1‑251 (MQ=255)                      aTTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAAATGGGTATGTCTGGAGCGGATTc                                                                                                                                                                                                                            >  2:409079/1‑250 (MQ=255)                      aTTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCa                                                                                                                                                                                                                           <  1:325564/251‑1 (MQ=255)                           gcgTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCAc                                                                                                                                                                                                                      <  2:269353/251‑1 (MQ=255)                                                        aGCACCCCCGTTAATAAGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGttt                                                                                                                                                                                         >  1:341012/1‑251 (MQ=255)                                                                  ttcatatGGGATGTAAAACACGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGATATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTAtggtctggt                                                                                                                                                                               >  2:458859/4‑251 (MQ=255)                                                                       atGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCt                                                                                                                                                                           >  1:433515/1‑250 (MQ=255)                                                                                            aaaaGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGGTAGGTGAc                                                                                                                                                      >  1:415440/1‑250 (MQ=255)                                                                                                                      caaCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCTCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAAtt                                                                                                                           <  1:312341/251‑1 (MQ=255)                                                                                                                       aaCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAAtt                                                                                                                           >  1:407811/1‑250 (MQ=255)                                                                                                                        aCTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGGGACAAGGTACTGGGTACTAAACGCGTAATTAt                                                                                                                         >  1:324042/1‑251 (MQ=255)                                                                                                                           gAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCt                                                                                                                      >  1:401419/1‑251 (MQ=255)                                                                                                                                                    aacaacCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGt                                                                                             <  2:407811/251‑1 (MQ=255)                                                                                                                                                          cGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGct                                                                                       >  1:290491/1‑251 (MQ=255)                                                                                                                                                            aaGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATg                                                                                         >  1:279161/1‑247 (MQ=255)                                                                                                                                                            aaGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATg                                                                                         <  2:279161/247‑1 (MQ=255)                                                                                                                                                                       tctcATCTTCTCGATTGCGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGATATCGGCGGCTGGTTAGGTGACAAGGTACTGGGCACTCAACGCGTAATTATGCTCGGCGCTATTGTGGTGGCGATTGGTTATGCGCTGGTTGCctggt                                                                          >  2:115945/1‑251 (MQ=255)                                                                                                                                                                                  cGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCacgac                                                               <  2:401419/251‑1 (MQ=255)                                                                                                                                                                                    aTTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGc                                                             >  1:490502/1‑251 (MQ=255)                                                                                                                                                                                           tATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATc                                                      >  1:131131/1‑251 (MQ=255)                                                                                                                                                                                              gggAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGtt                                                   >  1:171590/1‑251 (MQ=255)                                                                                                                                                                                                        tGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTTCCTGGTCTGGTCACGACGACGGTATCGTTTATATGGGtc                                         >  2:289794/1‑250 (MQ=255)                                                                                                                                                                                                         ggTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTAt                                        <  2:244973/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGTCGGTAACgg                <  2:433515/251‑1 (MQ=255)                                                                                                                                                                                                                                 aTGGCTGTTTACCTGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGTCGGTAACg                 <  2:171590/250‑1 (MQ=255)                                                                                                                                                                                                                                              tGGTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGGCGGGAACGGCCTGTTTAAAGCt   >  1:104255/1‑251 (MQ=255)                                                                                                                                                                                                                                                gTTAAATAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGTCGGTAACGGCCTGTTTAAAGCTa  >  1:397377/1‑250 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    GTTTCCCCTCCCACTATAACAATTGCGCGTATGTTTTTTATACATAACGCGAGAAAGCACCCCCGTTAATATGGGATGTAAAAAAAGAGGTAAAAGTGTCCACTGCAAACCAAAAACCAACTGAAAGCGTCAGTTTGAACGCTTTCAAACAACCGAAGGCGTTCTATCTCATCTTCTCGATTGAGTTATGGGAACGTTTTGGTTATTACGGCCTACAAGGAATTATGGCTGTTTACCTGGTTAAACAACTGGGTATGTCTGAAGCGGATTCAATCACCCTTTTCTCTTCCTTTAGTGCCCTGGTTTATGGTCTGGTCGCTATCGGCGGCTGGTTAGGTGACAAGGTACTGGGTACTAAACGCGTAATTATGCTCGGCGCTATTGTGCTGGCGATTGGTTATGCTCTGGTTGCCTGGTCTGGTCACGACGCCGGTATCGTTTATATGGGTATGGCGGCTATTGCGGTCGGTAACGGCCTGTTTAAAGCTA  >  NC_000913/1712674‑1713162

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