Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,216,822 C→A 5.6% P637P (CCG→CCT 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,216,8220CA 5.6% 300.8 / 2.5 88P637P (CCG→CCTprccarboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base C (54/29);  new base A (2/3);  total (56/32)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

CACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTT  >  minE/1216753‑1216892
                                                                     |                                                                      
cACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCg                                                                       >  1:2360117/1‑71 (MQ=255)
  ccGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGt                                                                     >  1:113797/1‑71 (MQ=255)
   cGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGtt                                                                    >  1:1363895/1‑71 (MQ=255)
    gTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGttt                                                                   <  1:2650658/71‑1 (MQ=255)
    gTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGttt                                                                   <  1:1155949/71‑1 (MQ=255)
       tcATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTAcc                                                                >  1:2226716/1‑71 (MQ=255)
       tcATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTAc                                                                 >  1:970846/1‑70 (MQ=255)
       tcATACAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTAcc                                                                >  1:4064465/1‑71 (MQ=255)
          tCCAGATAAGGATCCGGCTCCAGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTc                                                             <  1:1604929/71‑1 (MQ=255)
           ccAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTcg                                                            >  1:78710/1‑71 (MQ=255)
           ccAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCAGGTTTACCTTc                                                             >  1:3847565/1‑70 (MQ=255)
                tAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGttt                                                                   <  1:946584/59‑1 (MQ=255)
                tAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGtt                                                        <  1:259808/70‑1 (MQ=255)
                 aaGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGttt                                                       >  1:2793873/1‑70 (MQ=255)
                 aaGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGttt                                                       >  1:1694726/1‑70 (MQ=255)
                  aGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTa                                                      >  1:2486763/1‑70 (MQ=255)
                   ggATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTa                                                      >  1:892913/1‑69 (MQ=255)
                          gCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTcgcg                                                          <  1:280670/58‑1 (MQ=255)
                              cTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTAc                                                                 >  1:3895548/1‑47 (MQ=255)
                              cTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGcgttcgtt                                           <  1:1564455/69‑1 (MQ=255)
                              cTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCgttcgttc                                          <  1:866098/70‑1 (MQ=255)
                              cTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCa                                         <  1:1246512/71‑1 (MQ=255)
                              cTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCa                                         <  1:2592601/71‑1 (MQ=255)
                               tGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGc                                                  >  1:2121309/1‑61 (MQ=255)
                                 gTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGcgttcgtt                                           >  1:3134516/1‑66 (MQ=255)
                                 gTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCaaa                                       <  1:2199390/70‑1 (MQ=255)
                                 gTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAc                                      >  1:3362022/1‑71 (MQ=255)
                                 gTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAc                                      >  1:3968024/1‑71 (MQ=255)
                                  tAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAcg                                     <  1:580919/71‑1 (MQ=255)
                                   aaTCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAcgc                                    >  1:3930742/1‑71 (MQ=255)
                                   aaTCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAcgc                                    >  1:1110624/1‑71 (MQ=255)
                                   aaTCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAAcgc                                    >  1:1115275/1‑71 (MQ=255)
                                      cTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTaa                                                     <  1:3974677/51‑1 (MQ=255)
                                       tttCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTa                                                      >  1:1068552/1‑49 (MQ=255)
                                       tttCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCaaa                                       >  1:27872/1‑64 (MQ=255)
                                       tttCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACgcgc                                  >  1:1873397/1‑69 (MQ=255)
                                         tCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTa                                                                  >  1:3218124/1‑35 (MQ=255)
                                         tCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGcgttc                                              >  1:3238108/1‑55 (MQ=255)
                                         tCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACgcgc                                  <  1:1947579/67‑1 (MQ=255)
                                         tCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCa                                <  1:3052374/69‑1 (MQ=255)
                                           ggTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGAc                             <  1:1646680/70‑1 (MQ=255)
                                              aGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGt                         >  1:411659/1‑71 (MQ=255)
                                               gATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACgcgc                                  <  1:2444940/61‑1 (MQ=255)
                                               gATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTc                        >  1:4132107/1‑71 (MQ=255)
                                                 tcatcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTa                                                      <  1:4071102/39‑1 (MQ=255)
                                                 tcatcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCg                       >  1:587311/1‑70 (MQ=255)
                                                 tcatcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGc                      <  1:4121694/71‑1 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTa                                                      >  1:1251299/1‑38 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCaaa                                       >  1:153709/1‑53 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAgacgcgacgc                      >  1:3176532/1‑70 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGca                     >  1:1649308/1‑71 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGca                     >  1:1618329/1‑71 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGca                     >  1:323911/1‑71 (MQ=255)
                                                  catcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGc                      >  1:1574883/1‑70 (MQ=255)
                                                   atcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTaaa                                                    >  1:4123272/1‑39 (MQ=255)
                                                   atcCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGca                     <  1:2659215/70‑1 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAg                                                   >  1:2891457/1‑38 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGcgttc                                              >  1:73078/1‑43 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCa                                         <  1:3548886/48‑1 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGcatca                  >  1:1457228/1‑71 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGcatca                  >  1:2076528/1‑71 (MQ=255)
                                                     ccAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGc                      <  1:368664/67‑1 (MQ=255)
                                                     ccAGTTTCTTCAACTCAGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGAcgc                           <  1:124933/62‑1 (MQ=255)
                                                      cAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGa                              >  1:4117657/1‑58 (MQ=255)
                                                      cAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGa                              >  1:1105112/1‑58 (MQ=255)
                                                      cAGTTTATTCAACTCCGGTTTACCTTCGCGTTTAAAGCgttcgttc                                          >  1:828541/1‑46 (MQ=255)
                                                        gTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCt               >  1:824028/1‑71 (MQ=255)
                                                         tttCTTCAACTCCGGTTTACCTTCGCGTTTAAAGcgtt                                               <  1:3339781/38‑1 (MQ=255)
                                                         tttCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTaaaa                                       >  1:1056649/1‑46 (MQ=255)
                                                         tttCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCtt              <  1:309450/71‑1 (MQ=255)
                                                         tttCTTCAACTCAGGTTTACCTTCGCGTTTAAAGCGTTCGTTCaaa                                       >  1:625127/1‑46 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCtt              >  1:904124/1‑69 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTc             <  1:687988/70‑1 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCa            >  1:1001514/1‑71 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCa            >  1:353638/1‑71 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCa            >  1:128120/1‑71 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCa            >  1:2564361/1‑71 (MQ=255)
                                                           tcttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCa            >  1:2035752/1‑71 (MQ=255)
                                                             ttcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGa                              >  1:3951845/1‑51 (MQ=255)
                                                              tcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCatta         >  1:3474218/1‑71 (MQ=255)
                                                              tcAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCatt          >  1:1704864/1‑70 (MQ=255)
                                                               cAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCattat        >  1:3926712/1‑71 (MQ=255)
                                                               cAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCatta         <  1:2313273/70‑1 (MQ=255)
                                                               cAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCat           <  1:3084279/68‑1 (MQ=255)
                                                               cAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCat           <  1:3856097/68‑1 (MQ=255)
                                                               cAACTCAGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGACAGACGCGTCGCATCATCTTCat           <  1:1840076/68‑1 (MQ=255)
                                                               cAACTCAGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGACAGACGCGTCGCATCATCTTCat           <  1:1835352/68‑1 (MQ=255)
                                                                   tCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGcatca                  <  1:995824/57‑1 (MQ=255)
                                                                    ccGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGcatcat                 >  1:1955721/1‑57 (MQ=255)
                                                                    ccGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATtctc    >  1:2685296/1‑70 (MQ=255)
                                                                    ccGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTctct   >  1:1397616/1‑71 (MQ=255)
                                                                    ccGGTTTACCTTCGCGTTTAAAGAGTTCGTTCAAACGCGCCAg                               >  1:662717/1‑43 (MQ=255)
                                                                     cGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTctct   <  1:3109824/70‑1 (MQ=255)
                                                                     cGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCtt  <  1:1366477/71‑1 (MQ=255)
                                                                     |                                                                      
CACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTTCTTCAACTCCGGTTTACCTTCGCGTTTAAAGCGTTCGTTCAAACGCGCCAGACGCGTCGCATCATCTTCATTATTCTCTT  >  minE/1216753‑1216892

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: