Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,186,704 C→A 20.0% E389* (GAA→TAA)  U712_16400 ← 3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,186,7040CA20.0% 122.2 / 4.4 58E389* (GAA→TAA) U712_164003‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base C (28/17);  new base A (3/10);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 2.48e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCG  >  CP006881/3186578‑3186840
                                                                                                                              |                                                                                                                                           
tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGc                                                                                                                               <  7:145040/139‑1 (MQ=255)
tGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGc                                                                                                                               <  2:53663/139‑1 (MQ=255)
    cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑aa                                                                                                                                          <  6:186965/125‑1 (MQ=255)
    cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑aa                                                                                                                                          >  5:186965/1‑125 (MQ=255)
    cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTA‑AAATA‑CACCGGCCGC‑t                                                                                                                          >  8:334345/1‑139 (MQ=255)
     gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑Tc                                                                                                                         >  3:175864/1‑139 (MQ=255)
     gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCAATTTGTTTTTCCCCCTTTAGCATAATTC‑CAATA‑CACCTGCCGC‑Tc                                                                                                                         >  5:125909/1‑139 (MQ=255)
     gTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCCATTAAAAATCCGCAGCCATTTTTTACGTCAGCCGCTTTTTTTTCCATTTTTTTTTAACCCATTAGGTAAATTC‑AAAAA‑AACCCGcccgga                                                                                                                          >  6:211729/1‑135 (MQ=255)
              aTCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTGTTTGTACGGCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTccgc                                                                                                                <  5:140572/139‑1 (MQ=255)
              aTCCTCCGGTGGATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTccgc                                                                                                                <  4:142238/139‑1 (MQ=255)
               tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CCCCGGCCCC‑TCCCAATccgcc                                                                                                               >  3:319642/1‑139 (MQ=255)
               tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTccgcc                                                                                                               <  1:337518/139‑1 (MQ=255)
                    cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCCTTCCGCCGCCGa                                                                                                          >  7:337918/1‑139 (MQ=255)
                    cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGa                                                                                                          <  8:249522/139‑1 (MQ=255)
                     ggTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CCCCGGCCGC‑TCCCATTCCGCCGCCGAt                                                                                                         <  1:301676/139‑1 (MQ=255)
                      gtgtATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACAGGCCGC‑TCCACTTCCGCCGCCGATa                                                                                                        >  6:140572/1‑139 (MQ=255)
                         tataCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCAATTGTTTTTCCACCTTTAGCATAATTA‑AAATA‑CACCGGCCGCTTCCAATTCGGCCGCCCAAaca                                                                                                      >  7:6531/1‑139 (MQ=255)
                                ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑aaa                                                                                                                                         >  1:172750/1‑98 (MQ=255)
                                ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑aaa                                                                                                                                         <  2:172750/98‑1 (MQ=255)
                                ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AATTA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATAACCATGGGCa                                                                                              >  3:118583/1‑139 (MQ=255)
                                ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCCTTCCGCCGCCGATACACATGGtac                                                                                              >  4:274478/1‑137 (MQ=255)
                                ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCa                                                                                              >  8:300430/1‑139 (MQ=255)
                                 tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGc                                                                                                                               <  3:297293/106‑1 (MQ=255)
                                 tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGc                                                                                                                               >  4:297293/1‑106 (MQ=255)
                                 tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTA‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCAc                                                                                             >  3:41200/1‑139 (MQ=255)
                                   aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACCCATGGTCACGa                                                                                           >  3:324849/1‑139 (MQ=255)
                                        aTGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCCTTCCGCCGCCGATACACATGGTCACGACTcc                                                                                       >  8:56908/1‑138 (MQ=255)
                                                  cTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTACGCCGCCGATACCCATGGGCACGCCTCCAAATTGGTCGt                                                                            >  2:257072/1‑139 (MQ=255)
                                                   tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCCTTCCGCCGCCGCTACACATGGTCCCGACTCCAAATTGCTCGtt                                                                           >  1:42681/1‑139 (MQ=255)
                                                   tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGCCGACAAATTTCTCGtt                                                                           >  4:234688/1‑139 (MQ=255)
                                                   tCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTACCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACAAGGTCACGACTCCAAATTGCTCGtt                                                                           >  5:345258/1‑139 (MQ=255)
                                                                  ccccttttttGAACGCCA‑CCCCTTTTTTTTCCATTTTGTTTTCCCCCATTAATAAAAATA‑AAAAACCCCCGGCCGC‑CCCAATTCCGCCGCCTAAACACAGGGACACGACTCCAAATTGCTCGTTTCTCCGtttcatttc                                                            <  4:118583/134‑1 (MQ=255)
                                                                      cTTTTTGTACGTCCGCCGCTTTTTTTGCCATTTGTTTTCCCCCATTTCCAAAAATAA‑AAATA‑CCCCGGCCGC‑GCCCATCCCGCCGCCCATACACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                        <  6:356839/139‑1 (MQ=255)
                                                                      cTTTTTGAACGCCCCCCGCTTTTTTCCCCATTTGTTTTCCCCCCTTAACAAAAATAA‑AAAAC‑CCCCGGCCGC‑ACCCAGTCCGCCGCCTAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                        <  3:234688/139‑1 (MQ=255)
                                                                             tACGTCCGGGGCTTTTTTTGCCATTTGTTTTACCCAATTTATCTTAATTC‑AAATA‑CACCGGCCGC‑ACCCTTACCCCCGCCTATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                 <  6:345258/139‑1 (MQ=255)
                                                                             tACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCAATCCGCCGCCGATACACATGGTCACGACTCCAAATTGCGCCTTTCTCCGGTTGATTTCATGGAtgcgt                                                  >  5:356839/1‑138 (MQ=255)
                                                                              ccgccAGCCGCTTTTTTTCCCTTTTTTTTTACCCCATTTAGCAAAATTA‑AAAAC‑CCCCGCCCGC‑CCCCATTCCGCCGCCTAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                <  8:337918/135‑1 (MQ=255)
                                                                               cGTCAGCCTCTTTTTTTGCAATTTTTTTTTCCCCCTTTAGCATAATTA‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACAAATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                <  7:1515/138‑1 (MQ=255)
                                                                               cGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAAAA‑CACCGGCCGC‑GCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                <  7:56908/138‑1 (MQ=255)
                                                                                     ccGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCAATACGCCGCCCATACACATGGTACCGCCTCCAAATTGCTACTTTCTCCGTTTCATTTAATGTAtgagtgatacgtt                                          >  4:363467/1‑132 (MQ=255)
                                                                                     ccGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGAAAAATTA‑AAATA‑CACGGGCGGC‑TCCGATTCCGCGTCCTATACACTTGTACACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTa                                         <  5:140726/139‑1 (MQ=255)
                                                                                         tttttttGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTccgccgccg                                                                                                           >  5:262900/1‑69 (MQ=255)
                                                                                         tttttttGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTccgccgccg                                                                                                           <  6:262900/69‑1 (MQ=255)
                                                                                           tttttGCAATTTGTTTTACCCCCTTTAGCATAATTA‑AAAAA‑CCCCGGCCGC‑TCCCATTCCGCCGCCTAAACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCtt                                    <  3:374059/138‑1 (MQ=255)
                                                                                               tGCCATTTGTTTTTCCCCCTTTAGCATAATTA‑AAATA‑CACCGGCCGC‑TCCCTTTCCGCCGCCGATACAAATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTtgtg                                <  7:334345/138‑1 (MQ=255)
                                                                                                     ttGTTTTGCCCCCTTTCGAATAATCA‑AAAAC‑CCCCGGCCGC‑TCCCTCCCCGCCCCCCATACAACTGCCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTa                          <  4:41200/138‑1 (MQ=255)
                                                                                                      tttttatCCCCATTTATCTTATATC‑AAAAA‑TCCCGGCCGC‑CCACATTCCGCCTCCCATAACCCTGCTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                         <  2:194450/136‑1 (MQ=255)
                                                                                                      tGTTTTTCCCCCTTTAGCAAAATTC‑AAAAC‑CCCCGGCCGC‑CCCCATGCCGCCGCCGAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                         <  4:175864/138‑1 (MQ=255)
                                                                                                      tGTTTTCCCCCCTTTAGCATAATTC‑AAAAA‑CACCGCCCGC‑TCCCATTCCGCCGCCGATACACTTGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                         <  7:300430/138‑1 (MQ=255)
                                                                                                       ttttttCCCCCCTTAACAAAAAAAAAAAAA‑CCCCGCCCGC‑CCCAATTCCGCCCCCTAAACACTTGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTAc                         <  3:274478/137‑1 (MQ=255)
                                                                                                                        ataataC‑AAATA‑CACCGGCCGC‑TCCCATACCGCCGCCGAAACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAg      <  7:50264/139‑1 (MQ=255)
                                                                                                                          aaTTC‑AAATA‑CACCGGCCGC‑TGCCATTCCGCCGCCGCTACACATGGTCACGACTCCAAATTGCTCGGTTGTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTTCCCGTACAGCCCAGCGGGGGTCCAAgcg    >  2:354454/1‑139 (MQ=255)
                                                                                                                          aaTTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAAGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAgcg    >  7:84834/1‑139 (MQ=255)
                                                                                                                            ttC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGCTACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGGTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGcgcg  >  2:126863/1‑139 (MQ=255)
                                                                                                                            ttC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCCcgcgcg  >  8:243142/1‑139 (MQ=255)
                                                                                                                            ttC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAgcgc   >  2:17016/1‑138 (MQ=255)
                                                                                                                            ttC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGGTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACCGCCCAGCGGGTGTCCAcgcgcg  >  2:63592/1‑139 (MQ=255)
                                                                                                                            ttC‑AAATA‑CACCGGCCGC‑GCCCCTTCCGCCGCCGCTACACATGGGCCCGACGCCCAATAGCTCGTTTCTCCGGTTCGTTTCCTGGATGCGTGAGAGGGTAAGCTTTGGGCCCGTACAGCCAATCCGGTGTCCCcgcgcg  >  6:53973/1‑139 (MQ=255)
                                                                                                                              |                                                                                                                                           
TGCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCATTTGTTTTTCCCCCTTTAGCATAATTC‑AAATA‑CACCGGCCGC‑TCCCATTCCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCG  >  CP006881/3186578‑3186840

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