Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,217,503 A→G 62.1% F410F (TTT→TTC prc ← carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,217,5030AG62.1% 39.5 / 43.6 71F410F (TTT→TTCprccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base A (13/14);  new base G (32/12);  total (45/26)
Fisher's exact test for biased strand distribution p-value = 4.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.07e-01

TCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCC  >  minE/1217435‑1217547
                                                                    |                                            
tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1696575/69‑1 (MQ=255)
tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:722899/69‑1 (MQ=255)
tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1105453/69‑1 (MQ=255)
tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1425621/69‑1 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2555550/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2822046/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1621263/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1602530/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2661281/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1528188/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:883503/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1788703/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1440988/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:772065/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1932229/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3337067/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3341397/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1210925/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3381299/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2535444/1‑69 (MQ=255)
            tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2190429/1‑69 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3081916/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:2813141/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:2362679/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:900085/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3003160/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3073639/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1051372/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3090769/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3116489/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3260708/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3393470/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3411079/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:365180/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:410175/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1945540/55‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1358199/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1642763/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1390752/1‑55 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:2008149/1‑68 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:3198642/1‑68 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:735155/1‑68 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:10806/1‑68 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:177582/1‑68 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  <  1:2008486/67‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  <  1:2220293/67‑1 (MQ=255)
              ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  >  1:1616343/1‑67 (MQ=255)
                  cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:3280361/68‑1 (MQ=255)
                  cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:146929/68‑1 (MQ=255)
                  cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:36295/68‑1 (MQ=255)
                  cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:651067/68‑1 (MQ=255)
                  cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:1105411/68‑1 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGACCAGACCGGAGAGCGa                             >  1:467881/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACa                                          >  1:2335612/1‑49 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAg                                         >  1:2570805/1‑50 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:2090072/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:3298161/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1027734/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1062778/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1459989/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1084814/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:447620/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:2955909/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:557397/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:634805/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:196345/1‑62 (MQ=255)
                        ttGCCGTTGTTATCGCGGACCAGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1570744/1‑62 (MQ=255)
                           ccGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1807803/42‑1 (MQ=255)
                            cGTTGTTATCGCTGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTc                    <  1:2233749/67‑1 (MQ=255)
                            cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTc                    <  1:1390802/67‑1 (MQ=255)
                            cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTc                    <  1:2921053/67‑1 (MQ=255)
                            cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTc                    <  1:2643548/67‑1 (MQ=255)
                            cGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTc                    <  1:2541120/67‑1 (MQ=255)
                              ttgttATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGtt                <  1:2581977/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:2008691/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:722528/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1090129/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1395845/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:3083885/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:3079815/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:2880322/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1718920/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:2712975/69‑1 (MQ=255)
                                   aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1784674/69‑1 (MQ=255)
                                    tCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:3294870/68‑1 (MQ=255)
                                            cTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGcc  <  1:760826/69‑1 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:916520/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3341627/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3236948/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3197867/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3108223/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:224390/1‑52 (MQ=255)
                                                    tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:1314150/1‑52 (MQ=255)
                                                                gAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:2573026/40‑1 (MQ=255)
                                                                    |                                            
TCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCC  >  minE/1217435‑1217547

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: