Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,358,163 G→A A228V (GCG→GTG)  rhmA ← 2‑keto‑3‑deoxy‑L‑rhamnonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,358,1630GA97.1% 107.7 / ‑6.4 35A228V (GCG→GTG) rhmA2‑keto‑3‑deoxy‑L‑rhamnonate aldolase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (19/15);  minor base C (0/1);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGG  >  NC_000913/2358040‑2358288
                                                                                                                           |                                                                                                                             
tatCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATATAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                <  1:504945/139‑1 (MQ=255)
          tACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAgg                                                                                                      >  1:1001417/1‑139 (MQ=255)
              ttttATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGcc                                                                                                  <  1:595109/139‑1 (MQ=255)
               tttATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgc                                                                                                                    >  2:119473/1‑120 (MQ=255)
               tttATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgc                                                                                                                    <  1:119473/120‑1 (MQ=255)
                    gCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCATGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGcc                                                                                            >  1:941001/1‑139 (MQ=255)
                      gTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 <  2:325697/116‑1 (MQ=255)
                      gTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 >  1:325697/1‑116 (MQ=255)
                       tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                <  2:1240704/116‑1 (MQ=255)
                       tGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTgcg                                                                                                                >  1:1240704/1‑116 (MQ=255)
                        ggCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 >  2:1418956/1‑114 (MQ=255)
                        ggCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 <  1:1418956/114‑1 (MQ=255)
                         gCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCa                                                                                                                        <  2:733637/106‑1 (MQ=255)
                         gCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCa                                                                                                                        >  1:733637/1‑106 (MQ=255)
                             aTTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCa                                                                                   <  2:1057469/139‑1 (MQ=255)
                                    ccTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCtt                                                                            >  1:682917/1‑139 (MQ=255)
                                        aTTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTAcc                                                                        <  1:686130/139‑1 (MQ=255)
                                                  cGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACgg                                                              >  2:1560474/1‑139 (MQ=255)
                                                   gCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACC‑‑GCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATc                                                           >  1:1328257/1‑139 (MQ=255)
                                                    ccAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGAt                                                            >  1:1036447/1‑139 (MQ=255)
                                                       gTCGTTGATCCAGGGCATCGCTTAAGAGCATCGTGTCAACGCCAACAGCTACAAAGTTCGCTCCCTACCCCAGGCATTGCTGATCCATATCAGGAGCCAAAGCCAGAAAACCAGGCGCTTTACCCGCAGAACGGATccg                                                         <  1:1129474/139‑1 (MQ=255)
                                                                   gggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAg                                                                                                       <  1:1406787/81‑1 (MQ=255)
                                                                   gggCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAg                                                                                                       >  2:1406787/1‑81 (MQ=255)
                                                                        tCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCaa                                        >  2:241398/1‑139 (MQ=255)
                                                                        tCGCTGTAGAGCATCGTGTCAACGCCAAAAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCaa                                        >  2:997342/1‑139 (MQ=255)
                                                                         cGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 >  2:1419544/1‑65 (MQ=255)
                                                                         cGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATtgctgc                                                                                                                 <  1:1419544/65‑1 (MQ=255)
                                                                         cGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCaat                                       >  2:215211/1‑139 (MQ=255)
                                                                             gTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAAtt                                   <  2:798586/139‑1 (MQ=255)
                                                                                      cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCt                                                                             >  1:898464/1‑88 (MQ=255)
                                                                                      cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCt                                                                             <  2:898464/88‑1 (MQ=255)
                                                                                      cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCt                               >  2:1577265/1‑134 (MQ=255)
                                                                                      cGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCt                               <  1:1577265/134‑1 (MQ=255)
                                                                                                  aaCAGCGACAAAGTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCg              >  1:377359/1‑139 (MQ=255)
                                                                                                              gTTCGCTCCCCACACCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCggg  <  1:997342/139‑1 (MQ=255)
                                                                                                                           |                                                                                                                             
TATCAATAACTACCTTTTATGCGTGGCCCATTTTTGCCTGATTTAAACATCGCCAGTCGTTGATCCAGGGCATCGCTGTAGAGCATCGTGTCAACGCCAACAGCGACAAAGTTCGCTCCCCACGCCAGGCATTGCTGCGCCATATCAGGAGCCACAGCCAGAAAACCAGCCGCTTTACCCGCAGCACGGATCCGCCGAATACTGGTTTCAATAATTCGCTGCACTTCCGGGTGCCCGGCGTTATCCGGG  >  NC_000913/2358040‑2358288

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