Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,131,666 G→A 19.3% L428L (CTG→TTG)  deaD ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,131,6660GA19.3% 133.4 / 24.9 62L428L (CTG→TTG) deaDATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base G (29/21);  new base A (5/7);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTCAGCAGTGC  >  minE/2131597‑2131735
                                                                     |                                                                     
cACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGt                                                                      >  1:2255154/1‑71 (MQ=255)
   gTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGcc                                                                   >  1:38507/1‑71 (MQ=255)
   gTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGcc                                                                   >  1:1966840/1‑71 (MQ=255)
   gTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGcc                                                                   >  1:3728700/1‑71 (MQ=255)
          aCGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCa                                                            <  1:3902688/71‑1 (MQ=255)
              aTCGGCTCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGt                                                        >  1:847/1‑71 (MQ=255)
              aTCGGCGCATCTGGCGGTACGTTCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGt                                                                      >  1:94822/1‑57 (MQ=255)
              aTCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGt                                                                      >  1:2289604/1‑57 (MQ=255)
              aTCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGt                                                        >  1:1196590/1‑71 (MQ=255)
              aTCGGCGCATCTGGCGGTACGATAGAGTACGTTCACCCTGTGCCATTTTCAGCAGt                                                                      >  1:2600128/1‑56 (MQ=255)
               tCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGtt                                                       <  1:2931115/71‑1 (MQ=255)
               tCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGtt                                                       <  1:2244944/71‑1 (MQ=255)
               tCGGCCCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGtt                                                       <  1:940071/71‑1 (MQ=255)
                     cATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCgaga                                                 <  1:545980/71‑1 (MQ=255)
                     cATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCgaga                                                 <  1:1509129/71‑1 (MQ=255)
                      aTCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTc                                                     <  1:3214255/66‑1 (MQ=255)
                      aTCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTc                                                     >  1:3825005/1‑66 (MQ=255)
                       tCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGt                                                                      <  1:3091685/48‑1 (MQ=255)
                         tGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagcag                                                                       >  1:2992617/1‑45 (MQ=255)
                         tGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTcagcag                                                                       >  1:680809/1‑45 (MQ=255)
                          ggCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAg                                            >  1:3870425/1‑71 (MQ=255)
                           gCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATcc                                              <  1:167322/68‑1 (MQ=255)
                           gCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATcc                                              <  1:2057070/68‑1 (MQ=255)
                              gTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGctct                                        <  1:589651/71‑1 (MQ=255)
                              gTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGctct                                        <  1:2962949/71‑1 (MQ=255)
                              gTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGctct                                        >  1:2635340/1‑71 (MQ=255)
                              gTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGctct                                        >  1:3519060/1‑71 (MQ=255)
                                aCGATCAGAGTACGTTCACCCTGTGCCATTTTcagcag                                                                       >  1:273539/1‑38 (MQ=255)
                                aCGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCtt                                       >  1:3315015/1‑70 (MQ=255)
                                aCGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCtt                                       >  1:2604645/1‑70 (MQ=255)
                                  gATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCg                                                    >  1:3614054/1‑55 (MQ=255)
                                   aTCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCAcc                                   <  1:3035699/71‑1 (MQ=255)
                                     cAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGAt                                                >  1:2467453/1‑56 (MQ=255)
                                     cAGAGTACGTTCACCCTGTGCCATATTCAGCAGTGCCGCAGCCAGAGTTTCgag                                                  <  1:3666902/54‑1 (MQ=255)
                                       gagTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTc                                <  1:1157744/70‑1 (MQ=255)
                                        agTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTcag                              >  1:2554413/1‑71 (MQ=255)
                                          tACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGt                                                        <  1:1838530/43‑1 (MQ=255)
                                          tACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTcagca                            <  1:3132068/71‑1 (MQ=255)
                                           aCGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGc                                                              >  1:1838571/1‑36 (MQ=255)
                                            cGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGt                          >  1:698148/1‑71 (MQ=255)
                                             gTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTc                         <  1:2602995/71‑1 (MQ=255)
                                              ttCACCCTGTGCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCg                        <  1:1195112/71‑1 (MQ=255)
                                                  cccTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTc                                                     <  1:1052183/38‑1 (MQ=255)
                                                   ccTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGttt                                                      <  1:3946300/36‑1 (MQ=255)
                                                   acTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTc                                >  1:1936170/2‑58 (MQ=255)
                                                     tgtgCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGAt                                                >  1:3317491/1‑40 (MQ=255)
                                                     tgtgCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCgg                       >  1:3637630/1‑65 (MQ=255)
                                                       tgCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCa                                     <  1:2325357/49‑1 (MQ=255)
                                                       tgCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCt                     >  1:1952655/1‑65 (MQ=255)
                                                       tgCCATTTTCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAAt                 <  1:1594290/69‑1 (MQ=255)
                                                         ccATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTg             <  1:1297583/71‑1 (MQ=255)
                                                          cATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGc            <  1:3939210/71‑1 (MQ=255)
                                                           aTTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGCTCCAGCTCTTCAcc                                   >  1:3445668/1‑47 (MQ=255)
                                                            ttttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCt           >  1:287203/1‑70 (MQ=255)
                                                             tttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCtt                                       >  1:822317/1‑41 (MQ=255)
                                                             tttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCa                                     >  1:3186932/1‑43 (MQ=255)
                                                             tttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTcagc                             >  1:1004010/1‑51 (MQ=255)
                                                             tttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAAtt                >  1:955529/1‑64 (MQ=255)
                                                              ttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTcag        >  1:220314/1‑71 (MQ=255)
                                                              ttCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTca         <  1:3359401/70‑1 (MQ=255)
                                                              ttCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCAc                                    >  1:2269208/1‑43 (MQ=255)
                                                              ttCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCAc                                    >  1:3771663/1‑43 (MQ=255)
                                                               tCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTcag        <  1:2141248/70‑1 (MQ=255)
                                                               tCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAACAGTCGGCTGAATTTTGCTcag        >  1:3361875/1‑70 (MQ=255)
                                                               tCAGCAATGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTcag        <  1:2393026/70‑1 (MQ=255)
                                                                 agcagTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCa                                     >  1:2879556/1‑39 (MQ=255)
                                                                 agcagTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTcagca                            <  1:1870572/48‑1 (MQ=255)
                                                                    agTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTCAGCAGTgc  <  1:3586104/71‑1 (MQ=255)
                                                                     |                                                                     
CACGTTTCGGACGCATCGGCGCATCTGGCGGTACGATCAGAGTACGTTCACCCTGTGCCATTTTCAGCAGTGCCGCAGCCAGAGTTTCGAGATCCAGCTCTTCACCTTCAGCAGTCGGCTGAATTTTGCTCAGCAGTGC  >  minE/2131597‑2131735

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