Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,321,698 G→A 13.4% T93I (ACC→ATC)  cdh ← CDP‑diacylglycerol phosphotidylhydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,321,6980GA13.4% 141.6 / 12.3 52T93I (ACC→ATC) cdhCDP‑diacylglycerol phosphotidylhydrolase
Reads supporting (aligned to +/- strand):  ref base G (11/34);  new base A (3/4);  total (14/38)
Fisher's exact test for biased strand distribution p-value = 3.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

TTTTTTGCTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTGCAGTGGGCCATTAA  >  minE/2321634‑2321767
                                                                |                                                                     
ttttttGCTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGATca                                                                    >  1:643831/1‑68 (MQ=255)
 tttttGCTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTcaaca                                                                 >  1:597777/1‑70 (MQ=255)
    ttGCTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAg                                                              <  1:1055227/70‑1 (MQ=255)
       cTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACt                                                          <  1:868845/71‑1 (MQ=255)
         cATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTc                                                                     <  1:1898195/58‑1 (MQ=255)
          aTAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAg                                                              <  1:1994996/64‑1 (MQ=255)
          aTAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAg                                                              <  1:364095/64‑1 (MQ=255)
          aTAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAg                                                              <  1:1186693/64‑1 (MQ=255)
          aTAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACttt                                                        >  1:527648/1‑70 (MQ=255)
           tAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTcaaca                                                                 <  1:311917/60‑1 (MQ=255)
           tAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTc                                                       <  1:1484017/70‑1 (MQ=255)
             aaaTCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGt                                                    <  1:634787/71‑1 (MQ=255)
             aaaTCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGt                                                    <  1:1878162/71‑1 (MQ=255)
                      gcCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGATCAACAAAGGACTTTCAg                                                     >  1:1548325/1‑61 (MQ=255)
                        cTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatac                                         <  1:246963/71‑1 (MQ=255)
                        cTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGATCAACAAAGGACTTTCAGTACCGTTAatac                                         <  1:1687463/71‑1 (MQ=255)
                         tGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatac                                         >  1:1642558/1‑70 (MQ=255)
                           ccaagccaaCCAAAAGAAGTTCGGCGTTGAAGGATCGATCAACAAAGGACTTTCAGTACCGTTAatacga                                       >  1:558970/1‑70 (MQ=255)
                                ccaaccaaAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatacga                                       <  1:732519/65‑1 (MQ=255)
                                ccaaccaaAAGAAGTTCGGCGTTGAAGGATCGATCAACAAAGGACTTTCAGTACCGTTAATACGATACGtt                                 <  1:1125281/71‑1 (MQ=255)
                                   accaaAAGAAGTTCGGCGTTGAAGGATCGGTcaac                                                                  <  1:180962/35‑1 (MQ=255)
                                      aaaaGAAGTTCGGCGTTGAAGGATCTATCAACAAAGGACTTTCAGTACCGTTAatacgatac                                    <  1:1027673/62‑1 (MQ=255)
                                       aaaGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGtt                                 <  1:1920885/64‑1 (MQ=255)
                                       aaaGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGtt                                 <  1:2149632/64‑1 (MQ=255)
                                        aagaagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCaa                         <  1:1283887/71‑1 (MQ=255)
                                        aagaagTTCGGCGTTGAAGGATCGATCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCaa                         <  1:573516/71‑1 (MQ=255)
                                         agaagTTCGGCGTTGAAGGATCGGTCAACAAAGGa                                                            >  1:1592139/1‑35 (MQ=255)
                                         agaagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTAcc                                                 >  1:311565/1‑46 (MQ=255)
                                         agaagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTAcc                                                 >  1:2008223/1‑46 (MQ=255)
                                         agaagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCaa                         >  1:1839810/1‑70 (MQ=255)
                                          gaagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatacga                                       <  1:1895565/55‑1 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatacg                                        >  1:870667/1‑53 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCaa                         >  1:1376627/1‑68 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAg                      <  1:536455/71‑1 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAg                      <  1:1016894/71‑1 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAg                      <  1:478916/71‑1 (MQ=255)
                                           aagTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAg                      <  1:239706/71‑1 (MQ=255)
                                            agTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:302115/71‑1 (MQ=255)
                                            agTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:433389/71‑1 (MQ=255)
                                             gTTCGGCGTTGAAGGATCGGTCAACAACGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:1464055/70‑1 (MQ=255)
                                                 ggCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTg                <  1:186460/71‑1 (MQ=255)
                                                 ggCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTg                <  1:1506270/71‑1 (MQ=255)
                                                     ttGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAatac                                         <  1:1086514/42‑1 (MQ=255)
                                                     ttGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:1805898/62‑1 (MQ=255)
                                                     ttGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:2078006/62‑1 (MQ=255)
                                                     ttGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATAAGATACGTTGGCATCAACAGATATTGCAGt            <  1:647538/71‑1 (MQ=255)
                                                       gAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCa                             >  1:2115089/1‑52 (MQ=255)
                                                           gATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGc                              <  1:530885/47‑1 (MQ=255)
                                                           gATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGc                              <  1:182680/47‑1 (MQ=255)
                                                           gATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGc                              <  1:692898/47‑1 (MQ=255)
                                                           gATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGc                              <  1:7059/47‑1 (MQ=255)
                                                             tCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACa                       >  1:1185462/1‑52 (MQ=255)
                                                               ggTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAGCAGATATTGCAGTGGGCCATTaa  <  1:1938018/71‑1 (MQ=255)
                                                               ggTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTGCAGTGGGCCAtt    >  1:1186865/1‑69 (MQ=255)
                                                               ggTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTGCAGTGGGCCATTaa  <  1:1563513/71‑1 (MQ=255)
                                                               ggTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTGCAGTGGGCCATTaa  <  1:1824114/71‑1 (MQ=255)
                                                                gTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGa                     <  1:715948/51‑1 (MQ=255)
                                                                |                                                                     
TTTTTTGCTCATAAAATCACGCGCCTGCCAAGCCAACCAAAAGAAGTTCGGCGTTGAAGGATCGGTCAACAAAGGACTTTCAGTACCGTTAATACGATACGTTGGCATCAACAGATATTGCAGTGGGCCATTAA  >  minE/2321634‑2321767

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: