Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,697	G→T	11.3%	A164A (GCC→GCA)	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,697	0	G	T	11.3%	143.3 / 10.0	53	A164A (GCC→GCA)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (23/24);  new base T (3/3);  total (26/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.71e-02

GCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAA  >  NC_002947/1235552‑1235838                                                                                                                                                  |                                                                                                                                               gCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGcc                                                                                                                                             <  2:53261/149‑1 (MQ=255)           gcgGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATgg                                                                                                                                   <  4:62642/149‑1 (MQ=255)             gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGaa                                                                                                                                      >  4:115924/1‑144 (MQ=255)             gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGaa                                                                                                                                      <  3:115924/144‑1 (MQ=255)                    gTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATca                                                                                                                          <  6:102530/149‑1 (MQ=255)                         tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                     >  4:245581/1‑149 (MQ=255)                         tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                     >  6:203742/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                                                                                                                        <  6:177907/96‑3 (MQ=37)                                                                           agGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                        >  5:177907/1‑94 (MQ=37)                                                                            ggCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGgag                                                                 <  5:203742/149‑1 (MQ=255)                                                                                gATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAaca                                                             >  7:13012/1‑149 (MQ=255)                                                                                     gcagcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                       >  2:17224/1‑103 (MQ=255)                                                                                     gcagcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                       <  1:17224/103‑1 (MQ=255)                                                                                       agcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGc                                                      <  8:182862/149‑1 (MQ=255)                                                                                       agcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGc                                                      <  8:13012/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                         <  3:24023/132‑1 (MQ=255)                                                                                                                                     gggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCgca                         >  4:24023/1‑132 (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)                                                                                                                                                  |                                                                                                                                               GCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAA  >  NC_002947/1235552‑1235838

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