Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	470,579	A→T	14.3%	Y570F (TAT→TTT)	mdlA →	ABC transporter family protein MdlA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	470,579	0	A	T	14.3%	73.9 / 5.2	28	Y570F (TAT→TTT)	mdlA	ABC transporter family protein MdlA
Reads supporting (aligned to +/- strand):  ref base A (12/12);  new base T (2/2);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.98e-01

GGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCTATTTTATCGACAATATGGTAGCGAAAAATAACCTGCCGTTGAAA‑GTGGTTGC  >  NC_000913/470376‑470829                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            ggTGGACGGACTCACTGACCTCCAGTTCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGATCGGTAATCATCAGTGCCCATCGCCTTGCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAg                                                                                                                                                                                       <  1:312671/274‑1 (MQ=255)           cgCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTgg                                                                                                                                                                             <  1:85302/274‑1 (MQ=255)                   gCACCAGATCCTGCATAACCTGCGTCAGTGCGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGTCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGActc                                                                                                                                                                      >  1:237222/1‑273 (MQ=255)                     aCCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCa                                                                                                                                                                   >  1:285904/1‑274 (MQ=255)                               gCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGa                                                                                                                                                                                                    >  1:288447/1‑231 (MQ=255)                               gCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGAAGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGt                                                                                                                                                          >  1:127879/1‑273 (MQ=255)                                  tAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGc                                                                                                                                                                                                 <  2:288447/231‑1 (MQ=255)                                               tGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGcc                                                                                                                                          <  1:204301/273‑1 (MQ=255)                                                                           gTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGgttt                                                                                                                                                                                         >  2:195366/1‑198 (MQ=255)                                                                              cccATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGc                                                                                                                                                                                      <  1:195366/198‑1 (MQ=255)                                                                                                                           tGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCAGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCg                                                                                                                                                                          >  1:358335/1‑165 (MQ=255)                                                                                                                           tGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCAGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCg                                                                                                                                                                          <  2:358335/165‑1 (MQ=255)                                                                                                                                          tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTTTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGAtt                                                                                                                    >  1:125806/1‑204 (MQ=255)                                                                                                                                             cccAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTTTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCg                                                                                                                 <  2:125806/204‑1 (MQ=255)                                                                                                                                                aGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGa                                                                                                                                                                                                    >  1:266592/1‑118 (MQ=255)                                                                                                                                                  cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACtgtg                                                                                                                                                                              >  2:147566/1‑138 (MQ=255)                                                                                                                                                   gTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGc                                                                                                                                                                                                 <  2:266592/118‑1 (MQ=255)                                                                                                                                                     ggCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGcc                                                                                                                                                                           <  1:147566/138‑1 (MQ=255)                                                                                                                                                       ggaTCATGGGGTGCTGGAACAACAAGGCGGCTGGTATCGCGATCTGTATCGCGAGCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCTATTTTATCGACAATATGGTAg                                  <  2:83331/271‑1 (MQ=255)                                                                                                                                                             tgatgTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTg                                                                                                                                                                                                                                                   >  2:211671/1‑58 (MQ=255)                                                                                                                                                                tgtgCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAg                                                                                                                                                                                                                                                <  1:211671/58‑1 (MQ=255)                                                                                                                                                                     tGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTAAGTGGGCCGCTGCTTATCAGCTATTTTATCGACAATATgg                                     >  1:83331/1‑256 (MQ=255)                                                                                                                                                                         acaacaAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCggcggcggc                                                                                    >  2:304335/1‑205 (MQ=255)                                                                                                                                                                            acaaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGa                                                                                 <  1:304335/205‑1 (MQ=255)                                                                                                                                                                                agggggtGGTTTCGCGATATTTATCGCTTTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCTATTTTATCGACAATATGGTAGCGAAAAATAACCTGCCGTTGAAA‑GTg       <  1:128465/268‑1 (MQ=255)                                                                                                                                                                                      tGGTATCGCGATATGTATCGCTTTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCCGCAGTGCGTGGGCCGCTGCTTATACGCCATTTTTACGACAATATGGGAGCGGAAAAAACCCGACCGGTGAAACGTGGtgct  >  2:128465/1‑271 (MQ=255)                                                                                                                                                                                                    gTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAgcgc                                                                                                                                                              >  2:150961/1‑104 (MQ=255)                                                                                                                                                                                                       tCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                                                                           <  1:150961/104‑1 (MQ=255)                                                                                                                                                                                                            |                                                                                                                                                                                                                                                            GGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAACGGTAATCATCAGTGCCCATCGCCTTTCTGCACTGACGGAAGCCAGTGAAATTATTGTGATGCAGCACGGACATATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCTATTTTATCGACAATATGGTAGCGAAAAATAACCTGCCGTTGAAA‑GTGGTTGC  >  NC_000913/470376‑470829

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