Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,622,669 T→C 48.3% A659A (GCA→GCG maeB ← fused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,622,6690TC48.3% ‑3.7 / 50.6 58A659A (GCA→GCGmaeBfused malic enzyme predicted oxidoreductase and predicted phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base T (8/22);  new base C (14/14);  total (22/36)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTCGACGACGGGCAG  >  minE/1622604‑1622723
                                                                 |                                                      
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCTTGCATTTCACCATCAATCATCAGTTCTGGTGct                                                     <  1:765807/69‑2 (MQ=255)
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:3315921/69‑1 (MQ=255)
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:2637811/69‑1 (MQ=255)
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:2289683/69‑1 (MQ=255)
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:2095085/69‑1 (MQ=255)
tCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:1568488/69‑1 (MQ=255)
 cGTTGCGAAGCGCTTCCACCAGCGCTGCATCTCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCa                                                     <  1:1689988/68‑1 (MQ=255)
       gattCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1990363/65‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1046667/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:876570/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:661809/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:642557/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:377405/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:346821/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:2559381/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1773804/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1161134/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1303701/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1578146/68‑1 (MQ=255)
       gAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:1408644/68‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCCTCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:3110789/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:1747056/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:476329/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:472298/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:1368314/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:3312059/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:2934315/69‑1 (MQ=255)
        aaTCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:2699653/69‑1 (MQ=255)
                       cgcTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGCGCACGTTcc                                               <  1:247934/52‑1 (MQ=255)
                               cGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCa                      <  1:1721568/69‑1 (MQ=255)
                               cGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCa                      <  1:535483/69‑1 (MQ=255)
                               cGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCa                      <  1:508861/69‑1 (MQ=255)
                               cGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCa                      <  1:2684222/69‑1 (MQ=255)
                                       aTTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCt                                             <  1:2798192/38‑1 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:1417259/1‑69 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:888139/1‑69 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:1044296/1‑69 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:654242/1‑69 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:1378898/1‑69 (MQ=255)
                                          tCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:3125163/1‑69 (MQ=255)
                                          tCACAATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGc                       >  1:908545/1‑57 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGa                          >  1:2818083/1‑50 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:1525657/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:2890939/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:1681721/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:1599035/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:1463517/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:1439118/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:234532/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:2363110/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:3237907/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:3207708/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTc              >  1:3138881/1‑62 (MQ=255)
                                              cATCAATCATCAGTTCTGGCGCACGTTCCCTGACCAGTTACAGCGCCTGACGCATTTTGCTc              >  1:1777109/1‑62 (MQ=255)
                                              cATCAAACATCAGTTCTGGCGCACGTTCCCTGAGCAGTTCCAGCGCCTGACGCATTTTGCt               >  1:196340/1‑61 (MQ=255)
                                                    tcatcaGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTCGACGACGGGCAg  <  1:2567502/68‑1 (MQ=255)
                                                     catcaGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTcgac           >  1:1156044/1‑58 (MQ=255)
                                                     catcaGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTCGACGACGGGCAg  <  1:2343724/67‑1 (MQ=255)
                                                                 |                                                      
TCGTTGCGAATCGCTTCCACCAGCGCTGCATCGCCGTGCATTTCACCATCAATCATCAGTTCTGGTGCACGTTCCCTGACCAGTTCCAGCGCCTGACGCATTTTGCTCGACGACGGGCAG  >  minE/1622604‑1622723

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