Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,655,057 T→C T29A (ACG→GCG)  rspA ← bifunctional D‑altronate/D‑mannonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,655,0570TC96.9% 109.0 / NA 32T29A (ACG→GCG) rspAbifunctional D‑altronate/D‑mannonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (18/13);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGC  >  NC_000913/1654923‑1655185
                                                                                                                                      |                                                                                                                                
aCCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATa                                                                                                                              <  1:302871/139‑1 (MQ=255)
     tataGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTc                                                                                                                         >  1:889613/1‑139 (MQ=255)
            aaaCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtg                                                                                                                  >  1:627287/1‑139 (MQ=255)
              aCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                <  2:734501/139‑1 (MQ=255)
              aCTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                <  1:73750/139‑1 (MQ=255)
                 gCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                >  1:104032/1‑136 (MQ=255)
                 gCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAgtggt                                                                                                                <  2:104032/136‑1 (MQ=255)
                      atatCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTa                                                                                                        <  2:838466/139‑1 (MQ=255)
                        atCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAAt                                                                                                      >  2:578287/1‑139 (MQ=255)
                                          aTCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGACCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGcccc                                                                                    >  1:1001720/1‑139 (MQ=255)
                                               gACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGa                                                                                 <  1:381901/137‑1 (MQ=255)
                                               gACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGa                                                                                 >  2:381901/1‑137 (MQ=255)
                                               gACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACa                                                                               >  1:906432/1‑139 (MQ=255)
                                               gACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACa                                                                               >  1:659907/1‑139 (MQ=255)
                                                aCCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAg                                                                              <  2:627287/139‑1 (MQ=255)
                                                             ggACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACtt                                                                 >  2:218932/1‑139 (MQ=255)
                                                                   aGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGCAACAAAAACTTCAGCCt                                                           <  1:578287/139‑1 (MQ=255)
                                                                         tCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGa                                                     >  1:769165/1‑139 (MQ=255)
                                                                                        gCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGAtt                                      >  1:805109/1‑139 (MQ=255)
                                                                                          cACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTc                                                                                                                      <  1:969733/57‑1 (MQ=255)
                                                                                          cACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTc                                                                                                                      >  2:969733/1‑57 (MQ=255)
                                                                                                aaGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCAt                              <  1:721089/139‑1 (MQ=255)
                                                                                                      aCGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                       >  1:682260/1‑44 (MQ=255)
                                                                                                      aCGTCCATTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCt                                                                                                                       <  2:682260/44‑1 (MQ=255)
                                                                                                            aTTGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGa                  <  2:906432/139‑1 (MQ=255)
                                                                                                             ttGAGGGTGGCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGAt                 >  1:742003/1‑139 (MQ=255)
                                                                                                                      gCATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAAt        >  2:939619/1‑139 (MQ=255)
                                                                                                                       cATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                                >  2:194131/1‑82 (MQ=255)
                                                                                                                       cATCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTc                                                                <  1:194131/82‑1 (MQ=255)
                                                                                                                         tCCCCAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTa     <  2:659907/139‑1 (MQ=255)
                                                                                                                            ccAAGGCCCGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTATGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGc  >  2:593951/1‑139 (MQ=255)
                                                                                                                                   ccGCAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTa                                                        >  1:523172/1‑78 (MQ=255)
                                                                                                                                      |                                                                                                                                
ACCTTTATAGAAAAACTGCCAGATATCTTCGATACGGTGCGCATCGCGACCAATAAGCTGCGGACAAAGGTGATCCTGCAAATAAGAGGCCACGGAAAGCTCACGTCCATTGAGGGTGGCATCCCCAAGGCCCGTAATACCGTCCTCAGTGGTGATTTTTAATGTGACGAAATTACGCCCCGGACAGGTAACAAAAACTTCAGCCTTTACGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGC  >  NC_000913/1654923‑1655185

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