Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,309	A→G	L49L (CTA→CTG)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,309	0	A	G	91.7%	60.6 / ‑3.7	24	L49L (CTA→CTG)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (11/11);  total (12/12)
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.86e-01

AAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGA  >  NC_000913/12073‑12545                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               aaaTGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGc                                                                                                                                                                                                                                 <  1:124460/251‑1 (MQ=255)     gAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGAATCACAACCACATGATGAACGAATATAAAGTGGAGACGTTTAGATGGGTAAAATACTTGGTATCGACCTGGGTACTCCCAACTCTTGTGAAGAGATTAAGGGTGGAACCAATCCTCGTGTGCTGGAGCACGACGCAGGCGATCGCACCAAGCCTTCTATCATTGCATATACCAAAGATAGTGAAACTCTGG‑TTGGTCatctggctca                                                                                                                                                                                                                              >  1:357684/1‑241 (MQ=255)          tGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTc                                                                                                                                                                                                                         <  1:179346/250‑1 (MQ=255)          tGGGCAGTTGAAACCAGACGTTTCGCCCCTAGTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGATAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTc                                                                                                                                                                                                                         <  1:226531/250‑1 (MQ=255)          tGGGCAGTTGAAACCAGACGTTTCACCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTc                                                                                                                                                                                                                         <  1:317912/250‑1 (MQ=255)            ggCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAgg                                                                                                                                                                                                                      <  1:215022/251‑1 (MQ=255)                        cAGACGTTTCGCCCCTATTACAGACTCACAACCACCTGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGAAGGGACGAAcc                                                                                                                                                                                                          >  1:302641/1‑251 (MQ=255)                        cAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGTCGTTTGGCTGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTTTGCTGGAGAACGCCGACGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAACCTCTGG‑TTGTTCAGCCGACTAGACGTCCGGCCGTGACGAAcc                                                                                                                                                                                                          >  1:359941/1‑251 (MQ=255)                                                                  ccGAATAAATAGTGGAGACGTGTAGATGGGTAAAATAATTGGTATCGTCCTTGGTACTACCAACTCTTGTGGAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGGCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAACCGCCTGATTGgccgcc                                                                                                                                                               <  2:132076/251‑1 (MQ=255)                                                                               ggAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACg                                                                                                                                                   <  1:179098/250‑1 (MQ=255)                                                                                gagaCGTTTAGATGGGCAAAATAAATGGTATCGCCCTGGGTACTACCACCTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGATGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGACTATACCCAGGATGGTGAAACTCTAG‑TTGGTCAGACGGCTAAACGTCCGGCAGTGACGAACCAG‑AAAACAACCCTGTTTGCGATTCAACGCCTGATTGGACACCGCTTCCAGGACga                                                                                                                                                  >  2:219083/1‑250 (MQ=255)                                                                                             gggTAAAATAATTGGTATAGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTg                                                                                                                                    >  1:67164/1‑251 (MQ=255)                                                                                                        ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTATACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGCTTGGCCGCCGATTCCAGGACGGATAAGTActggt                                                                                                                                      >  2:5340/1‑234 (MQ=255)                                                                                                        ttGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGTGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCg                                                                                                                                      <  1:5340/238‑1 (MQ=255)                                                                                                                               caCTCGTGTGTAGCGATTATGGATTGCAACACTTCTCGTGTGCTGGAGAACGTCGTAGGCGATCGAATCTTGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCTGCTAAACGTCAGGCAGTGACGAACCCGCAAATCA‑CCCTGTTTGCTATTAAATGCCTGATTGGCAGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCCTCATGCCGTTCAAAATTATtgctgc                                                                                                  <  1:277433/250‑1 (MQ=255)                                                                                                                                                   ggatggCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATgg                                                                               <  1:48756/250‑1 (MQ=255)                                                                                                                                                                     tgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                              >  1:381721/2‑217 (MQ=255)                                                                                                                                                                     tgTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGAt                                                                                              <  2:381721/216‑1 (MQ=255)                                                                                                                                                                                   ccGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAACATTATTGCTGCTGATAACGGCGACGCATGGGCCGACGTTAAAGTCCAGAAAATGGCAccgccg                                              >  2:134485/1‑251 (MQ=255)                                                                                                                                                                                               gCACCACGCCTTCTATAATTGCCTATACCCAGGATGGTGAAACTCTGA‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGAAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTTCAGCGTGATGTTTCCATCATGCCGTTCAAACTTATTGCTGCTGAGAACGGAGACGCATGGTTCGCACTTAAAGGCCAGACAATGGCACCGCGGCAAATCTctgct                                  >  2:118430/1‑251 (MQ=255)                                                                                                                                                                                                     cGCCTTCTATCATTGCATATACCCAGGATGGTGTAACTCTGG‑TTGGTCAGCAGGCTAAACGTCAGGCAGTGACGAACCCGAAAAACA‑CCATGTTTGCGATTCACCGCATGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATAATGCCGTTCAAAATTATTGCTGCTGCTAACGGCGACGCATGGGTCGACGTTAACGGCCAGAAAATGGCACCGCCGCAGATCTCTGCTGAAg                              >  1:375766/1‑249 (MQ=255)                                                                                                                                                                                                                             aGGATGGTGAAAC‑CTGGGTTGGTCAGCCGGCAAAACGTCAGGCATGGACGACCCCGCAAAACC‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCACCTTGCCGTACAAAATTATGGCTGCTGATAACGGCGACGCAGGGGTCGAAGTTAAAGGCAAGAAAATGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCt    <  1:168419/251‑1 (MQ=255)                                                                                                                                                                                                                                aTGGAGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATATCTGCTGAAGAGCTGAAAAAAATGAAGAAAACCGCTGa  >  1:151618/1‑250 (MQ=255)                                                                                                                                                                                                                                   gTGAAACTCTGG‑TTGGTCAGCCGGCTAAACGTAAGGAAGTGACGAACCCGCAACACA‑CCCTGTTTGCGATTAAACGCCTGATTGGCCGCCGCTTCCAGGACGAAAAAGAAAAGCCTGACGTTTCCATCATGCAGTTAACAATTATTGCTGCTGATAACGGAGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCTGCAGATCTCTGATGAAGTGCGGAAAAAAATGAAGAAAACCGCTGa  >  1:368451/1‑247 (MQ=255)                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                               AAATGAAATTGGGCAGTTGAAACCAGACGTTTCGCCCCTATTACAGACTCACAACCACATGATGACCGAATATATAGTGGAGACGTTTAGATGGGTAAAATAATTGGTATCGACCTGGGTACTACCAACTCTTGTGTAGCGATTATGGATGGCACCACTCCTCGCGTGCTGGAGAACGCCGAAGGCGATCGCACCACGCCTTCTATCATTGCCTATACCCAGGATGGTGAAACTCTAG‑TTGGTCAGCCGGCTAAACGTCAGGCAGTGACGAACCCGCAAAACA‑CTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGA  >  NC_000913/12073‑12545

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