Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,559	T→G	17.9%	K213Q (AAG→CAG)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,559	0	T	G	17.9%	54.8 / 8.7	28	K213Q (AAG→CAG)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (13/10);  new base G (3/2);  total (16/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 = 9.70e-01

CCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACC  >  NC_002947/626424‑626707                                                                                                                                        |                                                                                                                                                       ccGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                          <  2:334640/150‑1 (MQ=255) ccGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGccaa                                                                                                                                          >  1:552640/1‑148 (MQ=255)  cGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAg                                                                                                                                         >  1:320226/1‑150 (MQ=255)   gTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAAgg                                                                                                                                        <  2:552640/150‑1 (MQ=255)    tCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGTTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCGGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGgc                                                                                                                                       <  2:260975/150‑1 (MQ=255)                 gCCGCTCGCGTACGGCCTGCCACTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCCGGCGCGGACTGACCCCAGGCCGTTCGGTCTTTCCCTGTTCAGTGATCCTTGGCCGTACGGCCCAGGCGCCCGCgcagcat                                                                                                                          <  2:257924/150‑2 (MQ=255)                             cGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcat                                                                                                              <  2:320226/150‑1 (MQ=255)                              ggCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                             <  2:338546/150‑1 (MQ=255)                                 cTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGATCCTGGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGc                                                                                                          >  1:496731/1‑150 (MQ=255)                                             tgGCGTGATCGCGTAACGCGCCAGATTACGCTTGCGCCGCCAGGGGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACCCGGCCAGCCGAGAGTCAGCGACACCAGGCACAgg                                                                                              >  2:227343/1‑150 (MQ=255)                                                       gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGccaa                                                                                                                                          <  2:468404/96‑3 (MQ=255)                                                       gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGccaa                                                                                                                                          >  1:468404/1‑94 (MQ=255)                                                       gggTAACGCGCAAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGc                                                                                    >  2:55733/1‑150 (MQ=255)                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                             <  1:380169/116‑1 (MQ=255)                                                                gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAatcatc                                                                                                             >  2:380169/1‑116 (MQ=255)                                                                   agagTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCTGCGACAGGATGATGATGCCCACCACCATCa                                                                        <  2:496731/150‑1 (MQ=255)                                                                       tACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGa                                                                    <  1:298873/150‑1 (MQ=255)                                                                          gCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAg                                                                                                                         <  2:33608/94‑1 (MQ=255)                                                                          gCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAg                                                                                                                         >  1:33608/1‑94 (MQ=255)                                                                                                            gttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGaagaag                             >  1:237407/1‑150 (MQ=255)                                                                                                                       aCCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGt                  >  2:430258/1‑150 (MQ=255)                                                                                                                         cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                >  2:217973/1‑150 (MQ=255)                                                                                                                         cTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGa                >  2:284344/1‑150 (MQ=255)                                                                                                                          tGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGAc               <  2:307410/150‑1 (MQ=255)                                                                                                                            ttCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa             >  2:484302/1‑150 (MQ=255)                                                                                                                               aGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgc          >  1:37304/1‑150 (MQ=255)                                                                                                                                gTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgcc         >  1:546716/1‑150 (MQ=255)                                                                                                                                      cTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACGATCAGCGACAGGTGATGATCGACTGCGCTTGCACG‑‑GCTGAATAAGTAGCCCAAGGTGACCAGGCCGCCGAc   >  1:442053/1‑150 (MQ=255)                                                                                                                                       ttGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGAcc  <  2:117111/150‑1 (MQ=255)                                                                                                                                        |                                                                                                                                                       CCGTCCAGCATTGGCAGCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCAGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGA‑TG‑AATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCCGACC  >  NC_002947/626424‑626707

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