Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 3,186,703 T→A 19.4% E389V (GAA→GTA) ‡ U712_16400 ← 3‑ketoacyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068813,186,7030TA19.4% 116.8 / 3.4 53E389V (GAA→GTA) ‡U712_164003‑ketoacyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base T (23/19);  new base A (2/9);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCG  >  CP006881/3186579‑3186840
                                                                                                                             |                                                                                                                                          
gCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGc                                                                                                                               <  1:152555/138‑1 (MQ=255)
   cGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCt                                                                                                                           >  1:15019/1‑139 (MQ=255)
             aTCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgc                                                                                                                <  1:37553/139‑1 (MQ=255)
              tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATaca                                                                                                                                    >  5:52541/1‑119 (MQ=255)
              tcctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATaca                                                                                                                                    <  6:52541/119‑1 (MQ=255)
               cctccGGTGTATTCTTTTGATCGGATGTGACATCCTCTATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgcct                                                                                                              <  4:287099/139‑2 (MQ=255)
                ctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgccg                                                                                                              <  4:101136/138‑1 (MQ=255)
                ctccGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgccg                                                                                                              <  8:211586/138‑1 (MQ=255)
                   cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgccgcc                                                                                                            <  2:272152/137‑1 (MQ=255)
                   cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgccgcc                                                                                                            <  1:225404/137‑1 (MQ=255)
                   cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCGCCCAGT‑ccgccgccc                                                                                                           >  1:146550/1‑137 (MQ=255)
                   cGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTAAAATACACCGGCCGCTCCCATT‑ccgccgccg                                                                                                           >  3:291635/1‑138 (MQ=255)
                    ggTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGGCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGAt                                                                                                         <  7:76400/139‑1 (MQ=255)
                     gtgtATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATATTTCATATACACCGGCCGCTCCCATT‑GCGCCGCCGATa                                                                                                        >  2:37553/1‑139 (MQ=255)
                     gtgtATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACCCCGGCCGCTCCCATT‑CCGCCGCCGATa                                                                                                        <  6:76726/139‑1 (MQ=255)
                     gtgtATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgc                                                                                                                >  2:240200/1‑131 (MQ=255)
                     gtgtATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑ccgc                                                                                                                <  1:240200/131‑1 (MQ=255)
                          taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACAt                                                                                                   >  8:68626/1‑139 (MQ=255)
                          taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCGTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGGCGCCGTGacacct                                                                                                   >  4:259350/1‑137 (MQ=255)
                          taCATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTATTGTACGTCAGCCGCTTTTTTTGCCA‑TTTTTTTTTCCCCCATTAGCATAAATAAAATACACCGGCCGCTCCCAGT‑CCGCCGCCGATaaccat                                                                                                   >  5:2027/1‑134 (MQ=255)
                               ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCT‑TTTGGTTTTCCCCCTTTAGCATAATTCAAATAAACCGGCCGCTGCCATT‑CCGCCGCCGATACACATGGtga                                                                                              >  5:95019/1‑137 (MQ=255)
                               ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCCTT‑CCGCCGCCGATACACATGGTCa                                                                                              >  6:264696/1‑139 (MQ=255)
                               ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACCTGGTCa                                                                                              >  2:8072/1‑139 (MQ=255)
                               ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTACCCTT‑CCGCCGCCGATACACCTGGTCa                                                                                              >  5:181910/1‑139 (MQ=255)
                               ttGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCCCCCATT‑CCGCCGCCGAGACACCAGGTcc                                                                                              >  1:139391/1‑138 (MQ=255)
                                tGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACCCATGGGCAc                                                                                             >  6:116120/1‑139 (MQ=255)
                                  aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCCAATACACCGGCCGCTCCCCTT‑CCGCCGCCCATACACATGGTCACgc                                                                                           >  4:289854/1‑138 (MQ=255)
                                  aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGa                                                                                           >  7:159356/1‑139 (MQ=255)
                                  aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATAAAAATACACCGGCCGCTACCCTT‑CCGCCGCCCCTACACAAgggcccgc                                                                                           >  1:143475/1‑133 (MQ=255)
                                  aTCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCAAAATTAAAATAAACAGGCCGCGCCCATT‑CAGCCGCCGATAAACAGGGTcgcgg                                                                                           >  7:279001/1‑138 (MQ=255)
                                       aTGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCg                                                                                                                             <  8:155453/101‑1 (MQ=255)
                                       aTGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCg                                                                                                                             >  7:155453/1‑101 (MQ=255)
                                       aTGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGACTcc                                                                                       >  3:279590/1‑138 (MQ=255)
                                             cATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTTACA‑TTTTTTTTTACCCCTTTAGCATAATAAAAATACACCGGCCGCTGCCAAT‑CCGCCGCCGATACACATGTTCCCGACAGCAAAtt                                                                                  >  6:100521/1‑137 (MQ=255)
                                                ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCC‑TTTGTTTTTGCCCCGTTAGCATAATTCAAAAACACCGGCCGCTCCCAAT‑ACGCCGCCGATACAAAAGGTGCCGCCTCCAAAATGGTCg                                                                             >  8:199178/1‑139 (MQ=255)
                                                ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTACCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCCATACACATAGTCCCGACACCAAATTgagag                                                                             >  4:236009/1‑135 (MQ=255)
                                                ccTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑GTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCAAT‑CCGCCGCCGAGACACAAGGGAACGGCTCCAAATTGCTCg                                                                             >  4:51688/1‑139 (MQ=255)
                                                                 ccgccTTTTTGTACGCCAGCCGCTTTTTTTGCCA‑TTTGTTTTACCCCCTTTAGCATAATACAAAAACCCCGGCCGCACCCATT‑CCGCCGCCGAAACACATGGACACGACTCCAAATTGCTCGTTTCTCCGtttcatttc                                                            <  1:8072/139‑1 (MQ=255)
                                                                     cTTTTTGTACGTCCGCCGCTTTTTTTGCCA‑TTTTTTTTCCCCCCTTTCGCAAAATTCAAAACCCCCGGCCGCACCCATT‑CCCCCGCCCATCCACAGGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATgg                                                        <  8:159356/139‑1 (MQ=255)
                                                                       ttttGTACGTCCGCCGCTTTTTTTGCAA‑TTTGTTTTTCCCACTTTAGCAAAATACAAATACACCGCCCGCTCCCATC‑CCGCCGCCTATACACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGa                                                       <  7:273487/138‑1 (MQ=255)
                                                                       ttttGTACGTCAGCCGCTTTTTTTGCCT‑TTTGTTTTTCCCCATTTACAATAATAAAAATACCCCGGCCGCCCCCATT‑CCGCCGCCGATACACTTGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGa                                                       <  7:68626/138‑1 (MQ=255)
                                                                        tttGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGAATAATTCAAATACCCCGGCCGCTCCAATC‑CCGCCGCCGATACACTTGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAt                                                      <  3:140818/138‑1 (MQ=255)
                                                                            aaCGTCCGCCGCTTTTTTTGCCT‑TTTGTTTTTCCCCCTTTCGAAAAATTAAAACACCCCGGCCGCCCCCATT‑CCGCCCCCGATACACTTGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                 <  4:146241/138‑1 (MQ=255)
                                                                             ccgccAGCCGC‑TTTTTTGCCATTTTGTTTTCCCCCCTTACGAATAATAAAAAAACCCCGCCCCCTCCAAAC‑CCCCCGCCAATACACAGGCCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                 <  2:24475/134‑1 (MQ=255)
                                                                             acgcccgCCGCTTTTTTTGCCT‑TTTGTGTCACCCCCTTACGCAAAATACAAAAACCCCGGCCGCACCCATT‑CCGCCGCCTATACACATGGTCACCACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtga                                                <  6:95019/133‑1 (MQ=255)
                                                                             aCGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTAACCATAATTAAAATCCACCGGCCGCTCCCATT‑CCGCCGCATATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtg                                                 <  7:248633/138‑1 (MQ=255)
                                                                              ctgcccccGATTTTTTTGCCA‑TTTGTTTTCCCCCGTTTAGAATAATACAAAACCCCCGGCCGCTCCCATG‑CCGCCGCCTATACAAATGGCCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGAtgagtgag                                               <  1:220668/133‑1 (MQ=255)
                                                                                   gCCGTTTTTTTTCCTA‑TTTGTTATGCCCCCTTTATCAAAATTAAAAAACACCGGCGGCTCCAATC‑CCGCCGCCGAAACACATGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAgtgt                                          <  7:67563/139‑1 (MQ=255)
                                                                                    ccGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTACCATAATACAAAACCCCCGGCCGCTCCCATT‑CCGCCGCCTATACACATGGCCACGACTCCAAATTGCTCGTTTCTCCGGTTCATTTCATGGATGAGTGAGAGTGTa                                         <  4:279590/139‑1 (MQ=255)
                                                                                                               cccTTTAGAAAAATAAAAAAACCCCGGCCGCGCCCATC‑CCGCCGCCTAACCACAGGGGCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg               <  6:181910/139‑1 (MQ=255)
                                                                                                               cccTATACCAAAATTAAAAAACACCGGCCG‑ACCCATTCCCCCCGGCCATAAATATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg               <  3:51688/139‑1 (MQ=255)
                                                                                                               ccatttACCATAATAAAACTACCCCGGCCGCTCCAATT‑CCGCCTCCTAAACACAGGGCCACGACTCCAAATTGCTCGATTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCgg               <  6:136473/136‑1 (MQ=255)
                                                                                                                        aTAATTAAATTCCCCCGCCCGCTCACATT‑CCTCCGCCCAAACACTGGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAg      <  3:102765/139‑1 (MQ=255)
                                                                                                                            ttCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAgcgc   >  6:133439/1‑138 (MQ=255)
                                                                                                                            ttCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGcgcg  >  1:203338/1‑139 (MQ=255)
                                                                                                                             |                                                                                                                                          
GCTCGTCAGTGAAATCCTCCGGTGTATACATTTGATCGTATGTGACATCCTCGATTAAAAATCCGCCGCCTTTTTGTACGTCAGCCGCTTTTTTTGCCA‑TTTGTTTTTCCCCCTTTAGCATAATTCAAATACACCGGCCGCTCCCATT‑CCGCCGCCGATACACATGGTCACGACTCCAAATTGCTCGTTTCTCCGTTTCATTTCATGGATGAGTGAGAGTGTAAGCTTTGTGCCCGTACAGCCAAGCGGGTGTCCAAGCGCG  >  CP006881/3186579‑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: ‑