Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,065,909	G→C	8.2%	A15P (GCC→CCC)	PP_2678 →	hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,065,909	0	G	C	8.2%	147.5 / 5.1	49	A15P (GCC→CCC)	PP_2678	hydrolase
Reads supporting (aligned to +/- strand):  ref base G (24/21);  new base C (2/2);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

GCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAGCGGTGTGGTGGT  >  NC_002947/3065768‑3066057                                                                                                                                               |                                                                                                                                                       gCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGCCGGCAGCCCCAGGc                                                                                                                                                 <  6:125941/149‑1 (MQ=255)     acacGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGAt                                                                                                                                             >  4:215220/1‑149 (MQ=255)       acGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCt                                                                                                                                           >  6:284132/1‑149 (MQ=255)                 ttGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCgg                                                                                                                                 >  8:131642/1‑149 (MQ=255)                    cGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCt                                                                                                                               >  5:57636/1‑148 (MQ=255)                    cGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCt                                                                                                                               <  3:246417/148‑1 (MQ=255)                    cGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTACCGGCCCAGGCTGATCTGGCGTACCGGCt                                                                                                                               <  6:351734/148‑1 (MQ=255)                            caGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCAc                                                                                                                            <  8:314354/143‑1 (MQ=255)                            caGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCAc                                                                                                                            >  7:314354/1‑143 (MQ=255)                            caGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCg                                                                                                                      >  1:352632/1‑149 (MQ=255)                               cGGCTACCAGTTCGAGGCGGCGCTGCCGTGAGCGCGTTGTCTGATTTCACACTGTCCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GTGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCa                                                                                                                   <  6:329899/149‑1 (MQ=255)                                 gCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGa                                                                                                                 <  7:131642/149‑1 (MQ=255)                                   aaCCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCGCCACCCgg                                                                                                                     >  1:140646/1‑143 (MQ=255)                                        gTTCGCGGCGGCGCTGCCCTGAGCGCGTTGGCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGa                                                                                                          <  4:261892/149‑1 (MQ=255)                                         ttCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAg                                                                                                         <  3:215220/149‑1 (MQ=255)                                                    gcTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGtt                                                                                              <  3:285314/149‑1 (MQ=255)                                                     cTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTg                                                                                             <  6:57636/149‑1 (MQ=255)                                                             tgcgcGTTGTCTGATTTCACCCTGAGCCGAGGCACCGTCATGCGTTGTATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAg                                                                                      <  2:140646/147‑1 (MQ=255)                                                                           tttCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGa                                                                                                          >  1:203998/1‑114 (MQ=255)                                                                           tttCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGa                                                                                                          <  2:203998/114‑1 (MQ=255)                                                                            ttCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGgca                                                                                   >  6:351268/1‑136 (MQ=255)                                                                            ttCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGgca                                                                                   <  5:351268/136‑1 (MQ=255)                                                                            ttCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑ttt                                                                    >  7:88479/1‑149 (MQ=255)                                                                              cACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTc                                                                   <  1:55900/148‑1 (MQ=255)                                                                                 ccTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCTGGGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGc                                                                 >  1:315242/1‑148 (MQ=255)                                                                                 ccTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGcc                                                                >  2:186043/1‑148 (MQ=255)                                                                                  cTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCAcccccc                                                                                                                       >  1:336051/1‑94 (MQ=255)                                                                                  cTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCAcccccc                                                                                                                       <  2:336051/94‑1 (MQ=255)                                                                                   tGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAgg                                                                                     <  3:83000/127‑1 (MQ=255)                                                                                   tGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAgg                                                                                     >  4:83000/1‑127 (MQ=255)                                                                                   tGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAgg                                                                                     <  7:185795/127‑1 (MQ=255)                                                                                   tGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAgg                                                                                     >  8:185795/1‑127 (MQ=255)                                                                                        cGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCaa                                                        >  3:325088/1‑149 (MQ=255)                                                                                           ggCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATgg                                                     <  2:315242/149‑1 (MQ=255)                                                                                                     tGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAAtgtg                                           <  2:352632/149‑1 (MQ=255)                                                                                                        gTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGt                                         >  1:401361/1‑148 (MQ=255)                                                                                                              tCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACAcc                                  >  6:210313/1‑149 (MQ=255)                                                                                                              tCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACAcc                                  >  7:252382/1‑149 (MQ=255)                                                                                                                     gTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCa                           >  1:15819/1‑149 (MQ=255)                                                                                                                            ctgTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGTCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCa                     <  5:284132/148‑1 (MQ=255)                                                                                                                            ctgTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTc                      <  2:15819/147‑1 (MQ=255)                                                                                                                              gtgCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGAGAGCACCGATGAGGTTTCGCCAAGGCCAATGGGGGCAATGTGGTCAACACCGCGTTCATCGTCAc                    >  1:334924/1‑149 (MQ=255)                                                                                                                                gCCT‑GGGCCGGCAGCCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAaca                                    >  4:23776/1‑129 (MQ=255)                                                                                                                                gCCT‑GGGCCGGCAGCCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAaca                                    <  3:23776/129‑1 (MQ=255)                                                                                                                                gCCT‑GGGCCGGCAGCCCCAGGCTGATCTGGCGTACCGGCTCAACCCCCGGCAGATCGCCGAGGGCACCTCGTTGCTGGAAGGCAGCACAGACAA‑‑TTTCGACAAGGCCAATGGCGGCAATGTGGTCATCACGGCGTTCATCGTCAccaa                 >  7:65242/1‑146 (MQ=255)                                                                                                                                  cT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAg              >  7:151066/1‑149 (MQ=255)                                                                                                                                     ggggCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGATGGCAGCACCGCCAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAGCg            >  3:27830/1‑149 (MQ=255)                                                                                                                                            ccGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGGAGATCGGCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCTACACCGCGTTCATCGTCACCGTCTGCTTTgtgg      <  4:325088/148‑1 (MQ=255)                                                                                                                                            ccGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAGCGGTgtggt     >  7:85087/1‑149 (MQ=255)                                                                                                                                               gCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAGCGGTgtggtggt  >  5:396875/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                       GCCAACACGGCTGCTGTTGCGGGACAACAGCGGCAACCAGTTCGAGGCGGCGCTGCCCTGAGCGCGTTGTCTGATTTCACCCTGACCCGAGGCACCGTCATGCGTTGGATCCTGCTGTTGCTGCTGTGCCT‑GGGGCTGCCGGCCCAGGCTGATCTGGCGTACCGGCTCACCCCCCGGCAGATCGCCGAGGGCACCTGGTTGCTGGAAGGCAGCACCGACAA‑‑TTTCGCCAAGGCCAATGGCGGCAATGTGGTCAACACCGCGTTCATCGTCACCGACAGCGGTGTGGTGGT  >  NC_002947/3065768‑3066057

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