Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,835,980	G→C	10.4%	W156C (TGG→TGC)	PP_3387 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,835,980	0	G	C	10.4%	141.0 / 6.7	48	W156C (TGG→TGC)	PP_3387	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/24);  new base C (2/3);  total (21/27)
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.72e-01

CTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAACTTGCC  >  NC_002947/3835836‑3836124                                                                                                                                                  |                                                                                                                                                 cTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGcca                                                                                                                                              >  7:199242/1‑148 (MQ=255)            ccAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGCCACGGGGGCCCAc                                                                                                                                   <  8:199242/149‑1 (MQ=255)             cAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGCGCCACGGGGGCCCAcc                                                                                                                                  <  6:245292/149‑1 (MQ=255)              aaCATCGTCACGGGCACGGACGGCCTCTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCAccc                                                                                                                                 <  7:63120/149‑1 (MQ=255)                 aTCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGgcgc                                                                                                                                     >  6:208142/1‑142 (MQ=255)                 aTCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGgcgc                                                                                                                                     <  5:208142/142‑1 (MQ=255)                     tCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAg                                                                                                                          <  7:210367/149‑1 (MQ=255)                        cGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCgg                                                                                                                       <  7:15009/149‑1 (MQ=255)                         gggCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGt                                                                                                                      <  3:17561/149‑1 (MQ=255)                              cggacggCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTtgtg                                                                                                                  >  3:163338/1‑148 (MQ=255)                                    gCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAg                                                                                                           <  4:86295/149‑1 (MQ=255)                                      cACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGat                                                                                                         <  2:284169/149‑1 (MQ=255)                                              ccTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCcggcg                                                                                                 <  8:44224/149‑1 (MQ=255)                                               cTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACCAGTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCCCCCGAGCAAGCGGTGGTGGATGCAGATAACcggcg                                                                                                 >  1:84892/1‑149 (MQ=255)                                                 gTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTg                                                                                              >  5:227272/1‑149 (MQ=255)                                                      gcTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCCCC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAg                                                                                         >  7:319324/1‑149 (MQ=255)                                                       cTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGc                                                                                        <  4:140980/149‑1 (MQ=255)                                                        ttCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGcc                                                                                       >  1:188776/1‑149 (MQ=255)                                                        ttCCACACCGCCGAGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCAATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGc                                                                                        >  4:132902/1‑148 (MQ=255)                                                                  ccgcGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACg                                                                             >  2:219774/1‑149 (MQ=255)                                                                              tAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCa                                                                 <  2:84892/149‑1 (MQ=255)                                                                                   cTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCAc                                                            <  4:163338/149‑1 (MQ=255)                                                                                    tGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTtgtg                                                                                                                  <  1:348550/94‑1 (MQ=255)                                                                                    tGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTtgtg                                                                                                                  >  2:348550/1‑94 (MQ=255)                                                                                    tGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACGGGTGCGCCACGGGGGCCCAccag                                                                                                                                <  3:118285/80‑3 (MQ=37)                                                                                    tGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACGGGTGCGCCACGGGGGCCCAccag                                                                                                                                >  4:118285/1‑78 (MQ=37)                                                                                      gACAGGCACC‑GTGAGCGGCCTGGGGTTAGCGAATGGCTGGCCGCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACg                                                           <  8:319324/147‑1 (MQ=255)                                                                                         aGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGaa                                                      >  3:275981/1‑149 (MQ=255)                                                                                           gCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAAc                                                     >  8:239251/1‑148 (MQ=255)                                                                                             aCC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACg                                                    >  4:279646/1‑147 (MQ=255)                                                                                                        gCCTGGTGTTAGCGAAAGGCTGGCCCCCATCCATCACAGGTGGGCCCCCGTGGCGCAACCGAGCAAGCGGATGAGGAAGCAGATATCCGGCGGTGATCCGCCTCAATCAACGCACTGTCCGCCAATCACGTTCAACGAAGCGTGTTggt                                        >  7:14482/1‑148 (MQ=255)                                                                                                        gCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTggg                                        >  4:64253/1‑149 (MQ=255)                                                                                                          cTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGcc                                      <  7:239251/149‑1 (MQ=255)                                                                                                          cTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGcc                                      >  5:126610/1‑149 (MQ=255)                                                                                                             ggcttAGCGAATGGCTGGCCCCCAGCCATCACTTGTGGGCCCCCGTGGCGCACCCGAGCAAGGCGTTGTGCATGGAGAAATGCGGCGGTGATCATCCTCAATCCACGCACTGTCCGCTAATCACGTTGAACGAAGCGTGATGTGCCg                                     <  8:14482/144‑1 (MQ=255)                                                                                                                    cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCgg                            <  6:227272/149‑1 (MQ=255)                                                                                                                    cGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCgg                            <  2:188776/149‑1 (MQ=255)                                                                                                                     gAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCggg                           <  1:219774/149‑1 (MQ=255)                                                                                                                        tggctggcCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGa                        <  6:126610/149‑1 (MQ=255)                                                                                                                                ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACt                                 >  4:135340/1‑132 (MQ=255)                                                                                                                                ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACt                                 <  3:135340/132‑1 (MQ=255)                                                                                                                                ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGa                <  5:35953/149‑1 (MQ=255)                                                                                                                                ccccATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGa                <  3:279646/149‑1 (MQ=255)                                                                                                                                        atcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCaa        >  7:73995/1‑149 (MQ=255)                                                                                                                                        atcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCa         >  7:323514/1‑148 (MQ=255)                                                                                                                                         tcACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAAc       <  3:64253/149‑1 (MQ=255)                                                                                                                                            cTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGa                <  6:140300/137‑1 (MQ=255)                                                                                                                                            cTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGa                >  5:140300/1‑137 (MQ=255)                                                                                                                                               gagggCCCCCGTGGCGCACCCCAGCAAGCGGTTGAGGATGCAGATATACGGCGGTGATCAGCCTCAAACCACGCACTGTCCGCCAATCACGTTGAACGAAGCGAGTTCGCCCGAAATGGCGGGTGACCCGCTGAGCTTTCAACTTGcc  >  8:139953/3‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CTTGAGCCACGCCAACATCGTCACGGGCACGGACGGCCACTCGTACCTGTTGCGCTTCCACACCGCCGCGGCTCTGCTAGTGCTGGACAGGCACC‑GTGAGCTGCCTGGGGTTAGCGAATGGCTGGCCCCCATCCATCACTGGTGGGCCCCCGTGGCGCACCCGAGCAAGCGGTTGTGGATGCAGATATCCGGCGGTGATCAGCCTCAATCCACGCACTGTCCGCCAATCACGTTGAACGAAGCGTGTTGGGCCGAACTGGCGGGTGACCCGCTGAGCTTTCAACTTGCC  >  NC_002947/3835836‑3836124

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: ‑