Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,656,546	T→A	18.0%	K65N (AAA→AAT)	hdeA ←	periplasmic acid stress chaperone HdeA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,656,546	0	T	A	18.0%	50.0 / 3.2	22	K65N (AAA→AAT)	hdeA	periplasmic acid stress chaperone HdeA
Reads supporting (aligned to +/- strand):  ref base T (9/9);  new base A (2/2);  total (11/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACAGGAGGAATAGCTGCTTTTAAAGAAGATATTCCGGCCTGTTCAATATGCAAGGAAGTACGATGTAACGCATTCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATATCCTCAACTATAAAGTGAAAGAGCCGTCACGAATCAATTTCGACACT  >  NC_000913/3656326‑3656794                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         aaCAGGAGGAATAGCTGCTTTTAAAGAAGATATTCCGGCCTGTTCAATATGCAAGGAAGTACGATGTAACGCATTCCCGGAATTACATATATTTCTTAATTTTGTCCCATGCGCCCTTAACTTTATCTTTATAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGAt                                                                                                                                                                                                    >  1:214674/1‑275 (MQ=255)                                                                      cGCATTCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGc                                                                                                                                                                                                             >  2:182748/1‑197 (MQ=255)                                                                          ttCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCTCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGa                                                                                                                                                                                                         <  1:182748/197‑1 (MQ=255)                                                                           tCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGATTATCTttgt                                                                                                                                                                                                                                                >  1:17763/1‑157 (MQ=255)                                                                           tCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCAACTTCTGGATTATCTttgt                                                                                                                                                                                                                                                >  1:17768/1‑157 (MQ=255)                                                                               ggAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGATTAGCTTTgttgtt                                                                                                                                                                                                                                            <  2:17763/157‑1 (MQ=255)                                                                               ggAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGATTAGCTTTgttgtt                                                                                                                                                                                                                                            <  2:17768/157‑1 (MQ=255)                                                                                                                                tttAAAGTTGGGTTGTTTATCCTTAGTACAAGCCTGAACGATAGCTGGGGGTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTcac                                                                                                                                                                           <  1:33238/173‑1 (MQ=255)                                                                                                                                                             caGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAAGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCa                                                                                                                            >  2:105458/2‑191 (MQ=255)                                                                                                                                                                 cTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAAGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCg                                                                                                                        <  1:105458/191‑1 (MQ=255)                                                                                                                                                                    aaCGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACtt                                                               >  2:165687/1‑245 (MQ=255)                                                                                                                                                                        aTAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACtttttt                                                           <  1:165687/245‑1 (MQ=255)                                                                                                                                                                         tAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGt                                                                                                                                                               >  1:15884/1‑144 (MQ=255)                                                                                                                                                                             tGGGGTTCCGGTTGCAATACCCTGAACATCTGAAACCGCATCTTCTTGTTTATCTTTGTTGTTCAGCGCTTCTGCAAAACCAACTGCAGTTTGCTGGAAGGATTCTTCCACAGCCAGGAAATCTTCACAGGTCCAGGGGTTGAc                                                                                                                                                           <  2:15884/144‑1 (MQ=255)                                                                                                                                                                                          gCAGTACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGc                                                                                                                       >  2:136386/1‑167 (MQ=255)                                                                                                                                                                                          gCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATATCCTCAACTATAAAGTGAAAGAGCCGTCACGAATCAAt           >  1:41094/1‑275 (MQ=255)                                                                                                                                                                                              tACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCa                                                                                                                   <  1:136386/167‑1 (MQ=255)                                                                                                                                                                                                   caaaCATCAAAAACCTCTTCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAAGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCTCAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATATCCTCAACTATAAAGTGAAAGAGCCGTCACGAATCAATTTCGACACt  <  2:41094/273‑1 (MQ=255)                                                                                                                                                                                                     aaCATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTg                                                                                                      >  1:259129/1‑173 (MQ=255)                                                                                                                                                                                                         gcTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTGTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCag                                                                                                  <  2:259129/172‑1 (MQ=255)                                                                                                                                                                                                          cTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTTTCAGCTGCttttt                                                                                                                                 >  1:153387/1‑141 (MQ=255)                                                                                                                                                                                                              aaCCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTgcgc                                                                                                                             <  2:153387/141‑1 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                         AACAGGAGGAATAGCTGCTTTTAAAGAAGATATTCCGGCCTGTTCAATATGCAAGGAAGTACGATGTAACGCATTCCCGGAATTACATATCTTTCTTAATTTTGTCCCATTCGCCTTTAACTTTATCTTTAAAGTTGGCTTGTTTATCCTGAGTACAAGCCTGAACGATAGCTGGGGTTACGGTTGCAATACCCTGAACATCTAAAACCGCATCTTCTGGTTTATCTTTGTTGTTCAGCGCTTCAGCAAAACCAACTGCAGTTGGCTGGAAGGATTCGTCCACAGCCAGGAAATCTTCACAGGTCCAGGAGTTGACCGGTTTTTTGTTATCAGCTGCTTTTTGCGCATCCGCTGCATTGCTCACAACTGGCAGAAGAAGCAGACCACCAAGAATAACGCCTAATACTTTTTTCATCGTAATATCCTCAACTATAAAGTGAAAGAGCCGTCACGAATCAATTTCGACACT  >  NC_000913/3656326‑3656794

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