Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	803,900	G→T	9.1%	S358I (AGC→ATC)	proB →	glutamate 5‑kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	803,900	0	G	T	9.1%	154.4 / 4.2	51	S358I (AGC→ATC)	proB	glutamate 5‑kinase
Reads supporting (aligned to +/- strand):  ref base G (20/26);  new base T (2/3);  total (22/29)
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

AGGGCAGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCGTGTTCAAGTTCGTCGGGC  >  NC_002947/803753‑804048                                                                                                                                                     |                                                                                                                                                     aGGGCAGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAgc                                                                                                                                                     >  7:277057/1‑149 (MQ=255)   ggCAGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAgcgc                                                                                                                                                   >  4:28621/1‑149 (MQ=255)    gCAGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGcc                                                                                                                                                  <  8:277057/149‑1 (MQ=255)      aGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGa                                                                                                                                                <  1:42804/149‑1 (MQ=255)         ttCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGcc                                                                                                                                             >  1:177572/1‑149 (MQ=255)                  ggcgATATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGGCTCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTg                                                                                                                                     <  1:250227/149‑1 (MQ=255)                   gcgAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTa                                                                                                                                   <  7:158875/149‑1 (MQ=255)                        aTGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACcg                                                                                                                              >  7:90346/1‑149 (MQ=255)                         tggtggTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACcgc                                                                                                                             <  3:28621/149‑1 (MQ=255)                         tggtggTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACcgc                                                                                                                             >  4:3855/1‑149 (MQ=255)                               tgtgCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAAc                                                                                                                       <  2:177572/149‑1 (MQ=255)                                 tgCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCt                                                                                                                     >  6:40630/1‑149 (MQ=255)                                  gCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTg                                                                                                                    <  3:199934/149‑1 (MQ=255)                                       ggCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCt                                                                                                               >  8:181750/1‑149 (MQ=255)                                        gCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTg                                                                                                              <  7:283776/149‑1 (MQ=255)                                          ccGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAgc                                                                                                                                                     >  8:242525/1‑108 (MQ=255)                                          ccGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAgc                                                                                                                                                     <  7:242525/108‑1 (MQ=255)                                              aCGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGa                                                                                                        <  3:43681/149‑1 (MQ=255)                                                ggTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGc                                                                                                      >  4:37778/1‑149 (MQ=255)                                                 gtgtTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGt                                                                                                                  >  8:240030/1‑136 (MQ=255)                                                 gtgtTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGt                                                                                                                  <  7:240030/136‑1 (MQ=255)                                                        aGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTga                                                                                               <  3:300874/148‑1 (MQ=255)                                                          tGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTgaaga                                                                                            <  6:246875/149‑1 (MQ=255)                                                            gCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAAc                                                                                          >  1:211483/1‑149 (MQ=255)                                                             cTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACg                                                                                         <  3:37778/149‑1 (MQ=255)                                                                   gCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCg                                                                                   <  3:3855/149‑1 (MQ=255)                                                                      tGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATATCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCtg                                                                                >  8:101677/1‑149 (MQ=255)                                                                      tGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATATCGCCGAGCCGGAACTTGTACACAGCGACAACCTGGTTCTGGTATGAGCCTGCTGAAGCACGTAATCGCtg                                                                                >  8:101679/1‑149 (MQ=255)                                                                      tGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCtg                                                                                >  2:321711/1‑149 (MQ=255)                                                                        gcgaaCTTCAGCGCGCTGTTGGCGGAGATGATCATTGG‑TCAGTCGTCCGAAGCCATCTAGTCGATCCTTGGCTTTAGCGCCGAGCCGGAACTTGTACAACGCGACAACCTTGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGgtgt                                                                               <  2:159651/145‑1 (MQ=255)                                                                         ccAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTgg                                                                             <  2:211483/149‑1 (MQ=255)                                                                           aaCTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCa                                                                           <  5:40630/149‑1 (MQ=255)                                                                                 aGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAgcgtgc                                                                     >  3:118468/1‑149 (MQ=255)                                                                                          gAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGcc                                                             <  8:90346/148‑1 (MQ=255)                                                                                          gAGGCGCAGAAGATCATGGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGcc                                                             >  6:318019/1‑148 (MQ=255)                                                                                                   aagaTCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTg                                                                                                              <  1:83466/90‑1 (MQ=255)                                                                                                   aagaTCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTg                                                                                                              >  2:83466/1‑90 (MQ=255)                                                                                                   aagaTCATTGG‑GCAGTCGTCCGCAGCCATCGAGTCGATCCTTGGCTATATCGCCGATCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAATAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTtggct                                                   <  7:101679/149‑1 (MQ=255)                                                                                                    agaTCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATATCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTtggct                                                   <  7:101677/148‑1 (MQ=255)                                                                                                     gaTCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTtggctgg                                                 >  6:120936/1‑149 (MQ=255)                                                                                                     gaTCAGTGG‑TCAGTCCTCCGAAGCCAACGAGTCGACCGTTGGCTATATCCCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTtggctgg                                                 <  8:144746/149‑1 (MQ=255)                                                                                                               tCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCACGCACGCTGccaca                                                          <  6:231561/131‑4 (MQ=255)                                                                                                               tCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCACGCACGCTGccaca                                                          >  5:231561/1‑128 (MQ=255)                                                                                                                 aGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAagag                                       <  7:181750/148‑1 (MQ=255)                                                                                                                  gtcgtcCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTCTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCCTGCTGCCGATGTTGTCTGGCGCCGAATAGGt                                     >  2:212557/1‑149 (MQ=255)                                                                                                                              cATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCtg                                                                                >  2:14315/1‑94 (MQ=255)                                                                                                                               aTCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTa                                                                                                           >  7:46588/1‑66 (MQ=255)                                                                                                                               aTCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTa                                                                                                           <  8:46588/66‑1 (MQ=255)                                                                                                                                tCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTtgg                                                     <  7:202561/119‑1 (MQ=255)                                                                                                                                   aGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCg                    <  1:321711/149‑1 (MQ=255)                                                                                                                                     tCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCgtg                  >  2:159998/1‑149 (MQ=255)                                                                                                                                       gATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCgtgt                 >  3:56989/1‑148 (MQ=255)                                                                                                                                           cTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGATGTCGGCCAGGTGTCTACCGTGTTCAAg            <  1:212557/149‑1 (MQ=255)                                                                                                                                                cTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCGTGTTCAAGTtcgt       <  4:118468/149‑1 (MQ=255)                                                                                                                                                  ataGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCGTGTTCAAGTtcgtcg     >  3:221744/1‑149 (MQ=255)                                                                                                                                                     gcgcCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCGTGTTCAAGTTCGTCGGGc  >  4:304609/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     AGGGCAGCTTCCGTCGCGGCGAGATGGTGGTGTGCGTCGGCCCGGACGGTGTTGAAGTGGCTCGCGGCCTGGCCAACTACAGCGCCCTGGAGGCGCAGAAGATCATTGG‑TCAGTCGTCCGAAGCCATCGAGTCGATCCTTGGCTATAGCGCCGAGCCGGAACTTGTACACCGCGACAACCTGGTGCTGGTATGAGCCTGCTGAAGAACGTAATCGCTGTGGCAGCGTGCGTGCTGCCGATGTTGGCTGGCGCCGAAGAGGTCGGCCAGGTGTCCACCGTGTTCAAGTTCGTCGGGC  >  NC_002947/803753‑804048

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