Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 405,082 A→G 11.9% I236I (ATT→ATC mrdA ← transpeptidase involved in peptidoglycan synthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE405,0820AG11.9% 105.1 / 3.5 58I236I (ATT→ATCmrdAtranspeptidase involved in peptidoglycan synthesis
Reads supporting (aligned to +/- strand):  ref base A (21/30);  new base G (3/4);  total (24/34)
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.27e-01

TGAGATCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACATCT  >  minE/405021‑405145
                                                             |                                                               
tGAGATCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAacac                                                          >  1:2549629/1‑69 (MQ=255)
     tCCAGCTTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:1085537/69‑1 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGccaa                                                     <  1:2715859/69‑3 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:1198741/69‑1 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:1531080/69‑1 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:1454525/69‑1 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:562947/69‑1 (MQ=255)
     tCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCa                                                     <  1:1745190/69‑1 (MQ=255)
           gTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGGATAACACGCCCACGgtt                                                >  1:631066/1‑68 (MQ=255)
           gTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGGATAACACGCCCACGgtt                                                >  1:1677129/1‑68 (MQ=255)
           gTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGGATAACACGCCCACGgtt                                                >  1:2897255/1‑68 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGgttgttt                                            <  1:596580/64‑2 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:3192724/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            >  1:2094713/1‑64 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:2163206/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:2288771/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:3370327/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:849967/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            >  1:1298884/1‑64 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:1189251/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            >  1:1177607/1‑64 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:1036159/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCCAATAACACGCCCACGGTTGTTa                                            <  1:3031150/64‑1 (MQ=255)
                   aaTATCGTGTCCGGCTTGCGGGGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:1326510/64‑1 (MQ=255)
                      atCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAc                                    <  1:398877/69‑1 (MQ=255)
                      atCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAc                                    <  1:368954/69‑1 (MQ=255)
                      atCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAc                                    <  1:371572/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAATTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:521812/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:3206307/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:662439/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:578613/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:658391/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAcc                                   <  1:1363169/69‑1 (MQ=255)
                       tCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAc                                    <  1:458853/68‑1 (MQ=255)
                         gtgtCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:3188415/58‑1 (MQ=255)
                         gtgtCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAAc                                    <  1:375501/66‑1 (MQ=255)
                                gCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTa                                            <  1:1207088/51‑1 (MQ=255)
                                          gtACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGAc                >  1:2842506/1‑69 (MQ=255)
                                             cttcttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCg              >  1:1020009/1‑68 (MQ=255)
                                             cttcttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCg              >  1:2818732/1‑68 (MQ=255)
                                             cttcttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCg              >  1:1602781/1‑68 (MQ=255)
                                             ctgcttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCg              >  1:1465087/4‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCaa                                     >  1:509772/1‑42 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:3369255/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:1010484/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:1048890/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:1090364/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:1206875/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:1641340/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:180950/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:2303891/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:2635348/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:2744662/1‑68 (MQ=255)
                                                cttTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTgc           >  1:2795868/1‑68 (MQ=255)
                                                        tGGCGGATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACGTCt  <  1:2448613/69‑1 (MQ=255)
                                                        tGGCGGATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACGTCt  <  1:686059/69‑1 (MQ=255)
                                                        tGGCGGATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACGTCt  <  1:3001841/69‑1 (MQ=255)
                                                        tGGCGGATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACGTCt  <  1:937109/69‑1 (MQ=255)
                                                             |                                                               
TGAGATCCAGCGTCAGGTAAATATCGTGTCCGGCTTGCGGTGGTACTTCTTTTAACTGGCGAATAACACGCCCACGGTTGTTAACTTCAACCTCTTCATAACCGGTCTGACCGTGCAGCACATCT  >  minE/405021‑405145

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: