Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,697	G→T	13.2%	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	13.2%	143.8 / 11.8	53	A164A (GCC→GCA)	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base T (3/4);  total (24/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 = 8.21e-01

GGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAG  >  NC_002947/1235551‑1235844                                                                                                                                                   |                                                                                                                                                    ggCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACgg                                                                                                                                                    <  6:380646/148‑1 (MQ=255)  gCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGc                                                                                                                                                   >  6:239204/1‑148 (MQ=255)     cGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCAtt                                                                                                                                               >  6:173186/1‑148 (MQ=255)              gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGTTTTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCa                                                                                                                                      <  1:123364/149‑1 (MQ=255)              gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATgg                                                                                                                                        >  3:209313/1‑147 (MQ=255)              gggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCa                                                                                                                                      <  5:239204/149‑1 (MQ=255)               ggCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAg                                                                                                                                     <  5:210948/149‑1 (MQ=255)                      ttCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATcac                                                                                                                              <  2:344103/149‑1 (MQ=255)                        cATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                             <  5:173186/148‑3 (MQ=255)                          tCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACt                                                                                                                          >  5:74507/1‑149 (MQ=255)                             cAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTc                                                                                                                        >  3:195309/1‑148 (MQ=255)                              aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                             >  7:100062/1‑140 (MQ=255)                              aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                             <  8:100062/142‑3 (MQ=255)                               caGCCCGTCAATGGGCCCTTGGCGCTGCCAATCACTGCAATCAGCAGGCCGAGCAGCAGCGGGGCCATGGCGCCGCGCGGAAGTCCGTTGTTGTGGTCGGTCAGGGGCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCaag                                                                                                                     <  5:257740/148‑2 (MQ=255)                                  ccAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                     >  2:5358/1‑146 (MQ=255)                                        aaTGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                            <  6:235266/149‑1 (MQ=255)                                        aaTGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGCGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCAAGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGc                                                                                                            >  5:286875/1‑149 (MQ=255)                                         aTGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGcc                                                                                                           >  4:70373/1‑149 (MQ=255)                                          tGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCct                                                                                                          <  5:276114/149‑1 (MQ=255)                                             gCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggc                                                                                                       <  3:305159/149‑1 (MQ=255)                                               ccATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCctggctg                                                                                                     <  4:108832/149‑1 (MQ=255)                                                 aTGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCtggt                                                                                                   >  8:344093/1‑149 (MQ=255)                                                  tGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCtggtg                                                                                                  <  4:62690/149‑1 (MQ=255)                                                       cTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAg                                                                                             <  1:5358/149‑1 (MQ=255)                                                           cAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGAc                                                                                         <  4:209313/149‑1 (MQ=255)                                                            aaTCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                     >  1:305619/1‑120 (MQ=255)                                                            aaTCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                     <  2:305619/120‑1 (MQ=255)                                                                  tGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCAGTTGTTGTCGTCGGTCAGGGCCATCATCACGGCGATCAGAATGGCAGTGATCACGACTTCAACCAGAAACGCCTGGCTGGTGGACACCGACGGGTGCg                                                                                  >  6:194103/1‑149 (MQ=255)                                                                    cAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCggg                                                                                <  6:286875/149‑1 (MQ=255)                                                                     aaTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGC‑‑‑TGGACAGCGACGGGTGCGGGTAgg                                                                            <  3:70373/149‑1 (MQ=255)                                                                     aaTCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGACAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGt                                                                               <  6:24698/149‑1 (MQ=255)                                                                       tCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAg                                                                             >  5:91615/1‑149 (MQ=255)                                                                          gcagGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTg                                                                          <  4:195309/149‑1 (MQ=255)                                                                              gCCGATCAGCAGCGGGGCCATGGCGCCGAGCGGCAGGCCGTTGTTGTCGTCTGTCAGGGCCATGATCACGGCAATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGgag                                                                       >  8:270914/1‑148 (MQ=255)                                                                                cGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGAt                                                                                                                                 <  1:164080/88‑1 (MQ=37)                                                                                cGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGAt                                                                                                                                 >  2:164080/1‑88 (MQ=37)                                                                                cGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGaa                                                                     >  4:22459/1‑148 (MQ=255)                                                                                  aTCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAacac                                                                  >  6:68882/1‑149 (MQ=255)                                                                                      gcagcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAg                                                              >  7:71495/1‑149 (MQ=255)                                                                                        agcGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGc                                                            <  3:106262/149‑1 (MQ=255)                                                                                             ggCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCt                                                       >  6:345229/1‑149 (MQ=255)                                                                                                   ggCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGc                                                 >  4:122609/1‑149 (MQ=255)                                                                                                     cgccgcGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTg                                               <  6:91615/149‑1 (MQ=255)                                                                                                      gccgcGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAAGCCAGCTCCAGGCTGc                                              >  7:317966/1‑149 (MQ=255)                                                                                                                    cGTTGTTGTCGTCTGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCgcgc                                <  7:270914/149‑1 (MQ=255)                                                                                                                            tcgtcgGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGc                        <  3:22459/149‑1 (MQ=255)                                                                                                                                 ggTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTggg                   >  3:139211/1‑149 (MQ=255)                                                                                                                                  gTCTGGGCCATGATCACTGCCATCAGAATGGCAGTGATCACCACTTCAACCAGACACGCGTGGCTGGTGGACTGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCAAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGc                  <  5:194103/149‑1 (MQ=255)                                                                                                                                     aGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTg               <  3:122609/149‑1 (MQ=255)                                                                                                                                     aGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAAGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCtt                <  8:317966/148‑1 (MQ=255)                                                                                                                                       ggCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGtt             >  1:70741/1‑149 (MQ=255)                                                                                                                                        gCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTc            <  6:74507/149‑1 (MQ=255)                                                                                                                                         ccATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTc            >  1:301355/1‑148 (MQ=255)                                                                                                                                                 aCGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaa   <  5:68882/149‑1 (MQ=255)                                                                                                                                                   ggCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCGTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaag  >  6:205378/1‑148 (MQ=255)                                                                                                                                                   ggCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaag  <  5:205378/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    GGCCCGAAATCGCGGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAG  >  NC_002947/1235551‑1235844

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