Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,372,554	G→A	V510V (GTG→GTA)	groL →	chaperonin GroEL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,372,554	0	G	A	96.7%	97.1 / NA	31	V510V (GTG→GTA)	groL	chaperonin GroEL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (17/12);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACC  >  NC_000913/4372425‑4372683                                                                                                                                  |                                                                                                                                  cacCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTAGTTGTGCCCTGCAGTAGGCAGGTTCTGtagct                                                                                                                                >  2:156843/1‑133 (MQ=255)   ccGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGtagctgg                                                                                                                              <  1:434130/133‑1 (MQ=255)    cGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGtagctggc                                                                                                                             <  2:308139/133‑1 (MQ=255)       tAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCtg                                                                                                                          <  1:126986/133‑1 (MQ=255)       tAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCtg                                                                                                                          >  2:439911/1‑133 (MQ=255)       tAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCtg                                                                                                                          >  2:126986/1‑133 (MQ=255)              ggcgACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATca                                                                                                                   >  2:484479/1‑133 (MQ=255)              ggcgACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATca                                                                                                                   >  1:69109/1‑133 (MQ=255)                          tACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCa                                                                                                       <  1:273344/133‑1 (MQ=255)                          tACGGTTACAACGCACCAACCGAAGAATACTGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCa                                                                                                       <  1:189403/133‑1 (MQ=255)                           aCGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCAt                                                                                                      <  2:69109/133‑1 (MQ=255)                                     cGCAGCAACCGAAGATTAAGGCTACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGAc                                                                                            <  1:252723/133‑1 (MQ=255)                                       cagcaACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTCGCTGGCCTAATGATCAGCACCGAATCCATGATTACCGACCt                                                                                          >  2:196703/1‑133 (MQ=255)                                       cagcaACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCt                                                                                          >  2:371936/1‑133 (MQ=255)                                         gcaACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGc                                                                                        >  1:427091/1‑133 (MQ=255)                                                   aaTACGGCCACATGATCGGCATAGGTATCCCGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGa                                                                              <  1:152823/133‑1 (MQ=255)                                                             cATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCAGCCGAAAAACGATTCAGCTGAc                                                                    >  1:202081/1‑133 (MQ=255)                                                                         gggTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCgctgct                                                        <  1:352771/133‑1 (MQ=255)                                                                              tCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAACCGATGCAGCTGTCTTAGGCGCAGCTGGCgg                                                   >  1:185956/1‑133 (MQ=255)                                                                                cTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTa                                                 <  2:504242/133‑1 (MQ=255)                                                                                    aTCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATggg                                             >  1:357402/1‑133 (MQ=255)                                                                                      ccaaccaaAGAAACTTGTTCTGCTGTGCAGTAGGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGggcg                                           <  2:202081/133‑1 (MQ=255)                                                                                        aaccaaAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGggcggc                                         >  1:289169/1‑133 (MQ=255)                                                                                         accaaAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGGGCTGCTGGCGGTATGGGCGGCa                                        >  1:186825/1‑133 (MQ=255)                                                                                         accaaAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCa                                        >  2:406014/1‑133 (MQ=255)                                                                                                      gTTCTGCTCTGCAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATggg                           <  1:406014/133‑1 (MQ=255)                                                                                                       ttCTGCTCTGTAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGggc                          >  1:283009/1‑133 (MQ=255)                                                                                                        tctgctctgcAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGG‑‑GGCATGGGTGACATGggcggc                       >  1:416350/1‑133 (MQ=255)                                                                                                           gctctgtAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCa                      <  2:283009/133‑1 (MQ=255)                                                                                                                 cAGTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGATGGCGGTAAGGGGGGCGTGGGTGGCATGGGGGGCATGATGt                >  2:185926/1‑133 (MQ=255)                                                                                                                   gTACGCAGCTTCTGTGGCTGGCCTTATGTTCACCACCGAATTCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTaa              <  2:186825/133‑1 (MQ=255)                                                                                                                   gTACGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTaa              <  1:371936/133‑1 (MQ=255)                                                                                                                      cGCAGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGGATGGGTGGCATCGGACGCATGCTGTAATTg           >  1:243057/1‑133 (MQ=255)                                                                                                                         aGCTTCTGTAGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCGTGGTTGCGATGGGCGGCATGACGTAATTGccc        >  2:169824/1‑133 (MQ=255)                                                                                                                               tgtaGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCAcc  <  2:427091/133‑1 (MQ=255)                                                                                                                                  |                                                                                                                                  CACCGTTAAAGGCGGCGACGGCAACTACGGTTACAACGCAGCAACCGAAGAATACGGCAACATGATCGACATGGGTATCCTGGATCCAACCAAAGTAACTCGTTCTGCTCTGCAGTACGCAGCTTCTGTGGCTGGCCTGATGATCACCACCGAATGCATGGTTACCGACCTGCCGAAAAACGATGCAGCTGACTTAGGCGCTGCTGGCGGTATGGGCGGCATGGGTGGCATGGGCGGCATGATGTAATTGCCCTGCACC  >  NC_000913/4372425‑4372683

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