Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,217,496 C→A 44.7% A413S (GCG→TCG) ‡ 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,4960CA44.7% 42.9 / 93.0 96A413S (GCG→TCG) ‡prccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base A (32/11);  total (57/39)
Fisher's exact test for biased strand distribution p-value = 1.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

TCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCC  >  minE/1217431‑1217547
                                                                 |                                                   
tCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAgg                                                  >  1:2368468/1‑69 (MQ=255)
tCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAgg                                                  >  1:864169/1‑69 (MQ=255)
    tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGaaaa                                              <  1:2625005/69‑3 (MQ=255)
    tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:722899/69‑1 (MQ=255)
    tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1425621/69‑1 (MQ=255)
    tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1105453/69‑1 (MQ=255)
    tCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1696575/69‑1 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1788703/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1602530/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1621263/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2822046/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:772065/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2661281/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1528188/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1932229/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2555550/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2535444/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:2190429/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3381299/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3341397/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1440988/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:3337067/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:1210925/1‑69 (MQ=255)
                tCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGgag                                  >  1:883503/1‑69 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:365180/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3393470/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:2362679/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3411079/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3260708/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1051372/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:2813141/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3116489/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:410175/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:3090769/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3003160/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3081916/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:3073639/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1390752/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1358199/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              <  1:1945540/55‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:900085/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATa                                              >  1:1642763/1‑55 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:2008149/1‑68 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:10806/1‑68 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:3198642/1‑68 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:735155/1‑68 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgaga                                 >  1:177582/1‑68 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  <  1:2220293/67‑1 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  >  1:1616343/1‑67 (MQ=255)
                  ttCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGgag                                  <  1:2008486/67‑1 (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:1105411/68‑1 (MQ=255)
                      cGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGa                             <  1:651067/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:2955909/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:3298161/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1084814/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1062778/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1459989/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:447620/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:2090072/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:557397/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:634805/1‑62 (MQ=255)
                            ttGCCGTTGTTATCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGa                             >  1:1027734/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:1090129/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:2712975/69‑1 (MQ=255)
                                       aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1395845/69‑1 (MQ=255)
                                       aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:722528/69‑1 (MQ=255)
                                       aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:2008691/69‑1 (MQ=255)
                                       aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1784674/69‑1 (MQ=255)
                                       aTCGCGGACCTGAACAATCGGGCCGGATGGAATGAACAGACCGGAGAGCGATACGGCTTCAGTCAACGc           <  1:1718920/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:1314150/1‑52 (MQ=255)
                                                        tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:224390/1‑52 (MQ=255)
                                                        tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3108223/1‑52 (MQ=255)
                                                        tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3197867/1‑52 (MQ=255)
                                                        tGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGc           >  1:3236948/1‑52 (MQ=255)
                                                                 |                                                   
TCCGTCGGTATCGCTATCTTCACGAACCTTGCCGTTGTTATCGCGGACCTGAACAATGGGACCCGCAGGAATAAACAGACCGGAGAGCGATACGGCTTCAGTTAACGCCCCACCGCC  >  minE/1217431‑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: