Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,235,391 G→A 8.0% P266P (CCC→CCT)  glpF ← aquaglyceroporin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,235,3910GA 8.0% 139.4 / 2.3 50P266P (CCC→CCT) glpFaquaglyceroporin
Reads supporting (aligned to +/- strand):  ref base G (24/22);  new base A (2/2);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

TCCAGGGCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGATTTCGCCCCAACCGGC  >  NC_002947/1235245‑1235531
                                                                                                                                                  |                                                                                                                                            
tCCAGGGCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCggcggg                                                                                                                                            <  8:42834/149‑1 (MQ=255)
  cAGGGCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCggcgggc                                                                                                                                           <  2:150515/148‑1 (MQ=255)
  cAGGGCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGCATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCggcgggc                                                                                                                                           <  6:156956/148‑1 (MQ=255)
     ggCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCa                                                                                                                                       >  8:199997/1‑149 (MQ=255)
       cGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCAtt                                                                                                                                     <  7:107463/149‑1 (MQ=255)
          tgatgTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGc                                                                                                                                  <  7:106273/149‑1 (MQ=255)
             tgTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAg                                                                                                                                <  8:111924/148‑1 (MQ=255)
                 gTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGc                                                                                                                           >  7:74596/1‑149 (MQ=255)
                                   tgtcATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGa                                                                                                          >  8:222632/1‑148 (MQ=255)
                                     tcATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAg                                                                                                       >  1:26563/1‑149 (MQ=255)
                                       aTGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCt                                                                                                     >  6:374549/1‑149 (MQ=255)
                                       aTGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCt                                                                                                     >  2:179405/1‑149 (MQ=255)
                                             gTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGc                                                                                               <  6:270879/149‑1 (MQ=255)
                                             gTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGc                                                                                               <  4:258520/149‑1 (MQ=255)
                                              tAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCa                                                                                              >  5:5580/1‑149 (MQ=255)
                                               aGGCCTTGCAAGAATAGGGGGGCAGGGCTGCGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCAc                                                                                             >  6:161019/1‑149 (MQ=255)
                                                ggCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCAcc                                                                                            >  3:143185/1‑149 (MQ=255)
                                                   cTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCAccc                                                                                           >  2:35026/1‑147 (MQ=255)
                                                     tGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCg                                                                                                   >  2:276447/1‑137 (MQ=255)
                                                     tGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCg                                                                                                   <  1:276447/137‑1 (MQ=255)
                                                      gCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGc                                                                                       <  4:143185/148‑1 (MQ=255)
                                                       cAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAg                                                                                     <  8:74596/149‑1 (MQ=255)
                                                        aaGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAgg                                                                                    >  2:91616/1‑149 (MQ=255)
                                                        aaGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAgg                                                                                    >  7:286994/1‑149 (MQ=255)
                                                        aaGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAgg                                                                                    >  2:185705/1‑149 (MQ=255)
                                                               ggggggCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCaa                                                                              <  6:5580/148‑1 (MQ=255)
                                                                        ggCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGc                                                                    <  2:96671/149‑1 (MQ=255)
                                                                              gCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATAc                                                              >  2:167432/1‑149 (MQ=255)
                                                                                 tGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCTCTGCACCCAAGCAGGCGCCAAGGATCGGTGc                                                                    <  7:150686/140‑1 (MQ=255)
                                                                                 tGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGACGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGc                                                                    >  8:150686/1‑140 (MQ=255)
                                                                                   aTGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGaa                                                         >  3:85873/1‑149 (MQ=255)
                                                                                   aTGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGAGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGaa                                                         <  5:374549/149‑1 (MQ=255)
                                                                                            tCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAAt                                                <  1:35026/149‑1 (MQ=255)
                                                                                               aTTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAgg                                             <  1:167432/149‑1 (MQ=255)
                                                                                               aTTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAgg                                             <  1:185705/149‑1 (MQ=255)
                                                                                                ttGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAggg                                            <  2:26563/149‑1 (MQ=255)
                                                                                                ttGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAggg                                            <  1:179405/149‑1 (MQ=255)
                                                                                                ttGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAggg                                            <  7:222632/149‑1 (MQ=255)
                                                                                                   gCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGa                                           >  5:367333/1‑147 (MQ=255)
                                                                                                   gCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGa                                           <  6:367333/147‑1 (MQ=255)
                                                                                                        aGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATccc                                    >  3:140319/1‑149 (MQ=255)
                                                                                                        aGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATcc                                     <  7:301834/148‑1 (MQ=255)
                                                                                                         gTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCAGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCTGTGCAAATACCGGAACCAGGAAATAGGGa                                           <  4:196480/141‑1 (MQ=255)
                                                                                                         gTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCAGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGa                                           >  7:82060/1‑141 (MQ=255)
                                                                                                         gTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCAGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGa                                           >  3:196480/1‑141 (MQ=255)
                                                                                                         gTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCAGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGa                                           <  8:82060/141‑1 (MQ=255)
                                                                                                                           tCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGa                  <  8:286994/148‑1 (MQ=255)
                                                                                                                             aTTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACGAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGAttt               >  8:142181/1‑149 (MQ=255)
                                                                                                                              tttGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGATTTc              >  3:90281/1‑149 (MQ=255)
                                                                                                                                    tcAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGATTTCGCCCCa        >  7:264705/1‑149 (MQ=255)
                                                                                                                                        gggTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGt                          >  2:96685/1‑127 (MQ=255)
                                                                                                                                        gggTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGt                          <  1:96685/127‑1 (MQ=255)
                                                                                                                                          gTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGATTTCGCCCCAACCGGc  >  1:163542/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
TCCAGGGCGATGATGTAGTTCTTATCCTGGGTGTCTGTCATGGTCGTAGGCCTTGCAAGAATAGGGGGGCAGGGCTGTGCGTGATGGTGGCATCAATTGGCTTGAGTATCGCCCTGGCGAATGTCATTTGTCTCAGGGGTCGCGGCGGGCGCCATTGGCAGGTTGCGTGCAATCAGGCCGCGATAGCTCGCTGCACCCAAGCAGGCGCCAAGGATCGGTGCAAATACCGGAACCAGGAAATAGGGAATATCCCGCCCGCCAGTGAAGGCGATTTCGCCCCAACCGGC  >  NC_002947/1235245‑1235531

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