Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,591,074 G→A L418L (CTG→TTG)  ugpB ← sn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,591,0740GA96.0% 65.4 / ‑6.1 25L418L (CTG→TTG) ugpBsn‑glycerol 3‑phosphate ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (8/16);  total (8/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCT  >  NC_000913/3590832‑3591324
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           
aTAACGGTGAGGATGAGCTGCGGCGCGAGCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGTGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCctgct                                                                                                                                                                                                                                                     <  1:148967/251‑1 (MQ=255)
   aCGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCtgtg                                                                                                                                                                                                                                                  <  1:347255/251‑1 (MQ=255)
             tgaGCGGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGATAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTCCACGCGATGCCAGGTCAGATAAGCGAAACGAATCCGACATTACACTGATTAAGACTTCGTCGCTTTCTCACAGCGGCGCAGCAAATGATTTCCCCGATCAACGGCGGTATCCAATGACTGCTGTGGTGTcttctt                                                                                                                                                                                                                                        >  2:304915/1‑251 (MQ=255)
              gaGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTcttctt                                                                                                                                                                                                                                        <  1:248363/250‑1 (MQ=255)
                gCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATAAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTAcc                                                                                                                                                                                                                                     <  1:90480/251‑1 (MQ=255)
                   gcggcgCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGt                                                                                                                                                                                                                                  <  1:111408/251‑1 (MQ=255)
                            ccAGCAGATAAGGAAGCCAGCGCGAGCGGAACACCAGACGGGATGATGAAATTCCGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCACAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCAACGCTCCACGGCGGATTCCAATGCCTGCTGGTGAGTCTTCTTACCgtttcacccgac                                                                                                                                                                                                                         >  2:317479/1‑240 (MQ=255)
                                          agccagcGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCAt                                                                                                                                                                                                           >  1:36711/1‑251 (MQ=255)
                                                            acCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCt                                                                                                                                                                                         >  1:148627/1‑251 (MQ=255)
                                                                                                         gATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAAcggcg                                                                                                                                           <  1:82687/251‑1 (MQ=255)
                                                                                                          aTGCGCTTCGCTTATCCGGCCTACACGCGGTGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATTCCTGCAGTGGTGTCTTCTTCCCGGTCCACAAGCTCTCCAGCTCTTCATCCAAAATCACGCGGATATACGGCATGTTGCCCACACTCACCCC‑TTTGGCTAACGGCAAcggcgg                                                                                                                                          >  2:114104/1‑251 (MQ=255)
                                                                                                                         ccGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTg                                                                                                                           <  1:317479/251‑1 (MQ=255)
                                                                                                                              cTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGc                                                                                                                                                                                       >  2:77178/1‑187 (MQ=255)
                                                                                                                              cTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGc                                                                                                                                                                                       <  1:77178/187‑1 (MQ=255)
                                                                                                                                    gcgATGCCAGGGCAGATAAGCGAAACGCCTCCGACATTACACTGATTAAGACTTCGGCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCgg                                                                                                                <  2:148627/251‑1 (MQ=255)
                                                                                                                                              ttCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCCACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCTGATCTTCATACCCAATCACGCCGAACTACGGAATGTTACCCAGAC‑CACCCCCTTTGGGGAACGGCAACAGCGGCTTATTcatcatatcacccgtcccggtttcgccccc                                                                                                      >  2:91980/2‑227 (MQ=255)
                                                                                                                                                             cGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGACCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGCCGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTAAGCATATGACGCGTAGCGGTATCCGCCCCTGGGGTTTTCTCAt                                                                                        >  1:21090/1‑250 (MQ=255)
                                                                                                                                                                          acacTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCa                                                                          <  1:233042/251‑1 (MQ=255)
                                                                                                                                                                           cacTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCATCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCa                                                                          <  1:337991/250‑1 (MQ=255)
                                                                                                                                                                                    aGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTTAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtcaggt                                                                <  1:28544/251‑1 (MQ=255)
                                                                                                                                                                                      aCTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGggtcaggtca                                                              <  1:127467/251‑1 (MQ=255)
                                                                                                                                                                                                     cAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGt                                                <  1:19576/250‑1 (MQ=255)
                                                                                                                                                                                                              cctgctacggATTTCCACGCTCACCGGTGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGa                                      <  1:172401/242‑1 (MQ=255)
                                                                                                                                                                                                                                  ccAACGGCGGTATCCAATGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCGCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCa                  <  1:357523/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cggcggTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCg              <  1:333089/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTTTGATGCCACTCGGCAGCGTTTTCt  >  2:79624/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                           
ATAACGGTGATGATGAGCTGCGGCGCGACCAGCAGATAAGGCAGCCAGCGCGAGCGGAACACCGGACGGGATGATGACATTACGGTTAATTCCTGAACTGTGCCGGATGCGCTTCGCTTATCCGGCCTACACGCGATGCCAGGTCAGATAAGCGAAACGCATCCGACATTACACTGATTAAGACTTCGTCGATTTCTCAAAGCGGCGCAGCAACTGATTTCCACGCTCAACGGCGGTATCCAGTGCCTGCTGTGGTGTCTTCTTACCGGTCCACACGCTCTCCAGCTCTTCATCCACAATCACGCGGATCTGCGGCATGTTGCCCAGACGCAGCCC‑TTTGGTGAACGGCAACGGCGGCTTATTCAGCATCTGACGCGTCGCGGTATCCGCCCCTGGGTTTTTCTCATAAAAGCCCTGCTCACGGGTCAGGTCATACGCTGCTTTGGTGATTGGCAGATAACCGGTTTTCTGATGCCACTCGGCAGCGTTTTCT  >  NC_000913/3590832‑3591324

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑