Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,980,744 A→T 10.0% intergenic (‑109/+20) ygeA ← / ← araE amino acid racemase/arabinose:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,980,7440AT10.0% 93.1 / 2.4 38intergenic (‑109/+20)ygeA/araEamino acid racemase/arabinose:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (21/13);  new base T (2/2);  total (23/15)
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.62e-01

AGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGC  >  NC_000913/2980602‑2980886
                                                                                                                                              |                                                                                                                                              
aGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACa                                                                                                                                                <  1:68382/143‑1 (MQ=255)
         ctccCAGCAACCCAGTTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTTAATGACGAACGAGACAATCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTGCTGTACATCCg                                                                                                                                       <  2:317010/143‑1 (MQ=255)
          tccCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGCCAAACGCCTCAACGGCCAAGTGCCCAATCACTATTAACGAAAAAAGGGCCGGATGTACAGGACATCCgg                                                                                                                                      >  2:915737/1‑143 (MQ=255)
           cccAGCAAACCAATTGTTTTCATTCCTATCTCCTGATGTGTGTTCCCGCAACAAGGGTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAAGGCCGGAGGGAGAGCACAGCCgga                                                                                                                                     >  2:712279/1‑142 (MQ=255)
              aGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCg                                                                                                                                  >  1:342907/1‑143 (MQ=255)
                   aCCATTTGTTTTCCTGCCTCTCGCCTGATGTGTGTCACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGaaa                                                                                                                             <  2:736480/143‑1 (MQ=255)
                           gttggaTTCCTCTCGCCTGATGTGCGTTACCGCAACAAGGCGGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCACAATAAACAAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGc                                                                                                                     <  2:168211/138‑1 (MQ=255)
                                 ttCCTCTCTCCTGAGGTGTGTTACCGCAACAAGGCGGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGatat                                                                                                               <  1:535085/143‑1 (MQ=255)
                                        ctcCAGACGTGTGCTACCGCAACAAGGCTGAATAAGGAACGAGGCAAACGCCCCAACGGCCAAGTGCCCAATCACAATAAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCaa                                                                                                        <  2:222902/143‑1 (MQ=255)
                                                         cGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCCGCCATcc                                                                                       >  1:838221/1‑142 (MQ=255)
                                                          gCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGCAAAAAGGGCCGGATATACAGCACATCCGCCCCGTGAAATCAGACGCCGACATTTCTCAACTTCTCGCCTGCAACCAg                                                                                      >  2:161914/1‑143 (MQ=255)
                                                             acaaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGttt                                                                                   <  1:856135/143‑1 (MQ=255)
                                                               aaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCCTCAGTTTGc                                                                                 >  1:297886/1‑143 (MQ=255)
                                                               aaGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCCCAACTTCTCGCCGGGCATCCGTTTGc                                                                                 >  2:589020/1‑143 (MQ=255)
                                                                  gCTGAATAAGGAACGAGACAAACGCCTCAACGTCCAAGTGCCCAATCTTTATTAACGAAAAAAGTGCCGGATGTACAGCACAGCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGtt                                                                              <  1:273428/143‑1 (MQ=255)
                                                                  gCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGtg                                                                              >  1:2411/1‑142 (MQ=255)
                                                                           ggAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGtt                                                                     >  1:533733/1‑143 (MQ=255)
                                                                            gAACGAGACAAACGCCTAAACGGCCAAGAGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGAGAAATCAGACGCCGATATTTCACAACGTCGCGCCTGCCATCAGTTTGAGATCGATATGTTc                                                                    >  2:724970/1‑143 (MQ=255)
                                                                               cGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGACATGTTCAAg                                                                 >  1:81936/1‑143 (MQ=255)
                                                                                  gaCAAACGCCTCAACGGCCAAGTGACCAATCGCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAAACAGACGCCGATATTTCTCAACTGCTCGCCTGCCAAAAGTTTGCGTTCCAGATGttaaacaaa                                                              >  2:94678/1‑136 (MQ=255)
                                                                                     aaaCGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGAc                                                           >  1:11660/1‑143 (MQ=255)
                                                                                                     aggTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGa                                           <  2:760489/141‑1 (MQ=255)
                                                                                                       gTGCCCACTCTCTATTAGCGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAAt                                         <  1:602660/143‑1 (MQ=255)
                                                                                                       gTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCGCGTGAAATCAGCCGACGGTATTTGTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAAt                                         <  2:550069/143‑1 (MQ=255)
                                                                                                         gctcAATCTCTATTAACGAAAAAAGGTTCGGATGTACTTCACATCCGGCCCGTGAAATCAGACGCCGATATTGCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGa                                       <  1:223735/140‑1 (MQ=255)
                                                                                                          cccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGCACTTTTTAGTTTACGGAATGAg                                      >  2:371540/1‑143 (MQ=255)
                                                                                                           ccAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGc                                     <  1:517091/143‑1 (MQ=255)
                                                                                                                    aTTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGGTTCCGGAATGAGCCAGAAAGTa                            >  1:43642/1‑143 (MQ=255)
                                                                                                                     ttAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGTCCTTTTTGGTTTCCGTAATCAGCCTGAATGTag                           >  1:229148/1‑142 (MQ=255)
                                                                                                                        gcgctgcaaGTGCGGGGTGTACAGCACATCCGGCCCGTGAGCTCAGACGCGGATGTTTCTCAACTTCTCGCCTGCCATCCGTTTGCGTTCGATAGGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGc                        <  2:573697/136‑1 (MQ=255)
                                                                                                                          gAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCACAAACCAATGccc                      >  2:353817/1‑143 (MQ=255)
                                                                                                                           aaaaaaGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCACCTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAGTGAGCGAGCAGGTACTGCTca                     >  1:182474/1‑143 (MQ=255)
                                                                                                                                    cccgATGTCCACCCCCTCCGGCCCGTGAAACCAGCCGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCaa            <  1:84021/140‑1 (MQ=255)
                                                                                                                                    ccGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAAGGCCCACAAACGcga            >  2:518047/1‑141 (MQ=255)
                                                                                                                                    ccGGATGTGCTGTACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTGTGCGTTCGATATGTTCCAGCGTGACATTTTTGGCTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGGaa            >  1:383482/1‑143 (MQ=255)
                                                                                                                                       gATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGCATGAGCCAGAAAGTAATGCCCCAAAACACACTGt         >  2:883108/1‑143 (MQ=255)
                                                                                                                                        aTGTACATCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGtt        <  1:389857/143‑1 (MQ=255)
                                                                                                                                        aTGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTGCGGAATGAGCCAGCCAGTACTGCCCACAAACGGCATGtt        >  1:56394/1‑143 (MQ=255)
                                                                                                                                            aCAGCACATCCGGCCCGTGAAATCAGAAGCCGATAATTATCAACTTCTCGCCCGCCATCAGGTTGCGTTCGATCTGTTCCACCGTCACATCTTTGGTTTCAGGAATCAGCCAGTAAGTAATGCCCACAACAGCAATgcgccgc    >  2:618394/1‑137 (MQ=255)
                                                                                                                                              aGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTGCTCGCCTGCCATCAGTTTGCGGTCGATATGTTCCAGCGTGGCATTTTTGGGTTCCGGAGTGAGCCAGAAAGTAAGGCCCACAAACGCAACGTTCAgcgc  >  2:73345/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                              
AGCTCATTCCTCCCAGCAAACCAATTGTTTTCATTCCTCTCTCCTGATGTGTGTTACCGCAACAAGGCTGAATAAGGAACGAGACAAACGCCTCAACGGCCAAGTGCCCAATCTCTATTAACGAAAAAAGGGCCGGATGTACAGCACATCCGGCCCGTGAAATCAGACGCCGATATTTCTCAACTTCTCGCCTGCCATCAGTTTGCGTTCGATATGTTCCAGCGTGACATTTTTGGTTTCCGGAATGAGCCAGAAAGTAATGCCCACAAACGCAATGTTCAGCGC  >  NC_000913/2980602‑2980886

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