Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,746,128	G→A	20.3%	G33D (GGC→GAC)	PP_5042 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,746,128	0	G	A	20.3%	100.4 / 19.4	51	G33D (GGC→GAC)	PP_5042	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/21);  new base A (5/5);  total (25/26)
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.61e-01

AAGGTTGTTCATAAATTGAACAATTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGCATCATCGAGCGCTTCAAC  >  NC_002947/5745980‑5746276                                                                                                                                                     |                                                                                                                                                      aaGGTTGTTCATAAATTGAACAATTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATgg                                                                                                                                                       >  1:133316/1‑149 (MQ=255)           aTAAATTGAACAATTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTc                                                                                                                                              >  4:206104/1‑148 (MQ=255)            tAAATTGAACAATTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGa                                                                                                                                            >  7:292541/1‑149 (MQ=255)                      aaTTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTg                                                                                                                                                    >  3:196991/1‑131 (MQ=255)                        ttCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCCCCTGGAGACATGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCAccc                                                                                                                                >  7:156366/1‑149 (MQ=255)                              tATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCg                                                                                                                                             >  4:16573/1‑130 (MQ=255)                              tATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCg                                                                                                                                             <  3:16573/130‑1 (MQ=255)                                  atcGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTgg                                                                                                                       <  6:117075/148‑1 (MQ=255)                                  atcGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACCTGCTGCACCCTGCCGGTgg                                                                                                                       <  6:148099/148‑1 (MQ=255)                                     gcgcgcCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAgg                                                                                                                    <  6:227788/148‑1 (MQ=255)                                       gcgcCCCAAGCCCGTGATCCACCGTATAAATATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGa                                                                                                                 <  8:156366/149‑1 (MQ=255)                                           cccAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTg                                                                                                                  >  2:109534/1‑144 (MQ=255)                                           cccAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGa                                                                                                              >  4:184152/1‑148 (MQ=255)                                            ccAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGa                                                                                                                      <  6:349585/139‑1 (MQ=255)                                            ccAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGa                                                                                                                      >  5:349585/1‑139 (MQ=255)                                             cAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAg                                                                                                                     <  8:209437/139‑1 (MQ=255)                                             cAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAg                                                                                                                     >  7:209437/1‑139 (MQ=255)                                               aGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAg                                                                                                         <  2:332159/149‑1 (MQ=255)                                               aGACCGGGATCCACCCTATAATGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACCTGTTGCACCCTGCCGGTGGAGGTGAAGATCAAg                                                                                                         <  8:292541/149‑1 (MQ=255)                                                   cGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATc                                                                                                    <  3:206104/149‑1 (MQ=255)                                                   cGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATc                                                                                                    <  2:133316/149‑1 (MQ=255)                                                   cGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATc                                                                                                    >  6:359240/1‑149 (MQ=255)                                                             ccTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTAt                                                                                          >  6:225189/1‑149 (MQ=255)                                                                 taAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGc                                                                                      >  5:233961/1‑149 (MQ=255)                                                                 taAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGc                                                                                      >  6:25446/1‑149 (MQ=255)                                                                 taAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCGGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGc                                                                                      <  5:25446/149‑1 (MQ=255)                                                                  aaaGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGaa                                                                                                                >  2:337089/1‑123 (MQ=255)                                                                  aaaGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGaa                                                                                                                <  1:337089/123‑1 (MQ=255)                                                                           aaCGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACctgc                                                                                                                                       >  8:251319/1‑91 (MQ=255)                                                                           aaCGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACctgc                                                                                                                                       <  7:251319/91‑1 (MQ=255)                                                                              gACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTg                                                                                                                                                    >  3:167266/1‑75 (MQ=255)                                                                              gACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTg                                                                                                                                                    <  4:167266/75‑1 (MQ=255)                                                                                   tGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATc                                                                                                    <  4:304407/117‑1 (MQ=255)                                                                                   tGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATc                                                                                                    >  3:304407/1‑117 (MQ=255)                                                                                        gccgcGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGa                                                               <  5:225189/149‑1 (MQ=255)                                                                                          cgcgGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAgg                                                                                                                    >  6:260994/1‑95 (MQ=255)                                                                                          cgcgGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAgg                                                                                                                    <  5:260994/95‑1 (MQ=255)                                                                                                cGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGa                                                               <  1:36254/141‑1 (MQ=255)                                                                                                cGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGa                                                               >  2:36254/1‑141 (MQ=255)                                                                                                       gACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGa                                                 >  6:13617/1‑148 (MQ=255)                                                                                                           gCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTg                                            <  1:109534/149‑1 (MQ=255)                                                                                                            ccTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTg                                                                                                  <  6:7054/94‑1 (MQ=255)                                                                                                             cTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTTAAGATCAAGGCCTCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTgg                                           <  1:102884/149‑1 (MQ=255)                                                                                                                     cTGGCAGCGCTACACCAGCAACATGTGCCATGTCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGcc                                  >  8:306341/1‑149 (MQ=255)                                                                                                                        gCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTg                                >  3:108407/1‑148 (MQ=255)                                                                                                                        gCAGCGCTACACCAGCAACATGTGCCATGACTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTg                                >  6:133181/1‑148 (MQ=255)                                                                                                                                 cacCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAaagaa                       <  3:184152/148‑1 (MQ=255)                                                                                                                                 cacCAGCAACATGTGCCATGACTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAaagaag                      <  5:359240/149‑1 (MQ=255)                                                                                                                                     aGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGca                   <  2:86327/148‑1 (MQ=255)                                                                                                                                     aGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGca                   >  1:86327/1‑148 (MQ=255)                                                                                                                                              gtgtagcgCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGCATCATCGAgc          <  3:347838/141‑1 (MQ=255)                                                                                                                                                  aTGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAg                                      <  6:312671/116‑1 (MQ=255)                                                                                                                                                  aTGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAg                                      >  5:312671/1‑116 (MQ=255)                                                                                                                                                     gCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATAAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGCATAATCGAGCGCTTCAAc  >  6:33389/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AAGGTTGTTCATAAATTGAACAATTCCGGTATCATCGCGCGCCCCAAGACCGTGATCCACCCTATAAAGATGAAAACGACCCTGATCGCCGCGGCCGAAGTCGACCGCCTGGAGACCTGGCAGCGCTACACCAGCAACATGTGCCATGGCTGCCATTCGACCTGCTGCACCCTGCCGGTGGAGGTGAAGATCAAGG‑ATCTGATCCGTATCGGCGTGGTCGACGAGTTCGAAAAAGACGAACCACCCAAGAACGTGGCCAAGCGCCTGCAAAAAGAAGGCATCATCGAGCGCTTCAAC  >  NC_002947/5745980‑5746276

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