Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,099 G→A 6.9% L93L (TTG→TTA yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,0990GA 6.9% 176.2 / 2.7 58L93L (TTG→TTAyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (30/24);  new base A (2/2);  total (32/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 = 9.12e-01

TTGGCAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGC  >  minE/130029‑130167
                                                                      |                                                                    
ttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTg                                                                      >  1:9277/1‑71 (MQ=255)
 tggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTgg                                                                     >  1:3061023/1‑71 (MQ=255)
   gtaaaGGGTTTATCGATAAAGGCCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTaggt                                                                   <  1:52198/69‑2 (MQ=255)
   gcAAAGGTTTTATCTATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTg                                                                      <  1:4323570/68‑1 (MQ=255)
      aaGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGAc                                                                >  1:1589090/1‑71 (MQ=255)
       aGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGAcc                                                               >  1:82188/1‑71 (MQ=255)
         gTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTc                                                             <  1:1014767/71‑1 (MQ=255)
            ttATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTcaa                                                           >  1:2349897/1‑70 (MQ=255)
                 gATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGcc                                                     <  1:4513521/71‑1 (MQ=255)
                 gATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGc                                                      >  1:3747990/1‑70 (MQ=255)
                   tAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGc                                                   >  1:3421374/1‑71 (MQ=255)
                        gTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCt                                                  <  1:2409300/67‑1 (MQ=255)
                        gTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAg                                               >  1:1386807/1‑70 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTaggt                                                                   <  1:1952741/47‑2 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTaa                                            <  1:1423164/70‑1 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTaa                                            <  1:981645/70‑1 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           <  1:4014585/71‑1 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           <  1:4447047/71‑1 (MQ=255)
                           aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           <  1:557386/71‑1 (MQ=255)
                            tctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTcaac                                                          <  1:3828747/55‑1 (MQ=255)
                             ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATTGTTTAGGTGACCTca                                                            <  1:2398037/52‑1 (MQ=255)
                             ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTAGGTGACCTca                                                            <  1:1269985/52‑1 (MQ=255)
                             ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           >  1:4361102/1‑69 (MQ=255)
                             ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           <  1:4348823/69‑1 (MQ=255)
                              tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTAGGTGACCTCAACAAGCCGCTGAGt                                              >  1:889589/1‑65 (MQ=255)
                              tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           >  1:3317072/1‑68 (MQ=255)
                              tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                           >  1:582571/1‑68 (MQ=255)
                                gAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTcaaca                                                         >  1:2644679/1‑52 (MQ=255)
                                gAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGaa                                      >  1:1489154/1‑71 (MQ=255)
                                   ccGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACtt                                   <  1:2953094/71‑1 (MQ=255)
                                       ttCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTa                                             <  1:3000970/57‑1 (MQ=255)
                                       ttCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGtt                               <  1:4281021/71‑1 (MQ=255)
                                       ttCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGtt                               <  1:4112166/71‑1 (MQ=255)
                                        tCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTa                              >  1:3364461/1‑71 (MQ=255)
                                          aGGGGCGACATAAAGTTGAAGACGGTTTGGGCGACCTcaa                                                           <  1:3111856/40‑1 (MQ=255)
                                           ggggCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCt                           <  1:4101462/71‑1 (MQ=255)
                                            gggCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTg                          <  1:113533/71‑1 (MQ=255)
                                                gACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGaaa                      <  1:1201703/71‑1 (MQ=255)
                                                gACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGaaa                      <  1:3496690/71‑1 (MQ=255)
                                                 aCAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTgg                         >  1:3924519/1‑67 (MQ=255)
                                                  cAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                    <  1:3979571/71‑1 (MQ=255)
                                                   aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGa                                                >  1:1792451/1‑42 (MQ=255)
                                                   aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                    >  1:4635588/1‑70 (MQ=255)
                                                   aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                    >  1:1146204/1‑70 (MQ=255)
                                                   aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                    >  1:3714219/1‑70 (MQ=255)
                                                    aaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATa                  >  1:1850052/1‑71 (MQ=255)
                                                      aaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTgg                         >  1:2663256/1‑62 (MQ=255)
                                                          tGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCgg            >  1:4307321/1‑71 (MQ=255)
                                                           gAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGt           >  1:3010345/1‑71 (MQ=255)
                                                           gAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGt           >  1:742447/1‑71 (MQ=255)
                                                             aGATGGTTTAGGTGACCTCAACAAGCCGCTGAGTAATCAGAACTTAATAAcc                            >  1:1387247/1‑52 (MQ=255)
                                                             aGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGa                                       >  1:614320/1‑41 (MQ=255)
                                                               aCGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTa                                  <  1:3941622/44‑1 (MQ=255)
                                                               aCGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCta        >  1:1863368/1‑70 (MQ=255)
                                                               aCGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCta        >  1:2493107/1‑70 (MQ=255)
                                                                 ggTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCtata      <  1:889037/70‑1 (MQ=255)
                                                                 ggTTTGCGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATaa     >  1:4388/1‑71 (MQ=255)
                                                                  gTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAg                     <  1:433865/54‑1 (MQ=255)
                                                                  gTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATaa     >  1:950932/1‑70 (MQ=255)
                                                                   tttGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTa                                  <  1:4230709/40‑1 (MQ=255)
                                                                    ttGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCt                           >  1:1156639/1‑46 (MQ=255)
                                                                    ttGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATa                  >  1:895202/1‑55 (MQ=255)
                                                                    ttGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAAcgc  >  1:2473090/1‑71 (MQ=255)
                                                                      |                                                                    
TTGGCAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGC  >  minE/130029‑130167

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: