Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,731,689 C→T Q187* (CAA→TAA)  xylF → xylose ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,731,6890CT97.5% 116.6 / NA 40Q187* (CAA→TAA) xylFxylose ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (16/23);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATATTGTTCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACCAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCG  >  NC_000913/3731554‑3731826
                                                                                                                                       |                                                                                                                                         
gATATTGTTCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCAGCAGATAACAGCGCCAACCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGt                                                                                                                                    <  1:487607/143‑1 (MQ=255)
  tatTGTTCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGCAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGCTg                                                                                                                                  >  1:400116/1‑143 (MQ=255)
    ttGTTCCGCTAGGTAATTCCTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGCACAAATGAAAGTGTTAAAGCCTTACGTTGATTCCGTAAAAATTAAAGTCGTTGGTGACTAATGGGTTGAt                                                                                                                                <  2:234994/143‑1 (MQ=255)
      gTTCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATgg                                                                                                                              >  1:324464/1‑143 (MQ=255)
       ttCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACACATGAAAGTATTAAAAACCTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGAACAATGGATTGATGGc                                                                                                                             >  2:605134/1‑143 (MQ=255)
        tCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCt                                                                                                                            <  2:324464/143‑1 (MQ=255)
         ccGCGAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTg                                                                                                                           <  2:397109/143‑1 (MQ=255)
         ccGCAAGGTAATTACTTCCTGATTGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTg                                                                                                                           <  2:280747/143‑1 (MQ=255)
         ccGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTg                                                                                                                           <  1:546139/143‑1 (MQ=255)
          cGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGAGAAAACCTCGTGGTGTGACCTTCCGTCGATTCCGGAAAACTTAAAGTCGTTGGAGACTAATGGATGGATGGCTgg                                                                                                                          >  1:609650/1‑143 (MQ=255)
              aGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTAc                                                                                                                      <  1:157897/143‑1 (MQ=255)
                  aaTTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGa                                                                                                                  <  1:664520/143‑1 (MQ=255)
                         tCCTGATGGGTGGCTGGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCa                                                                                                           <  2:581431/143‑1 (MQ=255)
                          tcTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCAt                                                                                                          >  2:375706/2‑143 (MQ=255)
                          ccTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGGTTCCGGAAAAATTAAAGTCGTTGGTGACTAAAGGGTTGATGGCTGGTTACCGGAAAACGCAt                                                                                                          >  1:393621/1‑143 (MQ=255)
                               tGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGaaa                                                                                                     >  1:727881/1‑143 (MQ=255)
                               tGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTCGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGaaa                                                                                                     <  2:175641/143‑1 (MQ=255)
                                 ggcggcTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAAt                                                                                                   >  1:136707/1‑143 (MQ=255)
                                   cggcTCGCCTGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTa                                                                                                 <  1:404706/143‑1 (MQ=255)
                                     gcgcggcGGTAGAGAACGAGGCCGAGCTGGTGAGCGCCGGACAAATGAAAGTGTTAAGTCCTTACGGTGATTCCGTAAAAATTACAGTCGTTAAGGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATg                                                                                               <  2:292072/137‑1 (MQ=255)
                                     gcTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATg                                                                                               <  1:500729/143‑1 (MQ=255)
                                     gcTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATg                                                                                               <  1:498979/143‑1 (MQ=255)
                                                    aacacGCCAAGCTGTTCCGCGCCGGACGGCTCACAGCTGAAAAACCCTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAAcc                                                                                <  1:142296/140‑1 (MQ=255)
                                                      aacGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGCTTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAGATTATGGAAAACGCGCTAACCGc                                                                              >  2:165911/1‑143 (MQ=255)
                                                      aacGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGATGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGc                                                                              <  1:638060/143‑1 (MQ=255)
                                                                tGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAGAACGCATTGAAAATTCTGGAAAGCGCGCTAAGCGCACAGAATAAc                                                                    >  1:505275/1‑143 (MQ=255)
                                                                               aaaTGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGATGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTa                                                     >  1:96579/1‑143 (MQ=255)
                                                                                   gAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAACAATAACAAAATTGCTGCTGTAGATg                                                 >  2:552392/1‑143 (MQ=255)
                                                                                       gtgtTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTc                                             <  1:375706/143‑1 (MQ=255)
                                                                                          ttAAAACCTTACGCTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCaaa                                          <  1:144080/143‑1 (MQ=255)
                                                                                           tAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAAc                                         >  2:572711/1‑143 (MQ=255)
                                                                                              aaCCTTTCGTTTATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGAt                                      <  2:670137/143‑1 (MQ=255)
                                                                                              aaCCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGATTAATGGGTTGATGGATGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGAt                                      <  2:534555/143‑1 (MQ=255)
                                                                                               aCCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACGAAATTGATGCTGTAGTTGCCTCAAACGATg                                     >  1:459241/1‑143 (MQ=255)
                                                                                                ccTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACGGGAAAACGCATTGAAAATTATGGAAAACGCGCTAAGCGCCGAGAATGACAAAATTGATGCTGAAGTTGCCTCAAACGCTGc                                    >  2:87831/1‑143 (MQ=255)
                                                                                                              cGGAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGACAACGCGCTACCCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTgggg                      <  1:32016/143‑1 (MQ=255)
                                                                                                                gAAAAATTAAAGTCGTTGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAACCGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCCAAAGATGCCACCGCAGGGGGGGCa                    >  2:736376/1‑143 (MQ=255)
                                                                                                                               ttGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAgc     <  1:437495/143‑1 (MQ=255)
                                                                                                                                tGGTGACTAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAgcg    <  1:587957/143‑1 (MQ=255)
                                                                                                                                  gTGACTAAAGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCACTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAgcgcg  >  2:387010/1‑143 (MQ=255)
                                                                                                                                       |                                                                                                                                         
GATATTGTTCCGCAAGGTAATTACTTCCTGATGGGCGGCTCGCCGGTAGATAACAACGCCAAGCTGTTCCGCGCCGGACAAATGAAAGTGTTAAAACCTTACGTTGATTCCGGAAAAATTAAAGTCGTTGGTGACCAATGGGTTGATGGCTGGTTACCGGAAAACGCATTGAAAATTATGGAAAACGCGCTAACCGCCAATAATAACAAAATTGATGCTGTAGTTGCCTCAAACGATGCCACCGCAGGTGGGGCAATTCAGGCATTAAGCGCG  >  NC_000913/3731554‑3731826

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