Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,983,421 T→C Q235R (CAG→CGG)  kduI ← 5‑dehydro‑4‑deoxy‑D‑glucuronate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,983,4210TC96.6% 202.8 / ‑0.4 59coding (704/837 nt)kduI5‑dehydro‑4‑deoxy‑D‑glucuronate isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/28);  minor base . (1/1);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAACTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTG  >  NC_000913/2983285‑2983552
                                                                                                                                        |                                                                                                                                   
caACTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCt                                                                                                                                  <  2:47223/140‑1 (MQ=255)
  aCTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACGCCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCg                                                                                                                                <  2:172679/140‑1 (MQ=255)
  aCTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCg                                                                                                                                <  2:663993/140‑1 (MQ=255)
  aCTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCg                                                                                                                                <  1:327037/140‑1 (MQ=255)
     aGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTa                                                                                                                             >  1:803458/1‑140 (MQ=255)
     aGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTa                                                                                                                             <  1:494538/140‑1 (MQ=255)
       cgcAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATg                                                                                                                           >  2:230626/1‑140 (MQ=255)
       cgcAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATg                                                                                                                           >  1:615218/1‑140 (MQ=255)
       cgcAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATg                                                                                                                           <  1:485525/140‑1 (MQ=255)
           aaTCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATc                                                                                                                       >  1:666613/1‑140 (MQ=255)
                       ccACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGAcgc                                                                                                           >  2:410560/1‑140 (MQ=255)
                            tGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCt                                                                                                     <  1:413030/140‑1 (MQ=255)
                                  atatCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCgg                                                                                                <  2:581189/140‑1 (MQ=255)
                                         caAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcccc                                                                                         >  1:169276/1‑140 (MQ=255)
                                         caAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcccc                                                                                         >  1:751798/1‑140 (MQ=255)
                                         caAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcccc                                                                                         >  1:392155/1‑140 (MQ=255)
                                         caAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGcccc                                                                                         >  2:224941/1‑140 (MQ=255)
                                           aaGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCcat                                                                                       >  2:378217/1‑140 (MQ=255)
                                               ccTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCAtg                                                                                   >  2:511560/1‑140 (MQ=255)
                                                      ttCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGaaaa                                                                            >  1:326436/1‑140 (MQ=255)
                                                        cACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                                                          <  2:659876/140‑1 (MQ=255)
                                                            gACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAgg                                                                      <  2:751798/140‑1 (MQ=255)
                                                             aCCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                                                          >  1:155289/1‑135 (MQ=255)
                                                             aCCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACg                                                                          <  2:155289/135‑1 (MQ=255)
                                                             aCCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAgg                                                                      >  2:306934/1‑139 (MQ=255)
                                                             aCCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAgg                                                                      <  1:306934/139‑1 (MQ=255)
                                                             aCCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGc                                                                     >  2:315552/1‑140 (MQ=255)
                                                                aTGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGtc                                                                  <  2:666613/140‑1 (MQ=255)
                                                                 tGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGtca                                                                 >  2:313518/1‑140 (MQ=255)
                                                                     ccAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTcatcat                                                             <  1:268502/140‑1 (MQ=255)
                                                                     ccAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTcatcat                                                             <  2:597222/140‑1 (MQ=255)
                                                                     ccAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTcatcat                                                             <  1:43923/140‑1 (MQ=255)
                                                                      cAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCatcatc                                                            <  1:559650/140‑1 (MQ=255)
                                                                      cAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCatcatc                                                            <  2:392155/140‑1 (MQ=255)
                                                                      cAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCatcatc                                                            <  2:743217/140‑1 (MQ=255)
                                                                        gATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCa                                                          >  2:250207/1‑140 (MQ=255)
                                                                                tataAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa                                                  >  1:742346/1‑140 (MQ=255)
                                                                                tataAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaa                                                  <  2:742346/140‑1 (MQ=255)
                                                                                   aaGCTTTTGTTCCGACACCAGAATGGATCGACCAGCTCGGGGAGATCACCGCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaataga                                              >  2:19046/1‑140 (MQ=255)
                                                                                   aaGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaatag                                               >  1:777772/1‑140 (MQ=255)
                                                                                    aGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaataga                                              >  1:61367/1‑140 (MQ=255)
                                                                                       tttGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaatagaaat                                           >  1:517974/1‑140 (MQ=255)
                                                                                        ttGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTgaaatagaaata                                          <  1:250207/140‑1 (MQ=255)
                                                                                          ggTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATaaa                                        >  1:384524/1‑140 (MQ=255)
                                                                                                gACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCa                                  >  1:294363/1‑140 (MQ=255)
                                                                                                           ggATCGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCg                       <  1:315552/140‑1 (MQ=255)
                                                                                                               cGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTggg                   >  1:782038/1‑140 (MQ=255)
                                                                                                               cGACCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTggg                   <  1:410560/140‑1 (MQ=255)
                                                                                                                 aCCAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGgtg                 >  1:685812/1‑140 (MQ=255)
                                                                                                                  ccAGCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGgtgt                >  1:788232/1‑140 (MQ=255)
                                                                                                                     gCTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGac             <  2:61367/140‑1 (MQ=255)
                                                                                                                      cTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGaca            >  2:606296/1‑140 (MQ=255)
                                                                                                                      cTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGaca            <  1:606296/140‑1 (MQ=255)
                                                                                                                      cTCGGGGAGATCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGaca            <  1:748533/140‑1 (MQ=255)
                                                                                                                          gggAGATCACCGCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCa       <  1:19046/140‑1 (MQ=255)
                                                                                                                             agaTCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATg     <  1:230626/140‑1 (MQ=255)
                                                                                                                               aTCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCa       >  2:637603/1‑136 (MQ=255)
                                                                                                                               aTCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCa       <  1:637603/136‑1 (MQ=255)
                                                                                                                               aTCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGgt   <  2:803458/140‑1 (MQ=255)
                                                                                                                                tCACCGCCCGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGgtg  >  1:518561/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                   
CAACTAGCGCAAATCTTTAACGGCCACATGGTCCATATCATCAAAGACCTGGTTTTCACCGACCATGCCCCAGATAAAGGTATAAGCTTTGGTTCCGACACCGGAATGGATCGACCAGCTCGGGGAGATCACCGCCTGCTCGTTATGCATCACAATATGACGCGTTTCTTGCGGCTGCCCCATCATGTGGAAAACGCAGGCGTCATCATCCATATTGAAATAGAAATAAACTTCCATCCGGCGCTCGTGGGTGTGACACGGCATGGTG  >  NC_000913/2983285‑2983552

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