Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	12,714	C→T	G184G (GGC→GGT)	dnaK →	chaperone protein DnaK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	12,714	0	C	T	93.3%	77.1 / ‑5.9	31	G184G (GGC→GGT)	dnaK	chaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (15/13);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGA  >  NC_000913/12488‑12953                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aaaTGGCACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGtt                                                                                                                                                                                                                          >  2:83185/1‑251 (MQ=255)   aTGGCACAGCCGCAGATCTCTGCTGAAGTGCTGAAAAACATGACGAAAACCGCTGAAGTTTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGtttt                                                                                                                                                                                                                        >  2:32130/1‑250 (MQ=255)       cACCGCCGCAGATCTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAGGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTAAGGTCTGGACAAAGGTACTGGCAACCGTACTACCGCGGTTTATGAc                                                                                                                                                                                                                    <  2:49087/251‑1 (MQ=255)             cgcagaactCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGggt                                                                                                                                                                                                              <  2:56033/243‑1 (MQ=255)                       tgctgAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATCCTTTAACGATGCTCCGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACtt                                                                                                                                                                                                    >  1:43923/1‑251 (MQ=255)                          tgAAGCGCTGAAAAAAATGAAGAAAACAGCTGAAGATTACCTGGGTGAACCGGAAACTGCCGCTGTTCTCACCGTACCGGCATACTTTAACGCTGATAAGCGTCAGGCAACCAAAGACGCAGGCCGTCTCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGCCCTGGGTGGTGGTACTTTCg                                                                                                                                                                                                 >  2:268042/1‑251 (MQ=255)                            aaGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGat                                                                                                                                                                                               >  1:34885/1‑251 (MQ=255)                                    aaaaaaaTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTTGGTGGTGGTACTTTCGATATTTCtat                                                                                                                                                                                       >  1:200577/1‑251 (MQ=255)                                       aaaaTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTTGGTGGTGGTACTTTCGATAGTTCtattat                                                                                                                                                                                    <  1:57301/251‑1 (MQ=255)                                         aaTGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGCAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGAGTCAGGAAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATAATCAACGAACCAACCGCAGCTGCGCTGGCTTACGGTCTGGACCAAGGTACTGGCAACCATACTATCGCGGTTTATGACATGGGTGGAGGTACTTTCGATATTTCTATTATCg                                                                                                                                                                                  >  2:142385/1‑251 (MQ=255)                                                 aaaCCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGACCCGACCGCAGCTGCGCTGGCTTACGGCCTGGACAACGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATcgac                                                                                                                                                                          >  1:117537/1‑251 (MQ=255)                                                            gATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACgg                                                                                                                                                               >  1:6393/1‑251 (MQ=255)                                                               tACCTGGGTGAACCGGTAAATGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCg                                                                                                                                                             >  1:100473/1‑250 (MQ=255)                                                                         aaCCGGTAACTGAAGCTGTTATCACCGTACCGGCAAACTTTAACGATTCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGTTCTGGAAGTAAAACGTATCATCAACGAACCTACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCTGTTGATGACCTGGGTGGTGGTACTGTGGATATTTCTATTATCTAAATCGACGAAGTTGACGGCGAAAAAACCTTc                                                                                                                                                  <  1:98714/251‑1 (MQ=255)                                                                                               cACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGt                                                                                                                             <  1:46016/250‑1 (MQ=255)                                                                                               cACCGTACCGGCAAACTTTAACGATGCGCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGCCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAATGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCATCGGTg                                                                                                                            <  2:140939/251‑1 (MQ=255)                                                                                                    tACCGGCATCATTTAACGAAGCTCAGCGTCAGGCAACCACAGCCGCAAGCCGTATAGCTGGCCTGGAAGTAAAACGTATAATAAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTAGAAGTTCTGGAAACCAACGGTGATAcc                                                                                                                       >  2:41783/1‑251 (MQ=255)                                                                                                                               gTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGc                                                                                            >  2:57876/1‑251 (MQ=255)                                                                                                                                                               ggTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTAGCGAAATCGACGAAGTTGACGGCGAAAAAACCGTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTGCGACAGCCGTCTGATCATCTACC‑TGGTTGAAGAATTc                                                             <  2:43923/249‑1 (MQ=255)                                                                                                                                                                    ggAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGCTATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCaagaa                                                        >  1:160212/1‑249 (MQ=255)                                                                                                                                                                          aaaaCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAg                                                 <  1:83185/250‑1 (MQ=255)                                                                                                                                                                             aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGc                                              >  2:41739/1‑250 (MQ=255)                                                                                                                                                                             aCGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGc                                              >  1:146876/1‑250 (MQ=255)                                                                                                                                                                              cGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCAt                                            <  1:32130/251‑1 (MQ=255)                                                                                                                                                                              cGTATCATCAACGAACCGACCGCAGCGGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTGTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCAt                                            <  1:106795/251‑1 (MQ=255)                                                                                                                                                                                          gAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCTGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTAGCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTGCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTg                                                                     <  1:116004/214‑1 (MQ=255)                                                                                                                                                                                          gAACCGACCGCAGCTGCGCTGGCTTACGGCCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGCACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGCGGTGAAGACTTCGACAGCCGTCTGCTCAACTACC‑TGGTTg                                                                     >  2:116004/1‑214 (MQ=255)                                                                                                                                                                                           aaCCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCTCGTTTTATGACCTGGGTGGTTGTACTTTCGATATTTCTATTATCGAAATCGTCGATGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGTCTTCGACAGCCGTCTGATCAACT‑TCATGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAAc                               <  1:35414/251‑1 (MQ=255)                                                                                                                                                                                                                cTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTgc                                    >  2:230724/1‑225 (MQ=255)                                                                                                                                                                                                                cTTACGGTCTGGACAAAGGTACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTgc                                    <  1:230724/225‑1 (MQ=255)                                                                                                                                                                                                                          ggACAAGGGGACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGAGTGGTGAAGACTTCGACAGCCGTCTGATCAACTACC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATACGCTGGCAATGCAGCGCCTGAaaga  >  2:230247/1‑249 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 AAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGCAGGCCGTATCGCTGGTCTGGAAGTAAAACGTATCATCAACGAACCGACCGCAGCTGCGCTGGCTTACGGTCTGGACAAAGGCACTGGCAACCGTACTATCGCGGTTTATGACCTGGGTGGTGGTACTTTCGATATTTCTATTATCGAAATCGACGAAGTTGACGGCGAAAAAACCTTCGAAGTTCTGGCAACCAACGGTGATACCCACCTGGGGGGTGAAGACTTCGACAGCCGTCTGATCAACTATC‑TGGTTGAAGAATTCAAGAAAGATCAGGGCATTGACCTGCGCAACGATCCGCTGGCAATGCAGCGCCTGAAAGA  >  NC_000913/12488‑12953

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