Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,238,644	G→T	8.0%	G52G (GGG→GGT)	PP_5701 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,238,644	0	G	T	8.0%	142.3 / 2.0	46	G52G (GGG→GGT)	PP_5701	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/21);  new base T (2/2);  total (23/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 = 1.00e+00

GCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCC  >  NC_002947/5238496‑5238786                                                                                                                                                     |                                                                                                                                                 gCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTTCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAggg                                                                                                                                                  <  6:269423/149‑1 (MQ=255) gCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAggg                                                                                                                                                  <  8:151258/149‑1 (MQ=255)   gACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGggcc                                                                                                                                                >  8:47546/1‑149 (MQ=255)    aCCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGggccg                                                                                                                                               <  8:93806/149‑1 (MQ=255)                aaGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCctcttcgcgg                                                                                                                                   >  7:47546‑M2/1‑139 (MQ=255)                     aCCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCg                                                                                                                              >  1:180163/1‑149 (MQ=255)                              cATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTccc                                                                                                                     <  3:277714/149‑1 (MQ=255)                              aaTGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGTCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTTCCGGAAGTGCCGGCCTCTTCGCGGGCTCGCCCGCTccc                                                                                                                     <  3:272951/148‑1 (MQ=255)                                tGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCca                                                                                                                    >  3:69519/1‑148 (MQ=255)                                     cGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCgg                                                                                                              <  6:101104/149‑1 (MQ=255)                                     cGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCgg                                                                                                              <  8:228623/149‑1 (MQ=255)                                        tGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGCAGGTCCGGCcctttcgcgggcacgcccgctcccacaggtactg                                                                                                           <  2:58292‑M2/149‑35 (MQ=255)                                         gCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAg                                                                                                                 <  7:300303/142‑1 (MQ=255)                                         gCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAg                                                                                                                 >  8:300303/1‑142 (MQ=255)                                          cATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTcac                                                                                                          <  2:180163/148‑1 (MQ=255)                                           aTGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCctcttcgcg                                                                                                                                    <  8:119369‑M2/121‑10 (MQ=255)                                           aTGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCctcttcgcg                                                                                                                                    >  7:119369‑M2/1‑112 (MQ=255)                                            tGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCca                                                                                                                    <  7:314837/136‑1 (MQ=255)                                            tGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCca                                                                                                                    >  8:314837/1‑136 (MQ=255)                                                cTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCcac                                                                                                                   >  3:216299/1‑133 (MQ=255)                                                cTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCcac                                                                                                                   <  4:216299/133‑1 (MQ=255)                                                cTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑agta                                                                                                  <  4:69519/149‑1 (MQ=255)                                                 tGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               <  1:179705/136‑1 (MQ=255)                                                               ccTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGca                                                                                    >  2:24053/1‑147 (MQ=255)                                                                   ttcttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               >  3:280537/1‑118 (MQ=255)                                                                   ttcttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCg                                                                                                               <  4:280537/118‑1 (MQ=255)                                                                   ttcttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGt                                                                                <  1:88514/147‑1 (MQ=255)                                                                     cttGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCt                                                                            >  8:197175/1‑149 (MQ=255)                                                                         aGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTgg                                                                        <  6:172870/149‑1 (MQ=255)                                                                          gTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGTCCGGCcctttcgcgggcacgcccgctcccacaggtactgtgaaaggcttgactagtgtggtgaccgctgtggg                                                                         <  2:327916‑M2/149‑69 (MQ=255)                                                                             ggCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGCAGGTCCGGCcctttcgcgggcacgcccgctcccacaggtactcgtgcaagtcttgactactgtggtgaccgctgtggga                                                                       >  1:58292‑M2/1‑78 (MQ=255)                                                                             ggCGCTGGGTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGTCCGGCcctttcgcgggcacgcccgctcccacaggtactgtgcaagtcttgactactgtggtgaccgctgtgggag                                                                       >  1:327916‑M2/1‑78 (MQ=255)                                                                                         ggTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCcgc                                                        <  7:197175/149‑1 (MQ=255)                                                                                               gATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCctcttcgcggg                                                                                                                                  >  2:323516‑M2/1‑60 (MQ=255)                                                                                               gATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGaaaag                                                  >  7:280093/1‑147 (MQ=255)                                                                                                 tGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGCCTTGCACAGTACCAGTGGGAGCGGGCGTGCCCGCGAAAGGGc                                                >  3:9245/1‑149 (MQ=255)                                                                                                 tGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGc                                                <  1:24053/149‑1 (MQ=255)                                                                                                   tCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGTCCg                                              >  4:236479/1‑149 (MQ=255)                                                                                                   tCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAAGGCCg                                              <  8:280093/149‑1 (MQ=255)                                                                                                    ccGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGggccggc                                                                                                                                             <  3:59049/55‑1 (MQ=255)                                                                                                    ccGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGggccggc                                                                                                                                             >  4:59049/1‑55 (MQ=255)                                                                                                     cGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGa                                            >  8:92849/1‑149 (MQ=255)                                                                                                                ggAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACa                                  >  8:18244/1‑148 (MQ=255)                                                                                                                  aaCTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAat                               <  7:92849/149‑1 (MQ=255)                                                                                                                     ttCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATc                            >  4:31267/1‑149 (MQ=255)                                                                                                                            tAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGca                                                                                    <  4:280803/86‑1 (MQ=255)                                                                                                                             aaGGGTGGTTACCGGTCCGGATGGGCCGGCCTCTTCGGGGGCTCGCCGGCTCCCACACCGGTCACCAC‑TTTAGT‑CAATAGTTTCAAATTCCCTGTGGGAGGGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGcc                     <  4:9245/148‑1 (MQ=255)                                                                                                                                    gTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCACAAGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATAACCCCGGCCATGCAcc              >  1:104585/1‑149 (MQ=255)                                                                                                                                          gTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGc                                                                                     >  6:87366/1‑71 (MQ=255)                                                                                                                                          gTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGc                                                                                     <  5:87366/71‑1 (MQ=255)                                                                                                                                          gTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAg       >  6:80444/1‑149 (MQ=255)                                                                                                                                                gAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGTGCAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAgccgcc  >  1:270664/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                 GCGACCGATTGCCGCAAGGCACCTTGAGGCATGGCACGATGCATGAACTGGAGGCGCTGCGGCCTCTTCTTGAGTCGGCGCTGGTTCAGGTGCAGATGTCCGCTGAGCCCAGGAACTTCGCACTAAGGCTGGTTAACGTGCCGGAAGGGCCGGCCTCTTCGCGGGCTCGCCCGCTCCCACAGCGGTCACCAC‑AGTAGT‑CAAGACTTGCACAGTACCTGTGGGAGCGGGCGTGCCCGCGAAAGGGCCGGACCTGCCGACAAATATCCCCGGCCATGCACCGTCGCAGCCGCC  >  NC_002947/5238496‑5238786

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