Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,690,616	C→G	7.5%	I12M (ATC→ATG)	gshB →	glutathione synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,690,616	0	C	G	7.5%	156.8 / 2.0	50	I12M (ATC→ATG)	gshB	glutathione synthetase
Reads supporting (aligned to +/- strand):  ref base C (16/29);  new base G (2/2);  total (19/31)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

GGGTGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGCTGAAGGTGT  >  NC_002947/5690472‑5690761                                                                                                                                                 |                                                                                                                                                  gggTGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGt                                                                                                                                               <  1:288968/149‑1 (MQ=255)   gTGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTc                                                                                                                                              >  4:195141/1‑148 (MQ=255)   gTGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCAAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTcc                                                                                                                                             <  3:212315/149‑1 (MQ=255)    tGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCa                                                                                                                                            <  5:209592/149‑1 (MQ=255)        cATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCAtctc                                                                                                                                        <  5:23461/149‑1 (MQ=255)          ttCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTcc                                                                                                                                       >  7:290329/1‑148 (MQ=255)              aTCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACaa                                                                                                                                  <  2:110659/149‑1 (MQ=255)                   ccGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACaagaag                                                                                                                              <  5:202326/148‑1 (MQ=255)                               aCCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGGAGGACAGCTCGCTgg                                                                                                                 >  7:213395/1‑149 (MQ=255)                                 cGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGcc                                                                                                               <  4:285431/149‑1 (MQ=255)                                   tCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCAt                                                                                                             >  1:274213/1‑149 (MQ=255)                                      gCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCAtgct                                                                                                          >  3:149594/1‑149 (MQ=255)                                                   cACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGc                                                                                                                        <  5:159945/122‑1 (MQ=255)                                                   cACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGc                                                                                                                        >  6:159945/1‑122 (MQ=255)                                                      ccTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCa                                                                                          <  3:195141/149‑1 (MQ=255)                                                       cTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCAc                                                                                         <  6:279326/149‑1 (MQ=255)                                                            ccACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGc                                                                                    <  3:131951/149‑1 (MQ=255)                                                            ccACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCATCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGc                                                                                    <  2:261324/149‑1 (MQ=255)                                                             cACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCt                                                                                   <  1:285974/149‑1 (MQ=255)                                                              aatattctttcGATCCACCGCGGCCACCTGACGCCAAGAGGACCACCCATGAGCGTTAGCCTCGGAATTGTGATCGACCCCATAGCGTCAATCTCCTACAAGAAGGACAGCTCACTGTCCATGCTGCAGGACGCCAAGGCACGCGCCTg                                                                                  >  3:273250/5‑149 (MQ=255)                                                              aCTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTg                                                                                  <  2:281702/149‑1 (MQ=255)                                                                  ttctttcGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGa                                                                                <  2:311227/147‑1 (MQ=255)                                                                  ttctttcGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGa                                                                                <  4:68192/147‑1 (MQ=255)                                                                         gATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTc                                                                        <  4:149594/148‑1 (MQ=255)                                                                          aTCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTa                                                                      >  1:17144/1‑149 (MQ=255)                                                                          aTCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTa                                                                      >  8:293669/1‑149 (MQ=255)                                                                              aCCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTccatg                                                                    >  4:234032/1‑143 (MQ=255)                                                                              aCCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATg                                                                  <  2:17144/149‑1 (MQ=255)                                                                                   ggCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACa                                                             <  7:298488/149‑1 (MQ=255)                                                                                     ccACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCCATGTAGAACAgg                                                           <  3:234032/149‑2 (MQ=255)                                                                                      cACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAACGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGc                                                           >  5:207043/1‑148 (MQ=255)                                                                                       aCCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCaa                                                         <  4:270203/149‑1 (MQ=255)                                                                                       aCCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCaa                                                         <  1:50571/149‑1 (MQ=255)                                                                                       aCCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCAACGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCaa                                                         <  3:181428/149‑1 (MQ=255)                                                                                            ccgccAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGGAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGAc                                                      <  8:213395/147‑1 (MQ=255)                                                                                               ccAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGActcctg                                                 >  8:31321/1‑144 (MQ=255)                                                                                                 aaGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATc                                               <  5:127473/149‑1 (MQ=255)                                                                                                    agGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAgg                                            <  8:290329/149‑1 (MQ=255)                                                                                                        ccaccCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCg                                         >  5:118116/1‑148 (MQ=255)                                                                                                          accCATGAGCGTGCGCCTCGGCATTGTGATGGCCCCCATGGCGTCCATCTCATACAAGAACGACAGCTCGCTGGGCATGCTGCTTGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGaag                                      <  4:63620/149‑1 (MQ=255)                                                                                                           cccATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAg                                                        >  4:126175/1‑130 (MQ=255)                                                                                                           cccATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAg                                                        <  3:126175/130‑1 (MQ=255)                                                                                                            ccATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGaagg                                     >  3:168771/1‑148 (MQ=255)                                                                                                                                   gagATGGACGCGATGGGGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAgg                                            >  2:199609/3‑118 (MQ=255)                                                                                                                                   gagATGGACGCGATGGGGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAgg                                            <  1:199609/116‑1 (MQ=255)                                                                                                                                     gATGGACGCGATGGGGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGc           >  2:293974/1‑149 (MQ=255)                                                                                                                                      aTGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGc           <  2:274213/148‑1 (MQ=255)                                                                                                                                          aCCCCATCGCGTCCATCTCCTACAAGCAGGACAGCGCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGCTGAAg      <  7:242264/149‑1 (MQ=255)                                                                                                                                            cccATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGCTGAAGgt    >  7:189982/1‑149 (MQ=255)                                                                                                                                              cATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGCTGAAGgtgt  >  2:241750/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GGGTGGCCATTCCATCTGCCGCGCCCCTTGACCGTCGGCGTTCGCGGCGCCACCCTGAGCCACTTTTCTTTCGATCCACCGCGGCCACCTGCCGCCAAGAGGACCACCCATGAGCGTTCGCCTCGGCATTGTGATGGACCCCATCGCGTCCATCTCCTACAAGAAGGACAGCTCGCTGGCCATGCTGCTGGCCGCCCAGGCACGCGGCTGGAGCCTGTTCTACATGGAACAGCAAGACCTGTATCAGGGCGAAGGCAAGGCCCGCGCCCGCATGCGCCCGCTGAAGGTGT  >  NC_002947/5690472‑5690761

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