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.6% 421.5 / NA 142G168D (GGC→GAC) rhaDrhamnulose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (72/68);  minor base G (0/2);  total (72/70)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 145 total aligned reads displayed.

TTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGG  >  NC_000913/4091278‑4091423
                                                                         |                                                                        
ttCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCt                                                                        >  2:4348492/1‑76 (MQ=255)
 tggTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCtt                                                                       <  2:3391190/74‑1 (MQ=255)
 tCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCtt                                                                       <  2:945901/76‑1 (MQ=255)
  cGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTc                                                                      <  2:2890189/76‑1 (MQ=255)
   gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcc                                                                     <  1:2293327/76‑1 (MQ=255)
   gTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCCGTCTTcc                                                                     >  2:1534400/1‑76 (MQ=255)
    tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                    >  2:4095283/1‑76 (MQ=255)
    tCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTccc                                                                    >  1:1626939/1‑76 (MQ=255)
     ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTcccc                                                                   >  2:2994348/1‑75 (MQ=255)
     ccGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCACACACTCGGTGCTGTCTTCCca                                                                   >  2:4857760/1‑76 (MQ=255)
       gTGCCGGGCTCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCATACACTCGGTGCTGTCTTCCcaca                                                                 <  1:2105153/76‑1 (MQ=255)
       gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                 >  2:3627762/1‑76 (MQ=255)
       gTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCcaca                                                                 >  1:1991956/1‑76 (MQ=255)
        ggCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAg                                                                <  1:1215191/75‑1 (MQ=255)
         gCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGt                                                               >  2:2021334/1‑76 (MQ=255)
          ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                              >  2:3293006/1‑76 (MQ=255)
          ccGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGtt                                                              >  1:84144/1‑76 (MQ=255)
            gggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                            <  1:1973380/76‑1 (MQ=255)
            gggCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTgg                                                            >  1:637728/1‑76 (MQ=255)
               cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                         <  1:2836913/76‑1 (MQ=255)
               cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                         >  2:1642387/1‑76 (MQ=255)
               cACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGa                                                         <  1:1204954/76‑1 (MQ=255)
                aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACCGTTGGCGAg                                                        >  1:4696080/1‑76 (MQ=255)
                aCCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAg                                                        >  1:3705083/1‑76 (MQ=255)
                 ccATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGt                                                       <  1:1655253/76‑1 (MQ=255)
                  cATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTg                                                      >  2:1668905/1‑76 (MQ=255)
                   aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                     <  2:3881104/76‑1 (MQ=255)
                   aTCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGa                                                     <  1:3221547/76‑1 (MQ=255)
                    tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                    <  2:1991955/76‑1 (MQ=255)
                    tCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGaa                                                    <  1:706869/76‑1 (MQ=255)
                     ccACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAg                                                   <  1:1878850/76‑1 (MQ=255)
                      cgcGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                  <  2:3186266/74‑1 (MQ=255)
                      cACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGa                                                  >  1:4820277/1‑76 (MQ=255)
                      cACGGCAAAATGCCAACGCCATCCGGGAACACCACCAGACACTCGCTGCTGTCTTCCCACACACGGCGACTGCAGa                                                  >  2:4400008/1‑76 (MQ=255)
                       aCGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGCTGCTGTCTTCCCACAGTTGGCGAGTGAAGAc                                                 >  1:2102277/1‑76 (MQ=255)
                         ggCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcg                                               >  1:1385783/1‑76 (MQ=255)
                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                              <  1:2327308/76‑1 (MQ=255)
                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                              >  2:2651594/1‑76 (MQ=255)
                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                              <  1:3953892/76‑1 (MQ=255)
                          gCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgc                                              <  1:4730776/76‑1 (MQ=255)
                           cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGACGACcgcg                                             >  2:1172791/1‑76 (MQ=255)
                           cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                             <  1:920438/76‑1 (MQ=255)
                           cAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACcgcg                                             <  1:186881/76‑1 (MQ=255)
                            aaaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGGGCTGTCTGCCCACCGTCGGCGAGTGAAGACCGCgg                                            >  2:423446/1‑76 (MQ=255)
                             aaaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgt                                           <  2:3872435/76‑1 (MQ=255)
                              aaTGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGgtg                                          >  1:2379949/1‑76 (MQ=255)
                                 gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                       <  1:2597507/76‑1 (MQ=255)
                                 gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                       >  2:3283792/1‑76 (MQ=255)
                                 gCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCg                                       <  1:4323519/76‑1 (MQ=255)
                                  ccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGt                                      <  1:3214350/76‑1 (MQ=255)
                                  ccAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTCAAGACCGCGGTGTCGt                                      >  2:3753598/1‑76 (MQ=255)
                                   cAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGtt                                     >  1:4318970/1‑76 (MQ=255)
                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                    >  2:4021957/1‑76 (MQ=255)
                                    aaCGCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGttt                                    <  2:4300096/76‑1 (MQ=255)
                                     aCGCCATCGGGGAATACCACCAGACTCTCGGTGCTGTCTTCCCACAGTTGGCGATTGAAGAGCGCGGTGTCGtttt                                   <  1:813468/76‑1 (MQ=255)
                                      cGCCATCCGGGAATACCACCAGACACTCGGCGCTGTCTTCCCCCAGTCGGCGAGTGAAGACCGCGGGGGCGTTTTc                                  >  2:42485/1‑76 (MQ=255)
                                       gCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                 <  2:2818572/76‑1 (MQ=255)
                                       gCCATCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCa                                 >  2:3143439/1‑76 (MQ=255)
                                        ccATCCCGGGATACCACCAGACACTCACCGCTGTCTTTCCACAGTTGGCGCGTGGAGACCGCGGTGTCGTATTcga                                >  2:1141307/1‑74 (MQ=255)
                                          ctCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              <  1:3372781/75‑1 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCACGt                              >  2:254983/1‑76 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              >  1:1389779/1‑76 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              >  1:2527818/1‑76 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              >  1:2152552/1‑76 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              >  1:3217672/1‑76 (MQ=255)
                                          aTCCGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGt                              <  1:3854169/76‑1 (MQ=255)
                                            ccGGGAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTAc                            <  2:1844995/76‑1 (MQ=255)
                                              gggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACAt                          >  2:2581700/1‑76 (MQ=255)
                                               ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATa                         <  1:2885207/76‑1 (MQ=255)
                                               ggAATACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGCGAAGCCCGCGGTGGCGTTTTCAAGTCCATa                         >  2:4592732/1‑76 (MQ=255)
                                                 aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                       >  2:4530791/1‑76 (MQ=255)
                                                 aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                       <  1:42086/76‑1 (MQ=255)
                                                 aaTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAgg                       >  1:3962105/1‑76 (MQ=255)
                                                  aTACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAggg                      >  2:3189002/1‑75 (MQ=255)
                                                   tACCACCAGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTg                     <  1:3283792/76‑1 (MQ=255)
                                                      caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                  <  1:4520325/76‑1 (MQ=255)
                                                      caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                  <  1:3779249/76‑1 (MQ=255)
                                                      caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                  <  1:1445212/76‑1 (MQ=255)
                                                      caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAgg                  >  1:481080/1‑76 (MQ=255)
                                                      caccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCCCGGTGTCGTTTTCCAGTACATAGGTGCgg                  >  2:603215/1‑76 (MQ=255)
                                                       cccaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                 <  2:1126836/75‑1 (MQ=255)
                                                       accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                 <  1:3293006/76‑1 (MQ=255)
                                                       accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                 >  2:936525/1‑76 (MQ=255)
                                                       accaGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAggg                 <  1:3778951/76‑1 (MQ=255)
                                                         cggACACGCGGTGGTGGCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg               <  2:280703/74‑1 (MQ=255)
                                                          aGACACTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGa              <  1:1098038/76‑1 (MQ=255)
                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg               >  1:1753833/1‑73 (MQ=255)
                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCg               <  2:1753832/73‑1 (MQ=255)
                                                            acacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATc            <  1:4775342/76‑1 (MQ=255)
                                                            acacTCGGTGCCGGCCACCCGCTGTCCGCGATTGAATACCGCGGTGTAGTTTTCAAATTCATAGGTGAGGGCGATc            <  1:162522/76‑1 (MQ=255)
                                                             cacTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCa           >  1:3482478/1‑76 (MQ=255)
                                                             cacTCAGTCCTGTCTGCCCACAGTTGGCGAGTCAAGCCCGCGGTGGCGTTCTCAAGAACATAGGTGCGGGCGATCa           >  2:272957/1‑76 (MQ=255)
                                                               cTCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAgg         <  2:3404069/76‑1 (MQ=255)
                                                                tCGGTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGt        <  1:3313709/76‑1 (MQ=255)
                                                                  ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTCCATAGGGGAGGGCGATCCGGTtg      >  2:3367510/1‑76 (MQ=255)
                                                                  ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg      >  1:3157474/1‑76 (MQ=255)
                                                                  ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg      <  1:1638806/76‑1 (MQ=255)
                                                                  ggTGCTGTCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtg      <  1:2163521/76‑1 (MQ=255)
                                                                     gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTtggtg   <  1:1543553/76‑1 (MQ=255)
                                                                     gctgtCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGGACATAGGTGAGGGCGATCGGGTtgggg   >  1:298341/1‑74 (MQ=255)
                                                                         |                                                                        
TTCGTCCGTGCCGGGCACCATCCACGGCAAAATGCCAACGCCATCCGGGAATACCACCAGACACTCGGTGCTGCCTTCCCACAGTTGGCGAGTGAAGACCGCGGTGTCGTTTTCAAGTACATAGGTGAGGGCGATCAGGTTGGTGG  >  NC_000913/4091278‑4091423

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