Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	748,089	G→T	10.4%	E48* (GAG→TAG)	PP_5449 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	748,089	0	G	T	10.4%	158.9 / 5.7	56	E48* (GAG→TAG)	PP_5449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (20/30);  new base T (4/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATGTTATGCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGTGATGAAGCGTAACGTTGATGAGCG  >  NC_002947/747941‑748231                                                                                                                                                     |                                                                                                                                                gaaTGTTATGCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGCCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGtgatg                                                                                                                                                 >  2:169736/1‑149 (MQ=255)       tATGCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCa                                                                                                                                            <  1:321562/148‑1 (MQ=255)       tATGCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCa                                                                                                                                            <  3:160817/148‑1 (MQ=255)          gCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCa                                                                                                                                            >  2:327320/1‑145 (MQ=255)             ttcGCTATTTAAATTTGCCTCAGTAGCTGTTGACGGTGACGTCATTAGTGGATGCGGTGGAGAGGATCTCACTGGGGCGTATAGATATAATGAGCAAATCTAAAAACGCGCAATGGTGATGAATATCCATGGAGATTAGGCAGAGAttt                                                                                                                                     >  2:332130/3‑149 (MQ=255)             tACGCTATTTAAATTTGCCTCAGTAGCTGTTGCAGGTGGCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATTGATATCTTGAGCAAATCTCAAAAGGCGCTATGGTGCTGAATATCCATGGTGATGAGGCAGAGAatt                                                                                                                                     <  7:291726/149‑1 (MQ=255)               cGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCTAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTc                                                                                                                                   <  1:88040/149‑1 (MQ=255)                 cTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTAATGATAGAATGCGCAAATCACAAAAGGCGCATGGTTGCTGAATATCCATGTTGAGGAGGCAGAGAGTTGCg                                                                                                                                  >  6:277343/1‑148 (MQ=255)                    tttAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCg                                                                                                                                  >  2:339992/1‑145 (MQ=255)                    tttAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGa                                                                                                                              <  4:183544/149‑1 (MQ=255)                     ttAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGa                                                                                                                              >  2:184956/1‑148 (MQ=255)                        aaTTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATAtt                                                                                                                          <  1:169736/149‑1 (MQ=255)                         aTTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTg                                                                                                                         >  1:290987/1‑149 (MQ=255)                            tGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTc                                                                                                                       <  4:79386/148‑1 (MQ=255)                                tCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATc                                                                                                                  >  6:125046/1‑149 (MQ=255)                                   gTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTgg                                                                                                               >  8:53464/1‑149 (MQ=255)                                    tAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGt                                                                                                              <  6:179878/149‑1 (MQ=255)                                      gCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGTAGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCAATGGTGATTATGCAGAGATTTTCGCAGATATTGTCGCATCTGGTAt                                                                                                            >  1:56490/1‑149 (MQ=255)                                      gCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTa                                                                                                             <  7:170887/148‑1 (MQ=255)                                      gCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGATTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTa                                                                                                             <  3:190283/148‑1 (MQ=255)                                               ggTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGt                                                                                                   <  5:125046/149‑1 (MQ=255)                                                 tGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTgg                                                                                                 <  2:46695/149‑1 (MQ=255)                                                      tCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAg                                                                                            <  1:327320/149‑1 (MQ=255)                                                      tCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAg                                                                                            >  3:186673/1‑149 (MQ=255)                                                            cTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGc                                                                                      <  1:184956/149‑1 (MQ=255)                                                                        gCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCaaaa                                                                           <  5:347832/148‑1 (MQ=255)                                                                        gCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCaaaa                                                                           <  1:339992/148‑1 (MQ=255)                                                                          gAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCTATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGt                                                                                                              <  1:77057/111‑1 (MQ=255)                                                                          gAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATTAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGa                                                                       >  4:185446/1‑149 (MQ=255)                                                                          gAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATTAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGa                                                                       >  4:179658/1‑149 (MQ=255)                                                                          gAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGt                                                                                                              >  2:77057/1‑111 (MQ=255)                                                                          gAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTCATGAGTCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCAGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGa                                                                       >  6:287257/1‑149 (MQ=255)                                                                          gAGGATTCCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGa                                                                       >  5:21163/1‑149 (MQ=255)                                                                             gATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCAGCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATgg                                                                    <  4:230550/149‑1 (MQ=255)                                                                              aTTTCACTGGGGCTTATACATATCATCAGCAAATCTCAAAACGCGCAATGGAGCTGAATATCCATGGGGATGAGGCAGAGATTTTCGCAGATTTTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGaggcg                                                                   <  8:166670/149‑5 (MQ=255)                                                                                          cTTATAGATATCATCACCAAATCACAAAAGCCCCAATGCTGCTTAATATCCATGGTCATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCtt                                                       <  8:228866/149‑1 (MQ=255)                                                                                             ataGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGAAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTtgat                                                    >  6:77384/1‑149 (MQ=255)                                                                                              taGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTtgat                                                    <  6:132335/148‑1 (MQ=255)                                                                                                       tGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATaaaaa                                            >  1:334386/1‑147 (MQ=255)                                                                                                         gttctaTCTCAAAAGGCGGAATGGTGCTGAATATGCATGGTGATGATGCAGAGATTTTCGCATATTTTGTCGCATCTTGTATTAATCCGGTGGGTAAGCTCAGCGTTTCGCAAAACGGATGGGAAGCTCATCCTTGATGATAAAAATAg                                         <  1:238986/144‑1 (MQ=255)                                                                                                               tctcAAAATGCGCAATGGTGCTGAATATCCATTGTGATGAGGCATAGATTTTCGCAGATATTGTCGCATCTTGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAACGGATGGGAAGCTCATCCTTGATGATAAAAATAg                                         <  1:233811/143‑1 (MQ=255)                                                                                                               tctcAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAg                                         >  2:233811/1‑142 (MQ=255)                                                                                                               tctcAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAg                                         >  2:238986/1‑142 (MQ=255)                                                                                                                ctcAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTAc                                 >  5:131940/1‑149 (MQ=255)                                                                                                                   aaacGGCGCAATGGTGCTGAATATCCATGGTTATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGAt                              <  6:37523/149‑1 (MQ=255)                                                                                                                   aaaaGGCGCAATGGTGCTGAATATCCATGGTGATTAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTa                                  <  3:185446/145‑1 (MQ=255)                                                                                                                   aaaaGGCGCAATGGTGCTGAATATCCATGGTGATTAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTa                                  <  3:179658/145‑1 (MQ=255)                                                                                                                   aaaaGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATaa                                               <  4:265852/132‑1 (MQ=255)                                                                                                                   aaaaGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATaa                                               >  3:265852/1‑132 (MQ=255)                                                                                                                     aaGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTg                            <  7:53464/149‑1 (MQ=255)                                                                                                                       ggCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCCCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGt                          <  2:56490/149‑1 (MQ=255)                                                                                                                          gcAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGtgat                       <  6:21163/149‑1 (MQ=255)                                                                                                                                  gCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGTGATGAAGCGt                >  1:157323/1‑148 (MQ=255)                                                                                                                                       atatCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCGCATCCTTGATGATAAAAATAGCTCGTTACGATTGGTGATGAAGCGTAACGtt          <  2:290987/149‑1 (MQ=255)                                                                                                                                         atCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATCGCTCGTTACGATTGGTGATGAAGCGTAACGTtga        >  8:117278/1‑149 (MQ=255)                                                                                                                                          tCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTc                                                                                 >  4:256441/1‑76 (MQ=255)                                                                                                                                          tCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTc                                                                                 <  3:256441/76‑1 (MQ=255)                                                                                                                                                 tgatgaGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGTGATGAAGCGTAACGTTGATGAGCg  >  1:191168/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                GAATGTTATGCGTACGCTATTTAAATTTGCCTCAGTAGCTGTTGCCGGTGCCGTCATTACTGGATGCGGTGGCGAGGATTTCACTGGGGCTTATAGATATCATGAGCAAATCTCAAAAGGCGCAATGGTGCTGAATATCCATGGTGATGAGGCAGAGATTTTCGCAGATATTGTCGCATCTGGTATTAAGCCGGTGGGTAAGCTCAGCGTTTCGCAAAA‑GGATGGGAAGCTCATCCTTGATGATAAAAATAGCTCGTTACGATTGGTGATGAAGCGTAACGTTGATGAGCG  >  NC_002947/747941‑748231

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