Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,235,704	A→C	8.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	8.9%	167.6 / 7.5	56	L162R (CTG→CGG) ‡	glpF	aquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base A (24/27);  new base C (2/3);  total (26/30)
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.46e-01

GGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAA  >  NC_002947/1235564‑1235846                                                                                                                                             |                                                                                                                                               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                                                                                                                               <  8:100062/142‑3 (MQ=255)                 aaaGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACTGCCATTCTGATGGCCGTGATcatg                                                                                                                               >  7:100062/1‑140 (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                                                                                                                       <  2:305619/120‑1 (MQ=255)                                               aaTCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAAc                                                                                                                       >  1:305619/1‑120 (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)                                                                                                                                        ccATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaagaa  >  1:264734/1‑148 (MQ=255)                                                                                                                                         cATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaagaa  >  8:17688/1‑147 (MQ=255)                                                                                                                                          aTCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCgaagaa  >  6:138128/1‑146 (MQ=255)                                                                                                                                             |                                                                                                                                               GGGCCGGGTTCATCGCAAAGCCAGTCAATGGGCCCATGGCGCTGCCAATCACTGCAATCAGCAGGCCGATCAGCAGCGGGGCCATGGCGCCGCGCGGCAGGCCGTTGTTGTCGTCGGTCAGGGCCATGATCACGGCCATCAGAATGGCAGTGATCACCACTTCAACCAGAAACGCCTGGCTGGTGGACAGCGACGGGTGCGGGTAGGTGGAGAACACCGAGGCCAGCTCCAGGCTGCCTTCGCTACCGCGCAGCATGGCGTGGGCTTGTTCGAAATCGAAGAA  >  NC_002947/1235564‑1235846

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