Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,712	A→C	10.9%	T159T (ACT→ACG)	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,712	0	A	C	10.9%	155.8 / 11.1	55	T159T (ACT→ACG)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base A (25/24);  new base C (3/3);  total (28/27)
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.82e-01

GTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTA  >  NC_002947/1235571‑1235858                                                                                                                                              |                                                                                                                                                    gTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATca                                                                                                                                              <  6:102530/149‑1 (MQ=255)      tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                                         >  6:203742/1‑149 (MQ=255)      tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                                         >  4:245581/1‑149 (MQ=255)         cAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGAt                                                                                                                                                >  1:119313/1‑139 (MQ=255)         cAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGAt                                                                                                                                                <  2:119313/139‑1 (MQ=255)          aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCaa                                                                                                                                     <  8:32787/149‑1 (MQ=255)            aGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAcc                                                                                                                                   <  5:124990/149‑1 (MQ=255)               cAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGa                                                                                                                                <  1:136648/149‑1 (MQ=255)                aGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGaa                                                                                                                               <  4:299905/149‑1 (MQ=255)                  tCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAAc                                                                                                                             >  7:182862/1‑149 (MQ=255)                         gCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggc                                                                                                                      >  5:37051/1‑149 (MQ=255)                          cccATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggct                                                                                                                     >  3:21473/1‑149 (MQ=255)                                gcgcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGAAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGa                                                                                                               <  5:218646/149‑1 (MQ=255)                                 cgcTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGAc                                                                                                              <  4:21473/149‑1 (MQ=255)                                      ccAATCACTGCAATCAGCAGGCCGTTCAGATGCTGGGCCATGTCGCCGCGAGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCTTCAGAATGGGAGTTATCACCACTTCAACGAGAAAGGCGTGTCTGGTGGACAGCg                                                                                                          <  8:18849/148‑1 (MQ=255)                                         aTCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑gg                                                                                                     >  2:38278/1‑149 (MQ=255)                                              tGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCg                                                                                                >  5:195513/1‑149 (MQ=255)                                                  aTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCtggt                                                                                                                  >  8:229281/1‑128 (MQ=255)                                                  aTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCtggt                                                                                                                  <  7:229281/128‑1 (MQ=255)                                                        agGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                                            >  5:177907/1‑94 (MQ=37)                                                        agGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                                            <  6:177907/96‑3 (MQ=37)                                                         ggCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGgag                                                                                     <  5:203742/149‑1 (MQ=255)                                                             gATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAaca                                                                                 >  7:13012/1‑149 (MQ=255)                                                                  gcagcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                                           <  1:17224/103‑1 (MQ=255)                                                                  gcagcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                                           >  2:17224/1‑103 (MQ=255)                                                                    agcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGc                                                                          <  8:13012/149‑1 (MQ=255)                                                                    agcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGc                                                                          <  8:182862/149‑1 (MQ=255)                                                                     gcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGc                                                                          >  1:12111/1‑148 (MQ=255)                                                                       ggggCCATGGCGCCGCGAGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAACCGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAg                                                                       >  4:19432/1‑149 (MQ=255)                                                                        gggCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACTGGGTGCGGGTAGGTGGAGAACACCGAGGCCAg                                                                       <  1:38278/149‑1 (MQ=255)                                                                           ccATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTc                                                                    >  4:109960/1‑148 (MQ=255)                                                                            cATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTc                                                                    >  3:175368/1‑147 (MQ=255)                                                                             aTGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCTGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACTCGTGGCTGGTGTACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAg                                                                 <  5:37628/149‑1 (MQ=255)                                                                             aTGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAg                                                                 >  4:256863/1‑149 (MQ=255)                                                                               ggCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGc                                                               >  5:167263/1‑149 (MQ=255)                                                                                gcgcCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCgg                                                                                               >  3:81028/1‑116 (MQ=255)                                                                                gcgcCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCgg                                                                                               <  4:81028/116‑1 (MQ=255)                                                                                  gccgcGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGc                                                            >  3:197108/1‑149 (MQ=255)                                                                                               ccGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACcgcg                                               >  3:196896/1‑149 (MQ=255)                                                                                                 gttgttGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGa                                                                                                         <  5:30179/90‑1 (MQ=255)                                                                                                 gttgttGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGa                                                                                                         >  6:30179/1‑90 (MQ=255)                                                                                                        tcgtcgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGc                                      <  6:37051/149‑1 (MQ=255)                                                                                                            cgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTgg                                  >  1:47717/1‑149 (MQ=255)                                                                                                                 aGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTg                             <  6:167263/149‑1 (MQ=255)                                                                                                                  gggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                                             >  4:24023/1‑132 (MQ=255)                                                                                                                  gggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                                             <  3:24023/132‑1 (MQ=255)                                                                                                                    gCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTc                          <  1:292218/149‑1 (MQ=255)                                                                                                                     ccATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCg                         <  3:245581/149‑1 (MQ=255)                                                                                                                      ggtGATCACTGCCATTCTGATGGCCGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCg                         >  4:89166/3‑148 (MQ=255)                                                                                                                               ggCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGaaa                      <  3:200667/142‑1 (MQ=255)                                                                                                                               ggCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACAGCGCAGCATGGCGTGGGCTTGTTCGaaa                      >  4:200667/1‑142 (MQ=255)                                                                                                                                   aTCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaagaa              >  3:157613/1‑146 (MQ=255)                                                                                                                                         aTGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGATAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCt     <  6:195513/149‑1 (MQ=255)                                                                                                                                         aTGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCt     <  4:175368/149‑1 (MQ=255)                                                                                                                                             cAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTa  >  6:206356/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                    GTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAACAGATTGCTGTA  >  NC_002947/1235571‑1235858

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