Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,742 A→G 30.0% intergenic (‑107/+22) ygeA ← / ← araE Asp/Glu_racemase family protein/arabinose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,980,7420AG30.0% 58.2 / 31.4 40intergenic (‑107/+22)ygeA/araEAsp/Glu_racemase family protein/arabinose transporter
Reads supporting (aligned to +/- strand):  ref base A (14/14);  new base G (5/7);  total (19/21)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.54e-01

GGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGCAGTG  >  NC_000913/2980594‑2980890
                                                                                                                                                    |                                                                                                                                                    
ggATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTa                                                                                                                                                      <  1:718133/149‑1 (MQ=255)
 gATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     <  1:1009725/149‑1 (MQ=255)
        aGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     >  2:1238743/1‑142 (MQ=255)
        aGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     <  1:1238743/142‑1 (MQ=255)
        aGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACAt                                                                                                                                              <  3:513961/149‑1 (MQ=255)
          cTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATcc                                                                                                                                            <  3:1009577/149‑1 (MQ=255)
          cTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATcc                                                                                                                                            <  1:1710996/149‑1 (MQ=255)
          cTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATcc                                                                                                                                            <  1:1694366/149‑1 (MQ=255)
                         aaaCCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCg                                                                                                                                           <  3:1510481/135‑1 (MQ=255)
                         aaaCCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCg                                                                                                                                           >  4:1510481/1‑135 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     >  1:1673684/1‑99 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     >  1:1246675/1‑99 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     <  2:1137615/99‑1 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     >  1:1137615/1‑99 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     <  2:1246675/99‑1 (MQ=255)
                                                   cTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTAc                                                                                                                                                     <  2:1673684/99‑1 (MQ=255)
                                                       tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCgg                                                                                                                                          >  1:812912/1‑106 (MQ=255)
                                                       tgtgtgTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCgg                                                                                                                                          <  2:812912/106‑1 (MQ=255)
                                                                          gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAg                                                                                                                                                   <  1:1459533/78‑1 (MQ=255)
                                                                          gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAg                                                                                                                                                   >  2:1459533/1‑78 (MQ=255)
                                                                             gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGa                                                                                                                            <  1:1074181/98‑1 (MQ=255)
                                                                             gAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGa                                                                                                                            >  2:1074181/1‑98 (MQ=255)
                                                                                                                    cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                                                                       <  2:1421208/96‑1 (MQ=255)
                                                                                                                    cAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                                                                       >  1:1421208/1‑96 (MQ=255)
                                                                                                                                           gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                                                                     >  2:1669941/1‑91 (MQ=255)
                                                                                                                                           gCCGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAg                                                                     <  1:1669941/91‑1 (MQ=255)
                                                                                                                                            ccGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACaa                      >  2:1732973/1‑137 (MQ=255)
                                                                                                                                            ccGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACaa                      <  1:1732973/137‑1 (MQ=255)
                                                                                                                                            ccGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGcc                                        <  4:861089/119‑1 (MQ=255)
                                                                                                                                            ccGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGcc                                        >  3:861089/1‑119 (MQ=255)
                                                                                                                                              ggATGTGCTGTACATCCGGCCCGTGAGATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCcaca                       <  4:439327/134‑1 (MQ=255)
                                                                                                                                              ggATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCcaca                       >  3:439327/1‑134 (MQ=255)
                                                                                                                                              ggATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGa                                                                >  1:1315648/1‑93 (MQ=255)
                                                                                                                                              ggATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGa                                                                <  2:1315648/93‑1 (MQ=255)
                                                                                                                                               gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGaaa                                   <  3:1666914/121‑1 (MQ=255)
                                                                                                                                               gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGaaa                                   >  4:1666914/1‑121 (MQ=255)
                                                                                                                                               gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAgcg       >  2:179817/1‑149 (MQ=255)
                                                                                                                                                aTGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGaa                                              <  3:1610645/109‑1 (MQ=255)
                                                                                                                                                aTGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGaa                                              >  4:1610645/1‑109 (MQ=255)
                                                                                                                                                  gTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGCAg    >  3:1241255/1‑149 (MQ=255)
                                                                                                                                                   tACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTGCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGCAgt   >  4:166427/1‑149 (MQ=255)
                                                                                                                                                    aCAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGCAgtg  >  2:743472/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGATTCCCAGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGCAGTG  >  NC_000913/2980594‑2980890

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: