Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,921,591	T→A	17.2%	intergenic (+24/+37)	PP_1720 → / ← PP_1721	alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,921,591	0	T	A	17.2%	60.5 / 3.6	28	intergenic (+24/+37)	PP_1720/PP_1721	alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (12/11);  new base A (3/2);  total (15/13)
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

GGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921442‑1921739                                                                                                                                                      |                                                                                                                                                     gggCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGtt                                                                                                                                                      <  2:451410/150‑1 (MQ=255)       ttGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTCCGCGGGGCTGGCCCCTTCGCGGGTATACCcg                                                                                                                                                <  1:438955/150‑1 (MQ=255)          tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGTTACCCGCGa                                                                                                                                             >  2:182057/1‑150 (MQ=255)              cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGGGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                         <  2:358161/150‑1 (MQ=255)                       tCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGcaca                                                                                                                                >  1:122438/1‑150 (MQ=255)                              cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATa                                                                                                                         >  1:174239/1‑150 (MQ=255)                                      tGGAGTTTGCCCAGGCGTCGTCTGCTGATTTAGAGTCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:494732/150‑1 (MQ=255)                                      tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:111598/150‑1 (MQ=255)                                              gCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctcca                                                                                                         <  1:120170/150‑1 (MQ=255)                                               cccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:493557/141‑1 (MQ=255)                                               cccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 >  1:493557/1‑141 (MQ=255)                                                           aTGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGc                                                                                            <  2:551944/150‑1 (MQ=255)                                                            tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGcc                                                                                           >  2:488528/1‑150 (MQ=255)                                                                   cTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAtt                                                                                    <  1:488528/150‑1 (MQ=255)                                                                    tGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCGGCTCAGAGTCGGCCACCCATTg                                                                                   >  1:416463/1‑150 (MQ=255)                                                                           caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCgg                                                                            >  2:267084/1‑150 (MQ=255)                                                                                              cTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGc                                                         <  2:174239/150‑1 (MQ=255)                                                                                                                       tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTcc                                                     >  2:253513/1‑129 (MQ=255)                                                                                                                       tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTcc                                                     <  1:253513/129‑1 (MQ=255)                                                                                                                         gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                                  >  2:382034/1‑130 (MQ=255)                                                                                                                         gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCggg                                                  <  1:382034/130‑1 (MQ=255)                                                                                                                                        gCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTAAGCTGCAACTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTccacaacaa               >  2:77690/1‑150 (MQ=255)                                                                                                                                           ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg            >  2:63214/1‑150 (MQ=255)                                                                                                                                            cTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCgg                                                   >  2:263711/1‑110 (MQ=255)                                                                                                                                            cTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCgg                                                   <  1:263711/110‑1 (MQ=255)                                                                                                                                            cTTCGCGGGTAAGCCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCACAGTCAGCCACCCATTGATCTCGGCATGCGCGCGCAACACTGCGTCCGGGTCCACCACATGCAGCGGATCCACTCTCAACAACAACGGc           >  2:366895/1‑150 (MQ=255)                                                                                                                                            cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCCCATTCGGGTGATCCACTCTCAACAACACCGGc           >  1:541297/1‑150 (MQ=255)                                                                                                                                             ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa          >  2:103605/1‑150 (MQ=255)                                                                                                                                              tCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCaa         <  1:267084/150‑1 (MQ=255)                                                                                                                                                     tttACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTc  <  1:182057/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GGGCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921442‑1921739

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