Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 729,603 A→C 25.0% P207P (CCA→CCC fabD → malonyl‑CoA‑[acyl‑carrier‑protein] transacylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE729,6030AC25.0% 62.8 / 36.6 52P207P (CCA→CCCfabDmalonyl‑CoA‑[acyl‑carrier‑protein] transacylase
Reads supporting (aligned to +/- strand):  ref base A (19/20);  new base C (7/6);  total (26/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT  >  minE/729535‑729672
                                                                    |                                                                     
gcggcgGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACcag                                                                      >  1:1406939/1‑70 (MQ=255)
gcggcgGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACcag                                                                      >  1:3338842/1‑70 (MQ=255)
          cAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACaa                                                            <  1:3268154/70‑1 (MQ=255)
           aaaaCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAAc                                                          >  1:2450902/1‑71 (MQ=255)
             aaCGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTg                                                        <  1:2129514/71‑1 (MQ=255)
               cgcgcgCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTgg                                                       >  1:1044125/1‑70 (MQ=255)
                 cgcgcTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAg                                                    <  1:3306606/71‑1 (MQ=255)
                      tGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACaa                                                            <  1:864027/58‑1 (MQ=255)
                       gacGTTACCAGTGAGCGTACCGTGTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAt                                              >  1:1392662/3‑71 (MQ=255)
                        ccGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAAtt                                             <  1:2746857/71‑1 (MQ=255)
                        ccGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAtt                                             <  1:2122452/71‑1 (MQ=255)
                         cGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTa                                            <  1:1812753/71‑1 (MQ=255)
                             aCCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCg                                         >  1:1238798/1‑70 (MQ=255)
                               cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa                                      <  1:1333608/71‑1 (MQ=255)
                               cAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa                                      <  1:2352206/71‑1 (MQ=255)
                                aGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa                                     >  1:2593168/1‑71 (MQ=255)
                                 gTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa                                     <  1:2515254/70‑1 (MQ=255)
                                  tGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAAtt                                             <  1:2169469/61‑1 (MQ=255)
                                   gAGCGTACCGTCTCACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGaa                                               >  1:3487844/1‑58 (MQ=255)
                                    aGCGTACCGTCTCACTGTGCGCTGATTAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 <  1:3452148/71‑1 (MQ=255)
                                    aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 <  1:3280600/71‑1 (MQ=255)
                                    aGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCAc                                 >  1:3185462/1‑71 (MQ=255)
                                      cGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa                                      <  1:2991794/64‑1 (MQ=255)
                                      cGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCa                                  <  1:1384179/68‑1 (MQ=255)
                                        tACCGTCTCACTGTGCGCTGATGAAACCCGCCGCCGACaaa                                                           >  1:2925660/1‑41 (MQ=38)
                                        tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGc                                          >  1:39123/1‑58 (MQ=255)
                                        tACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaa                                      <  1:1199241/62‑1 (MQ=255)
                                        tACCGTCTCACTGTGCGCTGATGAAACCAGAGCCGACAAACTGGCAGTAGAATTAGCGAAAAt                                    <  1:2549133/63‑1 (MQ=255)
                                               tcACTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCttt                             <  1:787061/64‑1 (MQ=255)
                                               tcACTGTGCGCTGATGAAACCAGCAGCCGACAAACTgg                                                       <  1:4068101/38‑1 (MQ=255)
                                               tcACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAAt                                    >  1:65853/1‑57 (MQ=255)
                                                 aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:2497391/1‑64 (MQ=255)
                                                 aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:1289678/1‑64 (MQ=255)
                                                 aCTGTGCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCAcc                     >  1:3199168/1‑70 (MQ=255)
                                                 aCTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGaaaa                                     >  1:1772401/1‑54 (MQ=255)
                                                    gtgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa                   >  1:3947650/1‑69 (MQ=255)
                                                     tgCGCTGATGAAACCCGCAGCCGACAAACTGGCAGTa                                                  <  1:1417049/37‑1 (MQ=255)
                                                     tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCa                       >  1:914748/1‑64 (MQ=255)
                                                     tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAAc                  <  1:3863467/69‑1 (MQ=255)
                                                     tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACa                 <  1:3969633/70‑1 (MQ=255)
                                                     tgCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAg                >  1:73272/1‑71 (MQ=255)
                                                        gcTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCaa                   >  1:1555058/1‑65 (MQ=255)
                                                         cTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTcc            >  1:209615/1‑71 (MQ=255)
                                                            atgaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtg          >  1:2966419/1‑70 (MQ=255)
                                                             tgaAACCCGCAGCCGACAAACTGGCAGTAGAATTAg                                           >  1:778030/1‑36 (MQ=255)
                                                              gaAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttg       >  1:1321206/1‑71 (MQ=255)
                                                               aaaCCCGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt      <  1:3945791/71‑1 (MQ=255)
                                                               aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTaa                           >  1:2130742/1‑50 (MQ=255)
                                                               aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt      >  1:872711/1‑71 (MQ=255)
                                                               aaaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtg          <  1:3028040/67‑1 (MQ=255)
                                                                aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtgttgt      <  1:2063715/70‑1 (MQ=255)
                                                                aaCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCtggtgt      >  1:4127857/1‑70 (MQ=255)
                                                                 aCCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGa    <  1:919787/71‑1 (MQ=255)
                                                                  ccAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGaa   >  1:2105943/1‑71 (MQ=255)
                                                                   cagcagCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAt  <  1:2810698/71‑1 (MQ=255)
                                                                    |                                                                     
GCGGCGGGCGCAAAACGCGCGCTGCCGTTACCAGTGAGCGTACCGTCTCACTGTGCGCTGATGAAACCAGCAGCCGACAAACTGGCAGTAGAATTAGCGAAAATCACCTTTAACGCACCAACAGTTCCTGTTGTGAAT  >  minE/729535‑729672

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