Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	6,239	G→T	A74E (GCG→GAG)	yaaA ←	DNA binding and peroxide stress response protein YaaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	6,239	0	G	T	96.3%	67.6 / NA	27	A74E (GCG→GAG)	yaaA	DNA binding and peroxide stress response protein YaaA
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (17/9);  total (17/10)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGT  >  NC_000913/5999‑6466                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     ttgttgGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTgc                                                                                                                                                                                                                             >  2:159905/1‑249 (MQ=255)      ggTGATGAGATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGGCGATAAGGCTGCATTAACTCGAGCGGGCGCAGTACGCCTTACAATCCGTAATGCATTCGCAAGTGCTGTTGGGCAAAAGCGAAGGCGTTTTCGCGGAAGGTTTCGGCCTGTAAGCCGTTGTAGACATGACCGGTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTg                                                                                                                                                                                                                      <  1:240420/251‑1 (MQ=255)         gattatatCTCCCCAGAATTGATACAGATCTTTCCCTCGTGCATTCTCACGACGGATCCCCATTCCCAGACGATAAGGCTGCAGTAAAGCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTTAAGGTTTCGGCCTGCAAGCCGTTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTTTCCGGCGTGAAA‑TCTGGc                                                                                                                                                                                                                    <  2:81947/246‑1 (MQ=255)              tatCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCCAGACGGATCACCATTTCCAGACGATAAAGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCCTTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAAGTTTCGGCCTGCACGACGGTGTAGCCATCACCTTTATACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGCGAAA‑TCTAGATGCCa                                                                                                                                                                                                               >  1:153797/1‑250 (MQ=255)                 ctcCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCCAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTc                                                                                                                                                                                                            >  2:21765/1‑250 (MQ=255)                                                     tctcAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATAGAGCGGGCGCAGTCCGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAGATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGCATCGCCTGGCGGGCATTCTCCGGAGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGGTGATACCCGCCAGTTTGTCg                                                                                                                                                                       >  1:117290/1‑251 (MQ=255)                                                         aaGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCCGTAAGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGACA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTg                                                                                                                                                                    >  2:138093/1‑250 (MQ=255)                                                            ccGGATCCCCTTTTCCAGACTATAAGGCTGCTTTAAATCGATCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCGGTTATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATgc                                                                                                                                                                <  1:161699/250‑1 (MQ=255)                                                                aTCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGCAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCAt                                                                                                                                                            >  1:131583/1‑251 (MQ=255)                                                                      aTTTCCAGACGATAAGGTTGCATTAAATCGAGGGGGCGCAGAACGCCATACAAGCCGGAAAGCATTCGCAAAGGCGGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGGAGACATCTCCTTGAAACGCCAGAATCGCCGGGCGGGCATTCTCCGGCGGGAAA‑TCGGGTTGCCAGTAAAGAAAGTGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGt                                                                                                                                                      <  1:173242/251‑1 (MQ=255)                                                                                                   gagggggcgcAGTCCGCCAAACAAGCCGGAAAGCTTTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTGCCGGCCGGCAGGCCGGTGAAGACATCACTTTTAAACGCCAGAAGCGCGTGGCGGGCATTCGCCGGCGTGAAA‑TCTGGCTGCCAGTCAGGAAAGCGGGCGGCGTTGATACCCGCCAGGTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCcgcg                                                                                                                         <  2:121161/247‑1 (MQ=255)                                                                                                       gggcgcAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCgcgcttc                                                                                                                     >  1:191617/1‑248 (MQ=255)                                                                                                            caGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATgg                                                                                                                 >  1:103399/2‑250 (MQ=255)                                                                                                                  gCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCCCCTTTACACGCCCGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGGCGCTGGTGCGCATCAGCGTGTTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGGTCAAct                                                                                                          >  2:71937/1‑251 (MQ=255)                                                                                                                          ggCCGGAAAGCATTCGAAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATCCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGaa                                                                                                  <  1:196181/250‑1 (MQ=255)                                                                                                                              ggAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGt                                                                                              >  1:194376/1‑251 (MQ=255)                                                                                                                                                    ggggaaaaaTCGAAATCTTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTTTAGACATCACCTTTAAACGCCAGAATCGCCGGGCGGGCATTCTCCGGCGGGAAA‑TCTGGCTGCAAGTTATGAAAGCGAGCGGCGTTGATTCCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGtat                                                                         <  1:265710/245‑1 (MQ=255)                                                                                                                                                     gggCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGtata                                                                        <  1:21765/250‑1 (MQ=255)                                                                                                                                                       gCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAgcg                                                                     >  2:101905/1‑251 (MQ=255)                                                                                                                                                                                gggTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCAGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGATTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTCACAGCTCCTGCAGCGTATAGCGCGTGGTGGTCAACGGGTTTTGGTa                                             >  2:26587/2‑250 (MQ=255)                                                                                                                                                                                       ggCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCg                                     >  1:176321/1‑251 (MQ=255)                                                                                                                                                                                           tGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGtttt                                 >  2:86065/1‑251 (MQ=255)                                                                                                                                                                                              aaGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGc                              >  2:77783/1‑251 (MQ=255)                                                                                                                                                                                                  cGGTGTATACAGCCCCTTTAACGGCCATAATCGCCTGGTGGGGATTCTCCGGC‑TGAAAATCTGGCTGCGAGTCAGGAAAGTGATCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACGGCTGGGAATTGTCTAACAGCTCCGGCAGCGTAGAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGt                          <  1:235174/251‑1 (MQ=255)                                                                                                                                                                                                          aCATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAC‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGCGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGaacta                    >  2:232199/1‑246 (MQ=255)                                                                                                                                                                                                                cacttAAACGCCCGAATCGCCGGGCGGGCTTTCTCCGGCGTAAAA‑GCTGGCTGCCATTCATGAAAGCGAGCGGGGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAg            <  2:164411/248‑1 (MQ=255)                                                                                                                                                                                                                          ccAGAATCGCCTGGCGGGCATTCTCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCTTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGt  >  1:6223/1‑251 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     TTGTTGGTGATGATATCTCCCCAGAATTGATACAGATCTTTCCCTCGGGCATTCTCAAGACGGATCCCCATTTCCAGACGATAAGGCTGCATTAAATCGAGCGGGCGGAGTACGCCATACAAGCCGGAAAGCATTCGCAAATGCTGTTGGGCAAAATCGAAATCGTCTTCGCTGAAGGTTTCGGCCTGCAAGCCGGTGTAGACATCACCTTTAAACGCCAGAATCGCCTGGCGGGCATTCGCCGGCGTGAAA‑TCTGGCTGCCAGTCATGAAAGCGAGCGGCGTTGATACCCGCCAGTTTGTCGCTGATGCGCATCAGCGTGCTAATCTGCGGAGGCGTCAGTTTCCGCGCCTCATGGATCAACTGCTGGGAATTGTCTAACAGCTCCGGCAGCGTATAGCGCGTGGTGGTCAACGGGCTTTGGTAATCAAGCGTTTTCGCAGGTGAAATAAGAATCAGCATATCCAGT  >  NC_000913/5999‑6466

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: ‑