Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 2,021,817 C→T V139I (GTT→ATT)  ykhA ← acyl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,021,8170CT97.6% 131.2 / ‑3.8 41V139I (GTT→ATT) ykhAacyl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/19);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

TTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGA  >  CP000730/2021698‑2021948
                                                                                                                       |                                                                                                                                   
tttttAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACAcc                                                                                                                <  1:134707/141‑1 (MQ=255)
  tttAGGCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAg                                                                                                              >  2:163578/1‑141 (MQ=255)
   ttAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTTTTTTTTTTTTTCTTGATATACACCAgg                                                                                                             >  2:304767/1‑141 (MQ=255)
        ctctAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTACTGATATACACCAGCTACAt                                                                                                        >  2:371581/1‑141 (MQ=255)
           tAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATATAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGtt                                                                                                     <  1:304471/141‑1 (MQ=255)
           tAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGtt                                                                                                     >  1:92334/1‑141 (MQ=255)
               atGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAgg                                                                                                 <  2:101715/141‑1 (MQ=255)
                tGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGt                                                                                                >  2:20600/1‑141 (MQ=255)
                 gTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTACATGTTTAGGtt                                                                                               <  1:221598/141‑1 (MQ=255)
                  ttGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGttt                                                                                              >  2:246427/1‑141 (MQ=255)
                  ttGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGttt                                                                                              >  1:45884/1‑141 (MQ=255)
                    gCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGc                                                                                            >  2:235672/1‑141 (MQ=255)
                      aCTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCTATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCt                                                                                          >  2:208730/1‑141 (MQ=255)
                         tGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtca                                                                                       <  1:140679/141‑1 (MQ=255)
                         tGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtca                                                                                       <  1:266565/141‑1 (MQ=255)
                          gAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtcat                                                                                      <  2:248192/141‑1 (MQ=255)
                           aGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTtcatc                                                                                     >  1:114033/1‑141 (MQ=255)
                              tAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCCATATCACCTTCTTGATATACACCAGGTACATGTTTAGGTTTGCCTTCatcatc                                                                                  <  2:23841/141‑1 (MQ=255)
                                 atatTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTaa                                                                               <  1:271884/141‑1 (MQ=255)
                                    ttgtaTTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAGCCATTTTTCAATATCATCTTCTGGATATGCACCAGGTACATGTTTAGGTTTGCCTTCATAATCTAACGc                                                                            <  2:113198/137‑1 (MQ=255)
                                      cAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGcaa                                                                          >  1:130550/1‑141 (MQ=255)
                                       aaTTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGcaac                                                                         <  1:285096/141‑1 (MQ=255)
                                             ttgtttACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGt                                                                   <  1:68463/141‑1 (MQ=255)
                                                   aCTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTa                                                             <  2:266484/141‑1 (MQ=255)
                                                    cTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCCGGTCGATGCTTTGGCTTACGTTCATCATCTAACGCAACAAATGTTAAGTaa                                                            >  1:163558/1‑141 (MQ=255)
                                                     tttCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAAc                                                           >  2:22427/1‑141 (MQ=255)
                                                             cGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCa                                                   >  1:5113/1‑141 (MQ=255)
                                                              gACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGccg                                                  >  1:256287/1‑139 (MQ=255)
                                                              gACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAg                                                  >  1:256312/1‑141 (MQ=255)
                                                                     ttttCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACCCCCGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAGCAAATGTTAAGTAACTTAATGCATCTAAATc                                           <  2:312990/141‑1 (MQ=255)
                                                                        tcttcGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATg                                        <  2:98430/141‑1 (MQ=255)
                                                                                 aaCACGTTGTGGTGCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGttatt                               <  1:4401/141‑1 (MQ=255)
                                                                                  acacGTAGTGGAGCTGGATCATAAAACCATTTTTCAATATGATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCTTTCATCCTCTAACGCAACAAATGTTAAATAACTTAATGCAGCTAAATCATGCTTGttatta                              >  2:60632/1‑141 (MQ=255)
                                                                                   cacGTTGTGGTGCTGTATCATAAAACCATTTTTCAATTTCATCTTCTGGATATACACCATGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTaa                             >  2:47800/1‑141 (MQ=255)
                                                                                           ggtgCTGTATCATAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCg                     >  1:296162/1‑141 (MQ=255)
                                                                                               cggTATCATAAAACCAGTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACTCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCaa                 <  2:73345/139‑1 (MQ=255)
                                                                                                      aTAAAACCCTTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAAt          <  1:321408/141‑1 (MQ=255)
                                                                                                      aTAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAAt          <  1:292074/141‑1 (MQ=255)
                                                                                                       tAAAACCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAAtt         >  1:296508/1‑141 (MQ=255)
                                                                                                        aaaaCCATTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAAttt        <  1:304767/141‑1 (MQ=255)
                                                                                                              aTTTTTCAATATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGa  >  1:327467/1‑141 (MQ=255)
                                                                                                                       |                                                                                                                                   
TTTTTAGTCTCTAATATGTTGCACTTGAGCTAAATATTCAATTGTTTGTTTACTTTCAATGCGACGTGCTTTTCTTCGTTCAACACGTTGTGGTGCTGTATCATAAAACCATTTTTCAACATCATCTTCTGGATATACACCAGGTACATGTTTAGGTTTGCCTTCATCATCTAACGCAACAAATGTTAAGTAACTTAATGCAGCTAAATCATGCTTGTTATTAAATACATCGTCAATTCTAATTTGAACGA  >  CP000730/2021698‑2021948

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