Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 130,093 C→T 9.5% D91D (GAC→GAT yadE → predicted polysaccharide deacetylase lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE130,0930CT 9.5% 184.2 / 7.8 63D91D (GAC→GATyadEpredicted polysaccharide deacetylase lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base T (2/4);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TTGGCTTTGGCAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTA  >  minE/130023‑130161
                                                                      |                                                                    
ttggctttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGAc                                                                      >  1:2866230/1‑71 (MQ=255)
  ggctttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACg                                                                     <  1:3981178/70‑1 (MQ=255)
   gctttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACgg                                                                    <  1:3446979/70‑1 (MQ=255)
   gctttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACg                                                                     <  1:530127/69‑1 (MQ=255)
    ctttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGtt                                                                  <  1:1770974/71‑1 (MQ=255)
      ttggcAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGtt                                                                  >  1:3778680/1‑69 (MQ=255)
      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)
                  ttATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGtt                                                                  >  1:4206516/1‑57 (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:981645/70‑1 (MQ=255)
                                 aTCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTaa                                      <  1:1423164/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)
                                 aTCTCGAGCCGGTTCAGGGGCGAAAAAAAGTTGAAGACGGtt                                                                  <  1:3169746/42‑1 (MQ=255)
                                  tctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTcaac                                                    <  1:3828747/55‑1 (MQ=255)
                                  tctcGAGCCGGTTCAGGGGCCACAAAAAGTTGAAGACGGtt                                                                  <  1:2310416/41‑1 (MQ=255)
                                   ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATTGTTTAGGTGACCTca                                                      <  1:2398037/52‑1 (MQ=255)
                                   ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTAGGTGACCTca                                                      <  1:1269985/52‑1 (MQ=255)
                                   ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                     <  1:4348823/69‑1 (MQ=255)
                                   ctcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                     >  1:4361102/1‑69 (MQ=255)
                                    tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGAc                                                                      >  1:2800304/1‑35 (MQ=255)
                                    tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGATGGTTTAGGTGACCTCAACAAGCCGCTGAGt                                        >  1:889589/1‑65 (MQ=255)
                                    tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                     >  1:582571/1‑68 (MQ=255)
                                    tcGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAAt                                     >  1:3317072/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:4112166/71‑1 (MQ=255)
                                             ttCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGtt                         <  1:4281021/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:4635588/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: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:2493107/1‑70 (MQ=255)
                                                                     aCGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCta  >  1:1863368/1‑70 (MQ=255)
                                                                      |                                                                    
TTGGCTTTGGCAAAGGTTTTATCGATAAAGGTCATCTCGAGCCGGTTCAGGGGCGACAAAAAGTTGAAGACGGTTTGGGCGACCTCAACAAGCCGCTGAGTAATCAGAACTTAGTTACCTGGAAAGATACGCCGGTCTA  >  minE/130023‑130161

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: