Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,091,351 C→T G168D (GGC→GAC)  rhaD ← rhamnulose‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,091,3510CT98.0% 424.2 / ‑3.3 151G168D (GGC→GAC) rhaDrhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (70/77);  total (73/78)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

Only 100 of 158 total aligned reads displayed.

TCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGG  >  NC_000913/4091279‑4091423
                                                                        |                                                                        
tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCtt                                                                       <  2:2223650/76‑1 (MQ=255)
 cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                      >  2:4111916/1‑76 (MQ=255)
   cccgggccgggCCCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                    <  2:625849/75‑1 (MQ=255)
    ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCca                                                                   <  2:1561890/76‑1 (MQ=255)
     ggTGCGGGGCCCCAGCCGCGGCAAAATGGCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcac                                                                  <  2:408853/75‑1 (MQ=255)
     cGTGCCGGGCACCATCCACGGCAAAATGCTAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcac                                                                  >  1:4983848/1‑76 (MQ=255)
       tGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                >  1:2294815/1‑76 (MQ=255)
        gcggggCAGCATCCGCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                               <  2:2406412/73‑1 (MQ=255)
        gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCCGACACTCGGTGCCGTCTTCCCACAGt                                                               >  2:5464616/1‑76 (MQ=255)
        gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                               >  1:5410241/1‑76 (MQ=255)
        gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                               >  1:4777087/1‑76 (MQ=255)
         ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                              <  1:5253162/76‑1 (MQ=255)
         ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                              <  1:2349918/76‑1 (MQ=255)
          cGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTg                                                             <  1:4049887/76‑1 (MQ=255)
            ggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCAACAGTTGGc                                                           <  2:2423627/76‑1 (MQ=255)
             gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACCGTTGGCg                                                          >  2:3781493/1‑76 (MQ=255)
             gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACCGTTGGCg                                                          >  1:1436518/1‑76 (MQ=255)
             gCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGCTGGCg                                                          >  2:165484/1‑76 (MQ=255)
              cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                         <  2:3859130/76‑1 (MQ=255)
               aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                        <  1:53895/76‑1 (MQ=255)
               aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGCTGGCggg                                                        >  2:597144/1‑74 (MQ=255)
                ccacccacGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                       <  2:5092443/76‑1 (MQ=255)
                ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                       <  2:4156771/76‑1 (MQ=255)
                ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                       <  1:1203383/76‑1 (MQ=255)
                  aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                     >  2:3518040/1‑76 (MQ=255)
                  aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                     >  1:3494457/1‑76 (MQ=255)
                   tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                    <  2:3094511/76‑1 (MQ=255)
                   tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                    >  1:3986468/1‑76 (MQ=255)
                    ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                   <  2:3592919/76‑1 (MQ=255)
                    ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                   <  2:951543/76‑1 (MQ=255)
                     cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                  <  1:3924049/76‑1 (MQ=255)
                     cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                  >  2:4468601/1‑76 (MQ=255)
                      ccGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                 <  1:2222843/75‑1 (MQ=255)
                       cGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                <  2:2091298/76‑1 (MQ=255)
                       cGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGAcc                                                <  1:299364/76‑1 (MQ=255)
                        ggCAAAATGCCAACGCCGTCCGGGAATACCACCAGACACTCTGTGCTGTCTTCCCCCAGTTGGCGAGTGAAGACcg                                               <  2:391477/76‑1 (MQ=255)
                        ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                               <  2:584495/76‑1 (MQ=255)
                        ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                               >  1:1205895/1‑76 (MQ=255)
                         gCAAAATGCCAACGCCATCCGGGAATCCCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                              <  2:1233639/76‑1 (MQ=255)
                         gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGTTGCTGTCTTCCCACCGTTGGCGAGCGAAGACcgc                                              >  1:1010668/1‑76 (MQ=255)
                          gtgaaTGCCAACGCCAGCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                             <  1:5786410/73‑1 (MQ=255)
                          cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                             >  2:5630846/1‑76 (MQ=255)
                            aaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgt                                           <  2:3634785/76‑1 (MQ=255)
                            aaaTGCAAGCGCCATCCGGGAAGACCACCAGACACTCGGTGCTGGCTTCCCACAGTTGGCGAGTGAAGACCGCGgt                                           <  2:3656865/76‑1 (MQ=255)
                             aaTGCCAACGCGATCGGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                          <  2:5733041/76‑1 (MQ=255)
                             aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                          <  1:943798/76‑1 (MQ=255)
                             aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                          <  2:2294815/76‑1 (MQ=255)
                             aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                          >  2:4134018/1‑76 (MQ=255)
                              aTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtgt                                         >  2:725957/1‑76 (MQ=255)
                               tGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                        >  1:3002164/1‑76 (MQ=255)
                               tGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTc                                        >  2:4812949/1‑76 (MQ=255)
                                gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                       <  1:3422758/76‑1 (MQ=255)
                                   aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                    >  2:1058321/1‑76 (MQ=255)
                                    aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTTGCGAGAGAAGACCGCGGCGTCGtttt                                   >  2:309329/1‑76 (MQ=255)
                                    aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtttt                                   >  1:2576597/1‑76 (MQ=255)
                                    aCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttct                                   >  1:3189988/1‑74 (MQ=255)
                                     cGCCATCCGGTAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTc                                  <  2:3942934/76‑1 (MQ=255)
                                        cATCCGGGAATACCACCAGACGCTGGGTGCGGTCTGCCCACAGTTGGCGAGTGAAGGCCGCGGTGTCGTTTTCAAg                               <  2:184904/76‑1 (MQ=255)
                                          tCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAgca                             >  1:631003/1‑74 (MQ=255)
                                           ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGGGGCGTTTTCAAGTAc                            >  2:2885575/1‑76 (MQ=255)
                                           ccGGGAATACCACCAGACACTCGGTGCTGCATTCCCACCGTTGGCGCGTCACGCCCGCGGTGCCGTTTTCACGTAc                            >  2:464954/1‑76 (MQ=255)
                                            cGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACa                           >  2:3414413/1‑76 (MQ=255)
                                             gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                          >  1:27359/1‑76 (MQ=255)
                                              ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                         <  1:3464760/76‑1 (MQ=255)
                                               gAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAg                        <  1:4810779/76‑1 (MQ=255)
                                                aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                       >  2:877791/1‑76 (MQ=255)
                                                 aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGt                      >  1:4957248/1‑76 (MQ=255)
                                                  tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                     >  2:4816125/1‑76 (MQ=255)
                                                  tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                     <  2:1361906/76‑1 (MQ=255)
                                                   accaccAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGa                    >  1:4360334/1‑76 (MQ=255)
                                                    ccaccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAg                   <  1:1657777/76‑1 (MQ=255)
                                                      cacaGACACTCGGTGCTGTCTTCCCAGAGTTGGCGAGTGAAGACCGCGTTGTCGTTTTCAAGTACATAGGTGAggg                 <  2:4522072/74‑1 (MQ=255)
                                                      accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                 <  1:3766216/76‑1 (MQ=255)
                                                       ccaTACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGCCCGCGGTGTCGTTTTCAAGTACCTAGGTGAGGGc                <  1:5019799/76‑1 (MQ=255)
                                                       ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGc                <  2:1550570/76‑1 (MQ=255)
                                                       ccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGc                <  2:5243770/76‑1 (MQ=255)
                                                        caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg               <  2:839847/76‑1 (MQ=255)
                                                        caGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg               <  2:4355124/76‑1 (MQ=255)
                                                         aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGGGAGGGCGa              >  1:2226531/1‑76 (MQ=255)
                                                         aGACACTCGGTGCGGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGa              >  1:1568580/1‑76 (MQ=255)
                                                           acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc            <  1:2383163/76‑1 (MQ=255)
                                                           acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc            <  2:4360333/76‑1 (MQ=255)
                                                           acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc            <  1:388862/76‑1 (MQ=255)
                                                             acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGACGACCGCGGTGTCGTTTTCAAGTACATAGGTGACGGCGATCAg          >  2:1255017/1‑76 (MQ=255)
                                                             acTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAg          >  2:2093021/1‑76 (MQ=255)
                                                               tctgTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt        <  2:1575054/73‑1 (MQ=255)
                                                               tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt        <  1:725957/76‑1 (MQ=255)
                                                               tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt        <  1:5644322/76‑1 (MQ=255)
                                                               tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt        >  1:1942989/1‑76 (MQ=255)
                                                               tCGGTGCTGTCTTCCCACAGGTGGCGAGTGAAGACCGCGGTGTCGTCTTCAAGTACCTCGGTGAGGGCGATCAggg        >  1:3232493/1‑75 (MQ=255)
                                                                cGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGtt       <  1:671497/76‑1 (MQ=255)
                                                                 ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg      >  1:3784231/1‑76 (MQ=255)
                                                                 ggTGCTGTCTTCCCACAGTGGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg      <  2:3420160/76‑1 (MQ=255)
                                                                  gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg     >  1:5461004/1‑76 (MQ=255)
                                                                  gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg     >  1:2615053/1‑76 (MQ=255)
                                                                  gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg     >  1:2035244/1‑76 (MQ=255)
                                                                  gTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtgg     <  1:107978/76‑1 (MQ=255)
                                                                   tgctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggt    >  1:2127728/1‑76 (MQ=255)
                                                                    gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg   >  2:690822/1‑76 (MQ=255)
                                                                     ctgcCTACCCACAGTCGGCGAGTGAAGACCGCGGTGTCGTTATCAAGCACATTGGGGAGGGAGACCAGGAtggggg  >  1:1772604/1‑76 (MQ=255)
                                                                        |                                                                        
TCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGG  >  NC_000913/4091279‑4091423

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