Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,069,991 G→A 10.2% 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,9910GA10.2% 115.5 / 4.3 39Q296Q (CAG→CAAcfacyclopropane fatty acyl phospholipid synthase
Reads supporting (aligned to +/- strand):  ref base G (13/22);  new base A (2/2);  total (15/24)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AATGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGATAC  >  minE/1069929‑1070059
                                                              |                                                                    
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTg                                                                  <  1:524713/67‑1 (MQ=255)
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTg                                                                  <  1:1406463/67‑1 (MQ=255)
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTg                                                                  <  1:2112223/67‑1 (MQ=255)
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCt                                                                <  1:3137192/69‑1 (MQ=255)
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCt                                                                <  1:2382301/69‑1 (MQ=255)
aaTGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCt                                                                >  1:3598672/1‑69 (MQ=255)
  tGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGt                                                            >  1:3589274/1‑71 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACTCCAGATTGCTCAGTc                                                           >  1:1144035/1‑70 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTcc                                                          <  1:2674974/71‑1 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTcc                                                          <  1:3478050/71‑1 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTcc                                                          <  1:3689329/71‑1 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTcc                                                          <  1:1200514/71‑1 (MQ=255)
    ttGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTc                                                           >  1:140588/1‑70 (MQ=255)
     tGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTcc                                                          <  1:1472395/70‑1 (MQ=255)
     tGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTc                                                           >  1:3501178/1‑69 (MQ=255)
      gATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAg                                                        >  1:417075/1‑71 (MQ=255)
      gATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCa                                                         <  1:2691102/70‑1 (MQ=255)
            tGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCa                                                         >  1:3625852/1‑64 (MQ=255)
            tGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAAcc                                                  >  1:2888989/1‑71 (MQ=255)
              gATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAAccc                                                 >  1:256486/1‑70 (MQ=255)
              gATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAAccc                                                 >  1:2618871/1‑70 (MQ=255)
                      aataTATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAAGCCCACTTtgt                                         <  1:857549/70‑1 (MQ=255)
                      aataTATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTtgt                                         <  1:1618320/70‑1 (MQ=255)
                      aataTATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTtgt                                         <  1:638818/70‑1 (MQ=255)
                      aataTATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTtgt                                         <  1:3299365/70‑1 (MQ=255)
                          tatTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGc                                                                 <  1:523152/42‑1 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGaa                                                    <  1:3198452/53‑1 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGaa                                                    >  1:157595/1‑53 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGaa                                  <  1:2241332/71‑1 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGaa                                  <  1:2203055/71‑1 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGaa                                  <  1:765835/71‑1 (MQ=255)
                            tttttCCGAACGGTTGCCTGCCCTCTGTACGCCAAATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGaa                                  <  1:1432966/71‑1 (MQ=255)
                                                 ccTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCt             <  1:2632528/71‑1 (MQ=255)
                                                       tACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCa                          >  1:3115183/1‑52 (MQ=255)
                                                       tACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCa                          >  1:509470/1‑52 (MQ=255)
                                                       tACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTAc       <  1:3191098/71‑1 (MQ=255)
                                                       tACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTAc       >  1:3714958/1‑71 (MQ=255)
                                                       tACGCCAGATAGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCa                          >  1:2505174/1‑52 (MQ=255)
                                                            cAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGAtac  <  1:518711/71‑1 (MQ=255)
                                                              gATTTCTCAGTCCAGCGAACCCCACTTTGTGATgg                                    <  1:2266942/35‑1 (MQ=255)
                                                              gATTGCTCAGTCCAGCGAACCCCACTTTGTGATgg                                    <  1:1227568/35‑1 (MQ=255)
                                                              gATTGCTCAGTCCAGCGAACCCCACTTTGTGATgg                                    <  1:1852031/35‑1 (MQ=255)
                                                              gATTGCTCAGTCCAGCGAACCCCACTTTGTGATgg                                    <  1:974162/35‑1 (MQ=255)
                                                              |                                                                    
AATGTTGATCCCTGGATTAATAAATATATTTTTCCGAACGGTTGCCTGCCCTCTGTACGCCAGATTGCTCAGTCCAGCGAACCCCACTTTGTGATGGAAGACTGGCATAACTTCGGTGCTGATTACGATAC  >  minE/1069929‑1070059

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: