Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,069,991 G→A 41.3% Q296Q (CAG→CAA cfa → cyclopropane fatty acyl phospholipid synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,069,9910GA41.3% 27.7 / 53.3 46Q296Q (CAG→CAAcfacyclopropane fatty acyl phospholipid synthase
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base A (6/13);  total (22/24)
Fisher's exact test for biased strand distribution p-value = 7.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGATA  >  minE/1069929‑1070058
                                                              |                                                                   
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGc                                                                <  1:1673151/68‑1 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTc                                                          >  1:3640914/1‑70 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTcc                                                         <  1:3948230/71‑1 (MQ=255)
     tGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTcc                                                         <  1:3622334/70‑1 (MQ=255)
      gATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAg                                                       >  1:2374054/1‑71 (MQ=255)
      gATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAg                                                       >  1:3695413/1‑71 (MQ=255)
        tCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCg                                                     <  1:3078362/71‑1 (MQ=255)
          ccTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGaa                                                   >  1:3206275/1‑71 (MQ=255)
               aTTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAAcccc                                               <  1:265505/70‑1 (MQ=255)
                       atatatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGa                                      <  1:2249932/71‑1 (MQ=255)
                       atatatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGa                                      <  1:2433358/71‑1 (MQ=255)
                         atatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGAt                                     >  1:1474823/1‑70 (MQ=255)
                         atatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATg                                    >  1:1250107/1‑71 (MQ=255)
                         atatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATg                                    >  1:2481617/1‑71 (MQ=255)
                                tCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGAc                              <  1:46872/70‑1 (MQ=255)
                                tCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACt                             >  1:927992/1‑71 (MQ=255)
                                    aaCGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTtgt                                        <  1:3765886/56‑1 (MQ=255)
                                     aCGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCAt                        <  1:2673283/71‑1 (MQ=255)
                                        gTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTTGCATAAc                     >  1:579101/1‑71 (MQ=255)
                                        gTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAAc                     >  1:2848247/1‑71 (MQ=255)
                                        gTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAAc                     <  1:2792848/71‑1 (MQ=255)
                                        gTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTAGTGATGGAAGACTGGCATAAc                     <  1:94969/71‑1 (MQ=255)
                                          tgcctgccCTCTGTACGCCAGATTGGTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACtt                   >  1:3477367/1‑71 (MQ=255)
                                          tgcctgccCTCTGTACGCCAGATTGCTCAGTCCAg                                                       >  1:3391092/1‑35 (MQ=255)
                                          tgcctgccCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTtgt                                        <  1:99536/50‑1 (MQ=255)
                                          tgcctgccCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACtt                   >  1:1768587/1‑71 (MQ=255)
                                           gcctgccCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACtt                   <  1:2066561/70‑1 (MQ=255)
                                           gcctgccCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACtt                   <  1:2006838/70‑1 (MQ=255)
                                           gcctgccCTCTGTACGCCAAATTGCGCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACtt                   <  1:2086606/70‑1 (MQ=255)
                                             ctgccCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCgg                <  1:2200208/71‑1 (MQ=255)
                                             ctgccCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCgg                <  1:732075/71‑1 (MQ=255)
                                             ctgccCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCgg                <  1:4228128/71‑1 (MQ=255)
                                             ctgccCTCTGTACGCCAGATTGCTCAGTCAAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCgg                <  1:930987/71‑1 (MQ=255)
                                                cccTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTg              <  1:2237695/70‑1 (MQ=255)
                                                cccTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTg              <  1:3652237/70‑1 (MQ=255)
                                                    ctGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGAt         >  1:1721024/1‑71 (MQ=255)
                                                    ctGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGAt         >  1:494177/1‑71 (MQ=255)
                                                    ctGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGAt         >  1:719611/1‑71 (MQ=255)
                                                    ctGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGa          >  1:476492/1‑70 (MQ=255)
                                                     tGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATTGAAGACTGGCATAACTTCGGTGCTGa          <  1:3822621/69‑1 (MQ=255)
                                                         cGCCAAATTACTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGa    <  1:1598425/71‑1 (MQ=255)
                                                           ccAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGAta  >  1:3154516/1‑71 (MQ=255)
                                                           ccAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGAta  >  1:1292952/1‑71 (MQ=255)
                                                            cAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACt                             >  1:2489130/1‑43 (MQ=255)
                                                            cAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACt                             >  1:1623052/1‑43 (MQ=255)
                                                            cAGATTGCCCAGTCCAGCGAACCCCACTTTGTGATGGAAGACt                             >  1:3550658/1‑43 (MQ=255)
                                                              |                                                                   
AATGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGATA  >  minE/1069929‑1070058

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: