Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,704	A→C	10.9%	L162R (CTG→CGG) ‡	glpF ←	aquaglyceroporin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,235,704	0	A	C	10.9%	151.7 / 11.1	55	L162R (CTG→CGG) ‡	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 = 8.52e-01

GCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAA  >  NC_002947/1235562‑1235846                                                                                                                                               |                                                                                                                                                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                                                                                   <  4:81028/116‑1 (MQ=255)                                                                                         gcgcCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCgg                                                                                   >  3:81028/1‑116 (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)                                                                                                                                               |                                                                                                                                                GCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAC‑GGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAA  >  NC_002947/1235562‑1235846

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