Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,102 C→T 6.6% G94G (GGC→GGT yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,1020CT 6.6% 184.8 / 2.7 60G94G (GGC→GGTyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (2/2);  total (32/28)
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.63e-01

AAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGCGCC  >  minE/130035‑130170
                                                                   |                                                                    
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)
                     aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTaa                                               <  1:1423164/70‑1 (MQ=255)
                     aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTaa                                               <  1:981645/70‑1 (MQ=255)
                     aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                              <  1:4014585/71‑1 (MQ=255)
                     aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                              <  1:557386/71‑1 (MQ=255)
                     aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                              <  1:4447047/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:1146204/1‑70 (MQ=255)
                                             aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                       >  1:3714219/1‑70 (MQ=255)
                                             aaaaaGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGa                       >  1:4635588/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:742447/1‑71 (MQ=255)
                                                     gAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGt              >  1:3010345/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)
                                                                 ggCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGCGcc  <  1:2594258/71‑1 (MQ=255)
                                                                 ggCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGCGcc  <  1:3976027/71‑1 (MQ=255)
                                                                   |                                                                    
AAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTATAACGCGCC  >  minE/130035‑130170

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: